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Presentación 


Este libro, PROBLEMAS RESUELTOS DE INGENIERÍA QUÍMICA, junto unos materiales 
docentes editados en el 2015 bajo el título IncemErRÍA Química (QUO916), está 
concebido para que el estudiantado de la Universitat Jaume I que cursa la asignatura 
con el mismo nombre en el Grado en Química durante el segundo semestre del 
segundo curso de la titulación (6 ECTS). 

Este proyecto de libro de problemas resueltos surge de la necesidad de tener un 
manual de referencia que se ajuste al programa de la asignatura, y que permita al 
alumnado optimizar su tiempo de estudio y obtener un rendimiento satisfactorio 
en esta asignatura. 

PROBLEMAS RESUELTOS DE INGENIERÍA QUÍMICA está estructurado en cuatro bloques, 
cada uno de ellos se corresponde con un tema de la asignatura, donde se plantea y 
resuelve la colección de problemas: 


+ Bloque I. Balances de materia en sistemas sin reacción química 

+ Bloque II. Balances de energía total en sistemas sin reacción química 

+ Bloque III. Balances de materia en sistemas con reacción química. Reactores 
ideales 

+ Bloque IV. Balances de energía en sistemas con reacción química. Reactores 
ideales 


Para cada enunciado, se presenta el diagrama del proceso, un resumen en el 
que se recogen y esquematizan los datos del enunciado, se plantean los balances 
de propiedad necesarios para la resolución del mismo, dando al alumnado las 
pautas necesarias para su resolución, y se incluye la solución gráfica y/o numérica 
correspondiente, con el fin de que el estudiantado pueda evaluar la resolución 
realizada. Al final del libro, en los «Anexos», se puede encontrar un formulario 
correspondiente a cada bloque. 

Este libro, al igual que los materiales anteriormente ciados, es el resultado de 
la experiencia de los autores en el conjunto de las operaciones unitarias y de los 
procesos químicos. Pretende ser una herramienta de estudio del alumnado que 
le permita conocer, analizar, modelizar y calcular sistemas sin y con reacción 


química. 
V. SANZ SOLANA, N. VICENTE AGUT 
Castellón 2023 
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BLOQUE | 


Balances de materia 
en sistemas sin reacción 
quimica 


Problema 3.1 


En una torre de destilación se introducen 100 kg/h de una mezcla que contiene 35 % de etanol y 
65 % de agua. El destilado contiene un 85 % de etanol y el residuo un 5 % de etanol. Calcular: 

(a) kg de destilado por cada kg de alimento introducido. 

(b) kg de residuo por cada kg de alimento introducido. 

(c) kg de destilado por cada kg de residuo obtenido. 


Datos y notas: 


El sistema funciona en régimen estacionario. 
Todos los porcentajes son en peso. 


DIAGRAMA DEL PROCESO 
El sistema puede representarse esquemáticamente mediante el siguiente diagrama: 


Corriente 2: destilado 


85 % etanol 
| 15 % agua 
Corriente 1: alimento 


my = 100 kgh ————— Destilador 


35 % etanol 
| Corriente 3: residuo 
5 % etanol 


65 % agua 
95 % agua 


El sistema funciona en régimen estacionario y no hay reacciones químicas. 

El proceso tiene una etapa (destilador) y dos componentes, el etanol (1) y el agua (2). 

Los balances de materia se van a realizar en unidades másicas (kg/h) dado que los datos de las 
corrientes están en peso. 


Según el enunciado, el caudal de alimento es: m,:=100 e 
r 


Las fracciones másicas de los componentes (kg componente/kg totales) de cada una de las 


Componente 2 
agua 


corrientes del sistema son (*): 
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(*) Las fracciones másicas siempre estarán referenciadas como ww», donde f hace referencia al 
componente e / a la corriente. 


SOLUCIÓN 


Planteando los balances de materia total y de componente en el destilador, se obtiene el siguiente 
sistema de ecuaciones: 


BM total: M, =M,+My (1) 


BM etanol: — m,-w,,=M,-*W,,+M3*W13 (2) 


Es un sistema de dos ecuaciones con dos incógnitas, que puede resolverse fácilmente: 
M, = M; — Mg 
sustituyendo (1) en (2): MW, = (m, —- my) 12 + Mz*W313 


m;+ (13, > W13) =M3* (w33 al W19) 


Ñ My * (011 —Wy k 
el caudal de residuo es: e A EA 
(w33 a w1>) 8 
kg 
y el caudal de destilado: m,:=m, —m¿=0.01 = 
Ss 
Los resultados se piden como las siguientes razones: 
Ma a , 
(a) —=0.375 kg destilado/kg alimento 
my 
m3 . : 
(b) —=0.625 kg residuo/kg alimento 
UL 
TM j j 
(c) —= uE kg destilado/kg residuo 
m3 
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Problema 3.2 


Un proceso de destilación que funciona en régimen estacionario recibe una alimentación continua de 
100 kmol/h de una mezcla que contiene un 35 % en moles de tolueno y un 65 % en moles de 
benceno. Del proceso se obtienen dos corrientes, una con una fracción molar de 0.99 de benceno y 
otra que contiene el 5 % del benceno de la alimentación. Hallar los caudales molares de ambas 
corrientes y la composición de la corriente rica en tolueno. Determinar el caudal másico y la 
composición en porcentaje en peso de la corriente de alimento. 


DIAGRAMA DEL PROCESO 


El sistema puede representarse esquemáticamente mediante el siguiente diagrama: 


Corriente 2: destilado 
95 % benceno 
1 % tolueno 


Corriente 1: alimento 


F, = 100 kmol/h Destilador 
35 % tolueno 
65 % benceno 
Corriente 3: residuo 
5 % benceno del alimento 


El sistema funciona en régimen estacionario y no hay reacciones químicas. 

El proceso tiene una etapa (destilador) y dos componentes, el benceno (1) y tolueno (2). 

Los balances de materia se van a realizar en unidades molares (kmol/h) dado que los datos de las 
corrientes están en kmol y se indican las fracciones molares de las corrientes. 


kmol 


Según el enunciado, el caudal de alimento es: F,:=100 a 
Tr 


Y las fracciones molares de los componentes (kmol componente/kmol totales) de cada una de las 
corrientes del sistema son: 


Componente 1 Componente 2 
(CsHs, benceno) [C7Hs, tolueno 


A moss | emos 7 
III IEC NTE 


La condición del problema es: F,+x,13=0.05-F,-2;1 
Para cualquier corriente: > 1 =1 
j 
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SOLUCIÓN 


Planteando los balances de materia en estado estacionario al destilador y teniendo en cuenta la 
condición del problema: 


BM total: F,=F,+F, (1) 
BM benceno PF, 0 =FysD y EFD (2) 
Condición del problema: F¿+21350.05-F,+2;, (3) 


que es un sistema de tres ecuaciones con tres incógnitas que se resuelve fácilmente: 


Sustituyendo (3) en (2) F¡+21¡=F3-2%12+0.05-F, +21, (4) 
. F,+2,,—0.05*F,.x kmol 
despejando F2 y F3 PE 
Li hr 
kmol 
Fi=F,=Fo=376 == 
hr 
La composición de la corriente rica en tolueno es: 
0.05-F,-x 
21322 = 0.086 


3 


L73:=1—2,3=0.914 


Para el cálculo del caudal másico y la composición del alimento es necesario conocer los pesos 
moleculares de los componentes: 


k 
Peso molecular benceno: PM, :=78 zd 
kmol 
k 
Peso molecular tolueno: PM,:=92 Y 
kmol 
El caudal del alimento será: 
kg ay Kg 
my := (F, +21, +PMy) + (E, +271+PM)) =2.303 m,=(8.29.10*) 2 
Ss hr 
y su composición másica: 
F; Li1 -PM; A 
Wy1 == 0.612 kg benceno/kg totales corriente 1 
my 
F; + La PM) . 
Wo := ———— = 0.388 kg tolueno/kg totales corriente 1 
my 
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Problema 3.3 


Un evaporador funciona en régimen estacionario y se alimenta continuamente con 25 toneladas/h 
de una disolución que contiene el 80 % de agua, el 10 % de NaOH y el 10 % de NaCl. Durante la 
evaporación, el NaCl precipita en forma cristalina y se sedimenta. La solución concentrada que 
abandona el evaporador contiene un 48 % de agua, un 50 % de NaOH y un 2 % de NaCl. Calcular: 

(a) kg de agua evaporados por hora. 

(b) kg de NaCl precipitados por hora. 

(c) kg de solución concentrada producidos por hora. 


DIAGRAMA DEL PROCESO 
El sistema se representa esquemáticamente mediante el siguiente diagrama: 


Corriente 2: 100 % H,0 


Corriente 1: alimento | Corriente 4: solución concentrada 
my = 25 t/h PA ; 48 % agua 
80 % agua Eva porador 50 % NaOH 
10 % NaOH 2 % NaCl 
10 % Nacl 


Corriente 3: precipitado 
100 % NaCl 


El sistema funciona en régimen estacionario y no hay reacciones químicas. 
El proceso tiene una etapa (evaporador) y tres componentes: H20 (1), NaOH (2) y NaCl (3). 
Los balances de materia se van a realizar en unidades másicas (kg/h). 


, , z tonne 
Según el enunciado, el caudal de alimento es: m,:=25 


= 25000 42 
r hr 
Y las fracciones másicas de los componentes (kg componente/kg totales) de cada una de las 
corrientes del sistema son: 


SOLUCIÓN 
Aplicando los balances al sistema: 
BM total M¡ =M,+Mz+Ma (1) 
BM H20 M,¿*W¡ =M2+ 0,2 + Mz*W,3 + My¿*W4 (2) 
BM NaOH My * Wo] = M * Way + Mg * Waz + My * Way (3) 
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El NaOH es un componente enlace, por lo que la ecuación (3) puede resolverse directamente: 


(5.10?) $9) 


Wo hr 


Mj¿* Wo] 
() me 


Nos queda un sistema de dos ecuaciones con dos incognitas, que se resuelve por cualquier método: 
Despejando m, de la ecuación (1): M¿=M, —Mz¿— My 
y sustituyendo en la ecuación (2): MW, = (m, —m¿—My) -W¡ + Mz*W¡3 + M¿*Wi4 

My *W11 =M, *W19 — Mgz*Wj) — My *W19 + Mg* Wiz + My Wa 
Despejando my: 


My *W¡¡ My *Wi¡9 + My¿* Wi) My w k 
(b) ma:= it 1 ta 40 01 4 E =dom) Kg 
W¡3—W1 hr 


Despejando m,: 


k 
ay K9 
(a) m,:=m, —mz=m,=(1.76-10 y paa 
hr 
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Problema 3.4 


Necesitamos secar 1000 kg/h de un sólido que posee una humedad del 40 %, hasta obtener un 
producto con una humedad del 20 %. ¿Qué cantidad de agua debemos eliminar? 


DIAGRAMA DEL PROCESO 


El sistema puede representarse esquemáticamente mediante el siguiente diagrama: 


Corriente V: agua evaporada 


| Corriente P: producto 
Secadero -——> 20 % humedad 


Corriente A: alimento 
ma = 1000 kg/h 
40 % humedad 


El sistema funciona en régimen estacionario y no hay reacciones químicas. 
El proceso tiene una etapa (secador) y dos componentes: agua (1) y sólido (2). 


Según el enunciado, el caudal de alimento es: m,a:=1000 2 
Tr 


Y las fracciones másicas de los componentes de cada una de las corrientes del sistema son: 


Componente i Componente 2 
(agua) sólido 


SOLUCIÓN 
La suma de fraciones másicas de cada corriente es 1. 


Aplicando balances al sistema: 


BM total: M¿=My+Mp (1) 
BM sólido: My *Wy¿ = My * Way + Mp Wap (2) 


Usando el sólido como componente enlace: 


Mpa*Wza kg 
mp == = 750 
Wap hr 
el caudal de agua evaporada será: 
kg 
My := My —Mp=250 — 
hr 
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Problema 3.5 


En un secadero continuo, consistente en dos hornos instalados en serie, se deseca un material 
húmedo que contiene un 10 % de sólido y un 90 % de agua. El producto a la salida del primer 
horno contiene un 20 % de agua y, en el segundo horno, se reduce su contenido en humedad a un 
2 %. Si la capacidad del secadero es de 1000 kg de sólido húmedo/día, calcular el peso de producto 
que diariamente sale de cada uno de los hornos. 


DIAGRAMA DEL PROCESO 
El sistema puede representarse esquemáticamente mediante el siguiente diagrama: 


Corriente 2: agua evaporada Corriente 4: agua evaporada 


Corriente 1: alimento | | Corriente 5: producto 
ATA] Horno 1 — Horno 2 H—> 2 % humedad 
m, = 1000 kg mat. húmedo/día o ei 


Corriente 3: 20 % agua 
80 % sólido 


El sistema funciona en régimen estacionario y no hay reacciones químicas. 
El proceso tiene dos etapas (Horno 1 y Horno 2) y dos componentes: agua (1) y sólido (2). 
Según el enunciado, el caudal de alimento es: m, :=1000 E 

ay 


Para plantear los balances existen tres entornos posibles: Horno 1, Horno 2 y el sistema completo. 
En cada entorno se pueden plantear dos ecuaciones linealmente independientes. 

Se deben elegir los dos entornos que convengan, normalmente empezando por la corriente mejor 
definida. 


Las fracciones másicas de los componentes de cada una de las corrientes del sistema son: 


SOLUCIÓN 


Balance en el Horno 1: 


BM total: M¡ =M,+My (1) 

BM sólido: M¡+W,, =M)- 013 + Mg + W33 (2) 
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Resolviendo el sistema: 
kg 


m3 = 
Wiz day 
day 
Balance en el Horno 2: 
BM total: Mz¿=M,¿+Ms (3) 
BM sólido: Mgz*W;¡3 = My *W¿+ Mg*Wj5 (4) 


Resolviendo el sistema: 


Mai: w k 
ma =102.041 3 a 
ay 
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Problema 3.6 


El ácido residual de un proceso de nitración contiene 23 % de ácido nítrico, 57 % de ácido sulfúrico 
y 20 % de agua en peso. Este ácido debe concentrarse, hasta que contenga un 27 % de HNO» y un 
60 % de H2SO4, mediante la adición de ácido sulfúrico concentrado con un contenido del 93 % de 
de H2S0» y de ácido nítrico concentrado con un 90 % de HNOz. 


Calcular los caudales másicos de ácido sulfúrico concentrado, ácido nítrico concentrado y ácido 
residual que deben combinarse para obtener 1000 kg/h de la mezcla deseada. El sistema funciona 
en régimen estacionario. 


DIAGRAMA DEL PROCESO 
El sistema puede representarse esquemáticamente mediante el siguiente diagrama: 
Corriente 1: 23 % HNOz 

| 57 % H>SO, 


20 % H,O 
Corriente 2: Corriente 4: 27 % HNO; 


93 % H,SO, 60 % H,SO, 
7% H,O Mezclador 13 % HO 


Corriente 3: 90 % HNO, 
10 % H,0 


El sistema funciona en régimen estacionario y no hay reacciones químicas. 
El proceso tiene una etapa (mezclador) y tres componentes: HNOz (1), H2SOs (2) y H20 (3). 


A Ñ k 
Según el enunciado, el caudal de producto deseado es: m,:=1000 se 
> 


Las fracciones másicas de los componentes de cada una de las corrientes del sistema son: 
Componente i Componente 2 Componente 3 
(HNO3) H2S04 (H20 


SOLUCIÓN 

Aplicando balances al mezclador: 

BM total: M¡ +M,+Mg= Ma (1) 
BM HNO:: MW +FMy+W19 + Mz" W¡3= My" Wi (2) 
BM H2SO4: M; * Wa] + Mo" Wa7 + Mg * Wyz = My * Way (3) 


Se obtienen un sistema de tres ecuaciones con tres incógnitas, que puede resolverse fácilmente: 


kg kg kg 

m, =416.9 — mo, =389.6 — m3=193.5 — 

hr hr hr 
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Problema 3.7 


Se desea recuperar la acetona que forma parte de una corriente gaseosa junto con aire. Las 
operaciones y los datos conocidos se detallan en el siguiente diagrama de flujo. 


Calcular los valores de los caudales másicos de todas las corrientes. 


DIAGRAMA DEL PROCESO 


Corriente 3: 
Aire seco 


Corriente 2: ma, = 1400 kg/h 


| | Condensador 


| Absorbedor 


Corriente 5 


Corriente 1: |] [| Rectificador 


Aire + acetona 


—— 


15 % en moles de acetona Corriente 4 


Corriente 6: 
Destilado 
99 % acetona 
1 % agua 
(en peso) 


Corriente 7: 
residuo 

5 % acetona 
95 % agua 
(en peso) 


El sistema funciona en régimen estacionario y no hay reacciones químicas. 
El proceso tiene tres etapas (absorbedor, rectificador y condensador) y tres componentes: 


acetona (1), agua (2) y aire (3). 
Según el enunciado, los caudales de las corrientes 2 y 3 son: 


kg 


my:=1200 — mg:=1400 
h 


Tr 


kg 


El problema se resolverá en unidades másicas, por lo que es necesario recalcular los datos dados en 


unidades molares: 


x2,1:=0.15 31 :=0.85 PM, :=58.08 PM; :=28.9 
211 PM, 31 * PM; 
w;]3= =0.26 Waz] i= =0.74 
211 * PM, +23, PM; 211 * PM, +23, PM; 
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Las fracciones másicas de los componentes de cada una de las corrientes del sistema son: 


Componente i Componente 2 Componente 3 
Acetona Agua Aire' 
E 
- 


SOLUCIÓN 


El sistema tiene, aparentemente, tres etapas; sin embargo, en el condensador solo se produce un 
cambio de energía, sin variación en la composición de las corrientes 5 y 6, por lo que se puede 
tratar como si tuviera 2 etapas. En cada entorno, se puede plantear tres ecuaciones 
independientes. 


Balance en el absorbedor: 


En este entorno es donde más información se tiene, por lo que hay que intentar empezar el 
problema por el absorbedor. Los balances son: 


BM total: M, + My, = Mz + My (1) 
BM aire: My + W3z] =Mg*Wsg (2) 
BM agua: My * Way = My * Wo (3) 


Además, en todas las corrientes siempre debe cumplirse que la suma de las fracciones de los 
componentes debe ser 1. En este problema, tendremos que en la corriente 4. 


Corriente 4: W¡¿+FWa+Ww34=1 (4) 


Con lo que tenemos 4 ecuaciones con cuatro incognitas, que resolvemos por cualquier método. 


k k 
m,:=1895.73 2 m,:=1695.73 eN 
Tr r 
mMit*wW 
Wa == =0.71 w,,'=1—w,,—3,=0.29 
4 


Balance en el rectificador-rondensador: 


Los balances son: 
BM total: M,¿=M¿+mM, (5) 
BM agua: My * Way = Mg * Waz + My * Way (6) 


Que se resuelve por cualquier método. 


kg kg 

m¿=437.17 m,=1258.56 2 

hr hr 
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Problema 3.8 

En un secadero continuo, que funciona en régimen estacionario, se introduce pulpa de papel que 
contiene un 71 % de agua en peso. A la salida de este se encuentra que se ha eliminado el 60 % 
del agua inicial. Calcular: 


(a) La composición de la pulpa de papel a la salida. 
(b) El peso de agua eliminada por cada kg de pulpa húmeda introducida. 


DIAGRAMA DEL PROCESO 


El sistema puede representarse esquemáticamente mediante el siguiente diagrama: 


Corriente 2: Condición: m,+w,, = 0.6-m,*w», 


Corriente 1: 71 % agua ————>| Secadero Corriente 3 


El sistema funciona en régimen estacionario y no hay reacciones químicas. 
El proceso tiene una etapa (Secadero) y dos componentes: pulpa de papel (1) y agua (2). 


Según el enunciado, la condición que debe cumplirse es: M,* Wa =0.6-M;¡ Way] 


A k 
Se toma como base de cálculo: m,:=100 Fa 
A 


Las fracciones másicas de los componentes de cada una de las corrientes del sistema son: 


Componente 1 Componente 2 
¡pulpa de papel) agua 


AT esoo | eso 


SOLUCIÓN 


Planteando los balances y la condición del problema: 


BM total: mM, =M2+My (1) 

BM pulpa de papel: M¡+W, =M3*W33 (2) 

Corriente 3: 13 +W3=1 (3) 

Condición problema: ma, 39 =0.6-m; + Wa] (4) 
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Con lo que tenemos 4 ecuaciones con cuatro incognitas, que resolvemos por cualquier método. 


mm =49.0 42 
hr 
k 

ma:=57.4 49 
hr 


(a) La composición de la pulpa de papel a la salida es: 


Wwoz:=1—w:13=0.49 


(b) La razón de kg agua eliminada por kg de pulpa húmeda introducida es: 


; ma 
Razón := —=0.43 
m; 
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Problema 3.9 


En una refinería de petróleo interesa rebajar el contenido en iso-pentano de una corriente que 
contiene inicialmente el 80 % de n-pentano yel 20 % de iso-pentano (en peso) hasta el 10 % de 
iso-pentano. Para ello se hace pasar a través de una unidad en la que se elimina por completo el 
iso-pentano y la corriente obtenida se mezcla con parte de la corriente de alimentación para 
conseguir la especificación deseada, de acuerdo con el siguiente esquema: 


Corriente W: Isopentano 


Columna de 


isopentano 


¿Qué cantidad de corriente de alimento se debe mezclar con la corriente producto sin pasar a través 
de la columna en la que se elimina el iso-pentano? 


Nota. Todos los porcentajes son el peso. 


El sistema funciona en régimen estacionario y no hay reacciones químicas. 
El proceso tiene tres etapas (divisor de la corriente F, columna, mezclador de las corrientes C y B) y 
dos componentes: n-pentano (1) y iso-pentano (2). 


Se toma como base de cálculo: mp:=100 =Ñ= 


SOLUCIÓN 


Balance global: 


BM total: Mp= My + Mp (1) 

BM n-pentano: Mp*W¡p= My * W¡y + Mp*W¡p (2) 
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Resolviendo el sistema: 
My = mp— Mp 
Mp*Wip= (mp—mg) *Wiy + Mg * W¡5=Mp* Wiy — Mp* W¡y + Mp*W¡p 


Mp*W¡p= (mp— M4) *Wiy + Mg*W¡¿=Mp*W¡y + Mp+ (wip= Wiw) 
_ Mp» (w,p—W1yy) kg kg 


Mp:= =88.89 — My:=mp=mp=11.11 — 
Wi¡p— Wiyw hr hr 


Balance en la columna: 


BM total ma= My + Mao (3) 
BM n-pentano: My "Wi 4 = Myy + Wi yy + Mor Wic (4) 
Resolviendo el sistema: 
(myy+ mg) *Wj4 = My * Wi yy + Mg" Wig 


My * (W,4 = Wiw) =Mg* (w35 = W14) 


si RA EA ia mia 2 


(wc —- W14) hr hr 


Balance en el punto de unión: 
Mp=M,4+Mg (5) 


k 
Mg:=Mp—M,= 44.44 > 


La razón de kg en el bypass por kg de alimento es: 


m 
Razón:=—=0.44 kg en el bypass/kg de alimento 
Mp 
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Problema 3.10 


Necesitamos separar un componente X de una mezcla (X+Y). Para ello empleamos un 
sistema de compresión-expansión adiabática tal como se muestra en la figura. 


P 


Compresor Válvula de 
expansión 


La mezcla (X+Y) se comprime y después se deja expansionar a través de una válvula condensando 
el componente X y quedando el componente Y en estado gaseoso. Sin embargo, parte del 
componente X es arrastrado por el componente Y, que al ser un producto valioso lo volvemos a 
introducir en el sistema. Si no verificamos una purga en la corriente de recirculación se nos irá 
acumulando continuamente el producto Y dentro del sistema. 


La corriente de alimentación, que se introduce en el sistema a razón de 100 kg/h, está formada por 
un 80 % de X y un 20 % de Y. Debido a las características del compresor, nos interesa que la 
concentración de Y a la entrada del compresor sea del 40 %. La concentración de X en la corriente 
de recirculación (D) resulta ser del 30 %. Calcular: 


(a) Cantidad de producto X obtenido por hora. 
(b) Caudal másico (kg/h) de la corriente de purga. 
(c) Cantidad de gas que ha de comprimir el compresor (kg/h). 


El sistema funciona en régimen estacionario y no hay reacciones químicas. 
El proceso tiene tres etapas (mezclador de las corrientes R y F, separador y purga) y dos 
componentes: X y Y. 


Caudal másico corriente alimento (F): mp:=100 


Las fracciones másicas de los componentes de cada una de las corrientes del sistema son: 


| Corriente_ | Componente X_ | Componente Y | 
O TJ oso] moyp=02o 
A oo | eraz0n0 
Roo | ooo 


INCAS INT "SCETI INSI 
030 | poo 
O 0 
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SOLUCIÓN 


Balance mezclador de corrientes R y F: 


BM total 
BM de X 


resolviendo el sistema: 


(c) 


Balance en el separador: 


BM total 
BM de X 


resolviendo el sistema: 


(a) 


Balance en la purga: 


BM total 


(b) 


Ma =Mp+ Mr 


Ma *Wxy¿ =Mp* Wxp + Mp*Wxp 


(mp+ Mp) *Wx¿ = Mp*Wxp + Mp*Wxp 


Mp+* (Wwxa a Wxp) =Mp* (wxp e Wxa) 


_ Mp+ (Wx4—Wxp) —66.67 kg 


(wr = Wxa) hr 


mp: 


m,:=mp+mp=166.67 — 


Mg =Mp+ My 


Ma¿*Wxa¿ = Mp *Wxp + Mxyy * Wxy 


Mp =Ma — My 
M¿*Wxp4= (Ma —Myy) *W5 EM yy" xy 


ma (Wwxa Sl Wxp) = My * (20 xy al Wxp) 


My := (e (xa —"wxp)) ="1.43 Eg 
(Wxyy — xp) hr 
k 
Mp'= Mp — My = 95.24 se 
Mp=Mp+Mp 
k 
Mp:=Mp—Mp= 28.57 e 


La razón de kg en el bypass por kg de alimento es: 
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(2) 


(3) 
(4) 


m 
Razón :=—2=0.95 kg en el bypass/kg de alimento 


Mp 


28 
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Problema 3.11 


Un sistema continuo de evaporación-cristalización funciona en régimen estacionario según el 
esquema siguiente: 


Corriente 4: agua evaporada 


Corriente 5: 
cristales de KNO,, 
conteniendo 
4 kg H,0/100 (kg 

KNOz + H,0) 


Corriente 1: alimento Evaporador ——=: Cristalizador => 


Corriente 3 


La corriente alimento es de 10 toneladas/h, y la concentración de KNOz en la solución alimento es 
del 20 % en peso, la solución 3 está saturada (0.6 kg KNOz/kg agua), y la solución 6 tiene una 
concentración del 50 % en peso de KNOz, calcular: 


(a) kg de agua eliminados por hora en el evaporador. 
(b) Caudal másico (kg/h) de solución en las corrientes 6 y 3. 


El sistema funciona en régimen estacionario y no hay reacciones químicas. 
El proceso tiene tres etapas (Mezclador corrientes 1 y 3, Evaporador, Cristalizador) y dos 
componentes: KNOs (1) y H20 (2). 


Según el enunciado la corriente alimento es: m, := 10000 aa 
de 
Cálculo de las concentraciones según los datos del enunciado: 
] 0.6 
Corriente 3: Wi3 e Waz:=1—w,3=0.625 
Corriente 5: was := 0.04 w15:=1—w35;=0.96 


Las fracciones másicas de los componentes de cada una de las corrientes del sistema son: 
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SOLUCIÓN 
Balance global: 


BM total mM¡ =M,¿+My (1) 

BM KNOz M¡*W] =M,*W¡4+ M5 + W35 (2) 
Resolviendo el sistema 

M¿=M1 — My M, +11 = (M, —M5) 014 + M5 +35 


m.+ (w, = W14) =M5* (w35 - W14) 


az m]+ (1,1 014) kg 
a (w,5—34) hr 


(a) Los kg de agua eliminados por hora en el evaporador son: 


m,:=m, —m¿="7917 Es 
Balance en el cristalizador: 
BM total M¿= Mz + Ms (3) 
BM de X Mg *W¡¿= Mg + W¡3 + Mg -*Wj5 (4) 


Resolviendo el sistema 
(my +my5) *W]=Mg*W13 + M5" Wj5 


mg (w,6 > W13) =Mg+ (w,5 > W16) 


(b) El caudal másico (kg/h) de solución en las corrientes 6 y 3 es: 


m;+- (w —w ) k 
5 15 16 g 
my 50 7667 +2 
(w,¿—w33) hr 
kg 
Y 
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Problema 3.12 


En un tanque perfectamente agitado que contiene inicialmente 500 L de una disolución de KOH al 
80 % en peso, se introduce una corriente de agua pura con un caudal de 40 L/min, al mismo tiempo 
que se extrae una corriente con un caudal de 60 L/min. Determinar las ecuaciones que nos dan la 
composición y la masa de disolución contenida en el tanque en función del tiempo. La densidad de 
las disoluciones se puede considerar que coincide con la del agua pura (1 kg/L). 


DIAGRAMA DEL PROCESO 


El diagrama queda definido del modo siguiente: 


Q, = 40 L/min 


El sistema funciona en régimen no estacionario. 


Densidad de la disolución: 


Volumen y concentración iniciales en el tanque: 


Caudal y concentración de entrada: 


Caudal y concentración a la salida: 


p:=1000 
V¿:=500 
w:=0.8 
Qy¡:=40 
w,:=0 

Qy,:=60 
Wa =W 


Q), = 60 L/min 
kg 
a 
L Mo:=Vo¿*p=500 kg 
L k 
] my:=Q,p=40 2 
min min 
L k 
- my:=Q,p=60 
min min 


SOLUCIÓN 


El balance siempre se plantea en términos másicos, pero considerando la densidad como factor 
común, también se puede plantear en términos volumétricos. 


dM 
BM total mimo (QQ + Ge 


V 
QQ +50 


dV 


Reordenando E TÚR Ori O 


d 
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dV = 
dt 


Necesitaremos esta ecuación en 
el balance de componente 


Vicente Sanz Solana, Nuria Vicente Agut 
DOI: http://dx.doi.org/10.6035/Sapiential90 


Vv t 
Separando variables fi dV=(Q11 Qro)" J 1 di 
VA 0 


Integrando V-Vo= (Q11 Qro) *t 


Reordenando V (t):=Vo+(Qr, —Q;») et V(5-min)=0.4 m* 


Necesitaremos esta ecuación en el balance de componente 


Planteando ahora el balance de componente en estado no estacionario: 


M 
BM componente: my) my un O y, Oya pr O =0 
d(Vw 
Qr*W-Qy; w,+ ná o 


OJO, los productos Q:w y V'w no tienen sentido físico y son consecuencia de la eliminación de la 
densidad en la ecuación 


desarrollando la derivada: Qi) w—-Qy¡ w,¡+V- — +w = =0 
. dw 
sustituyendo Qr2w+ (Vo+(Qr —Qo) *t)* q (Q11-Qro) =0 
dw 
reordenando (Vo + (Qí, E Qra) ki t) 7 we (Qu + Qra) -Q/¿w=-w-*Q,1 
w t 
: 1 1 
separando variables ———— dw= dt 
=Q yw Vo+ (Q1, eS Qr>) «t 
wo 0 
- V,+ - .t 
integrando Es, mE |= - «In |— (Qu Q12) 
Qr wo QQ Vo 


Vo 


Qu 
E ar (Q 5 Qr2) .t | | (Q1= Q2) ) 
== Wo | —_—__—>———. 


w(5-min)=0.51 
Vo 


5 
OS 
+ 
A 

| 


despejando 
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Problema 3.13 


Un depósito tiene una capacidad máxima de 400 L y se encuentra inicialmente medio lleno. La 
tubería de descarga que sale de su base es muy larga, de forma que ofrece resistencia significativa 
al flujo de agua a su través. La fuerza que hace que fluya el agua es la debida a la altura de agua en 
el tanque y, de hecho, el caudal de salida es proporcional a dicha altura. Puesto que la altura es 
proporcional al volumen de agua que hay en el tanque, el caudal del agua de salida es: 


Q, 1 = 8 L/min 


Q,, 2 = K"V, donde k = 0.01 min"? 


El caudal de agua que entra en el tanque es constante. A la vista de la información adicional de la 
figura, decir si la cantidad de agua en el tanque aumentará, disminuirá o permanecerá constante. Si 
hay variación, calcular el tiempo necesario para que el tanque se llene o se vacíe por completo. 


El sistema funciona en régimen no estacionario. 


Volumen inicial en el tanque: V,:=200 L 
Volumen máximo del tanque: V mar i=400 E 
Caudal de entrada: Q1,:=8 ma 
main 
Caudal de salida: Qu.:=k-V k:=0.01-min 
SOLUCIÓN 


A la vista de la información proporcionada, el caudal de entrada es superior al caudal de salida, por 
lo que se llenará, así pues se puede calcular el tiempo que tardará en alcanzar el volumen máximo. 


En un sistema de densidad constante, el balance total de materia puede expresarse: 


d 
Qr2- Q11 + V=0 
di 


Teniendo en cuenta la expresión del caudal de salida y operando: 


d 
k-V-Qp¡ +—V=0 
dt 
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== 


me —k.V 
Integrando y despejando el tiempo: t(v)= Al sn e k a) 
Lo **Vvo 


t(Vinaz) =40.55 min 


despejando el volumen: V(t) o (Or, = (Qí1 —k-Vo)* et") 


VA O 
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Problema 3.14 


Un tanque perfectamente agitado contiene inicialmente 100 kg de una disolución al 50 % de Na2SO4 
en agua. En un instante dado se introduce en el tanque una corriente de disolución de Na2SOs al 15 
% con un caudal de 12 kg/min, al mismo tiempo que se extrae del tanque una corriente con caudal 
constante de 10 kg/min. Calcular: 


(a) Una expresión que permita calcular la masa de disolución presente en el tanque en un 
instante determinado. 

(b) Una expresión que permita calcular la concentración de Na2SOs en el tanque en un instante 
determinado. 

(c) La concentración de Na2SOs en el tanque a los 10 minutos de comenzar la operación. 


Datos y notas: 
Todos los porcentajes son en peso. 


m;, = 12 kg/min 
15% Na,SO, 


m, = 10 kg/min 


El sistema funciona en régimen no estacionario. 


Masa y concentración inicial en el tanque: M,:=100 kg w,:=0.50 
7 kg 
Caudal y concentración de entrada: m,:=12 —= w,:=0.15 
min 
Caudal y concentración de salida: m,:=10 9 Wy=Ww 
main 
SOLUCIÓN 


La masa total en el tanque vendrá dada por el balance total de materia: 


d Necesitaremos esta ecuación en 
el balance de componente 


S M t 
—————» fram= [ 1d: 
M_—Ma 0 


Necesitaremos esta ecuación en 


a M=M),+ (m, — mp) -t 

(a) ar 2 el balance de componente 
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El cálculo de la concentracion con el tiempo vendrá dado por el balance de componente: 


E A 
dt 


Sustituyendo wz = w y desarrollando la derivada 


mac aj m)=0 


sustituyendo las dos euaciones anteriores: 


maca (Mot (ra) 


Reordenando y simplificando: 


(os (rm) 0) S 0] mo mamma mn (0) 


dt 
w t 
z 1 1 
Separando variables: du = dt 
my + (w, 1) My + (m, — mp) +t 
wo 0 


Integrando: 


Simplificando y reordenando: 


m + (w, —)) Ma 
—m, 
Mo + (mm) «ty 0 
mi (w, —w) = mi lá (rw, wp) > 
Mo 
—m; 
Mo+ (m, —m,) .t ta 
w,—w= (w, —W) * 
Mo 
—my 
M + (m —m )"t eat 
: 0 1 2 
Despejando: (b) w (t):=w, — (w, —w) * 
Mo 
(c) w(10 min) =0.27 
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BLOQUE || 


Balances de energía total 
en sistemas sin reacción 
quimica 


Problema 4.1 


Se bombea agua desde el fondo de un pozo situado a 45 metros bajo tierra, con un caudal de 
900 litros/hora a un tanque de almacenamiento abierto a la atmósfera, situado a 10 metros sobre 
tierra. Para evitar que el agua se congele durante el invierno, una pequeña resistencia suministra 
9000 W al agua. Durante el transporte desde el pozo al depósito, el sistema pierde calor hacia los 
alrededores a una velocidad de 7000 W. Para elevar el agua se utiliza una bomba con una potencia 
de 1500 W., De la energía consumida por la bomba, el 55% es aportada al agua y el resto se disipa 
en forma de calor hacia los alrededores. Si la temperatura del agua en el pozo es de 2 “C, calcular la 
temperatura a la que entra el agua al tanque de almacenamiento. 


Datos y notas: 


Despreciar la variación en la energía cinética. Calor específico del agua: 4200 J/(kg'K). 


Fluido: 
Calor específico agua: 
Densidad agua: 


Caudal volumétrico: 


Caudal másico: 
Altura inicial: 
Altura final: 

Incremento de altura: 


La temperatura en el pozo es: 
Calor suministrado: 


Calor perdido: 


Agua 
c,:=4200 
kg * 
A ea 
3 
os E = bso E 
hr Ss 


m'"=0,-p=025 E 
Ss 
2, :=-45 m 


z,:=10m 


Az:=z,-2/=55 m 


T,:=2 *C 
qy:=9 kW 
qpi==7 kW 


El calor aportado al sistema es positivo, mientras que el calor perdido por 


el sistema es negativo 


Potencia nominal de la bomba: 


Rendimiento de la bomba: 


Potencia real bomba: 


La variación de energía cinética es despreciable: 
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Nominal := 1500 W 
y :=0.55 
Wieal == Wknominal n= 825 W 


2 
a 
2 10 
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SOLUCIÓN 


El balance de energía total, en estado estacionario, cuando el caudal másico de entrada y salida 
son iguales y no hay cambios de fase es: 


2 
cy AT +g dz Ed == qe + Wes siendo: dep =L Mia dE 
2 0. m m 


Sustituyendo términos, el balance queda: 


q,+ dp E Wieal 
m m 


Cp * (T,—T;) +8 (2,-2,)= 


y despejando T2 
+9,+0, «(2,2 
T,:= 7, 4 20490 + Wes 8922721) 46 ec 
me Cp Cp 
AS T,= T,=2.6 AC 
Se comprueba que la temperatura apenas ha cambiado. 
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Problema 4.2 


Un fluido cuya capacidad calorífica específica es de 2,09 kJ/(kg-9C) entra en un cambiador de calor a 
una temperatura de 5 *C. La temperatura de salida es de 30 *C. El tubo de salida está situado a 7 m 
sobre el tubo de entrada. La velocidad del fluido a la entrada es de 3 m/s y a la salida de 30 m/s. El 
sistema funciona en régimen estacionario. El régimen de circulación a la entrada y a la salida es 
turbulento. Calcular: 


(a) Cuántos kJ serán necesarios para calentar el fluido de 5 *C a 30 *C. 

(b) Cuántos kJ serán necesarios para elevar el material. 

(c) Cuántos kJ serán necesarios para aumentar la velocidad del fluido. 

(d) Qué porcentaje de la variación total de energía que experimenta el fluido corresponde a la 
energía cinética y a la energía potencial. 

(e) Calcular el calor que hay que suministrar al sistema para que éste funcione en las 
condiciones mencionadas. 


DATOS 
Calor específico del fluido: C,*=2090 E 
g . 
Temperatura de entrada T,:=5 *C 
Temperatura de salida T,:=30 *C 
Diferencia de alturas: Az:=7 m 
Régimen turbulento a:=1 
Velocidad de entrada: v¡:=3 cl 
Ss 
Velocidad de salida: v,:=30 pd 
Ss 
Estado estacionario 
SOLUCIÓN 


En estado estacionario, con una única corriente de entrada y otra de salida, el balance de energía 
total es: 


2 


Aiagadak «A St Wasp siendo: 1» wm, =% 
2 a m m 


Se deben calcular cada uno de los términos del balance de energía 


(a) La entalpía específica a una temperatura, si no hay cambios de fase y considerando el calor 
específico constante, es: 


h=h* +0, (1-T*) 
y el incremento de entalpía se puede relacionar con el incremento de temperatura: 


Ah=h,—h,=(n' +c)* (7,7) (10 +c,* (7,7%) 


Ah:=c, + (T,—T¡) =(5.225 + 10*) A 


kg 
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(b) El incremento de energía potencial es: 


Aepo*=8 + 47=68.6 cd 
kg 
(c) El incremento de energía cinética es: 
2) 2 
V —V 
Le. A E 
2-0 kg 
(d) El cambio de energía total es: 
ay y 
Lesoto :=4h+4epo1+ Aci =(5.276 101) Le 


por lo que las razones de cada tipo de energía respecto al cambio de energía total son: 


E )=00% Gn )=0.13% ¡E )0ssx 


Á total Á total Ctotal 


(e) Para determinar el calor necesario, se vuelve a plantear el balance de energía: 


Qesp :=4h Hero FC on == (5.276 ó 101 ) e 
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Problema 4.3 


Un gas perfecto, de peso molecular 20, circula turbulentamente a través de una conducción de 
hierro de 2 pulgadas de diámetro interno. El gas, que circula con un caudal másico de 1360 kg/h, 
entra en la conducción a una presión absoluta de 68 atmósferas y a una temperatura de 68 “C, 
abandonando la misma a la presión atmosférica. 


(a) Considerando circulación isotérmica, calcular las velocidades de entrada y salida del gas, así 
como el calor suministrado al mismo en kJ/h. 
(b) Considerando circulación adiabática, calcular la temperatura de salida del gas. 


Datos y notas: 
El calor específico molar del gas a presión constante, puede considerarse constante e igual a 
29.3 kJ/(kmol"9C). 


DATOS 
Gas perfecto: R:=8.314 db PM:=20 £ 
mol + K mol 
Cal if a nn J Debe convertirse a másico 
alor específico molar: Com *=29.3 Pa para usarse en el BE 
sE Las Com 3 J 
Calor específico másico cy := == (1.5 +10 ) — 
Regimen turbulento: aj:=1 a7):=1 
Diámetro conducción: D:=2 in= (5.1 . 10%) m 
Superficie conducción: iso = (2-10) m? 
Caudal másico: m:=1360 Eg 
hr 
Caudal molar: F:= "189 Po 
PM s 
Presión de entrada: P,:=68 atm=(6.9 + 10%) Pa 
Presión de salida: P,:=1 atm=(1-10*) Pa 
Temperatura de entrada: T,:=68 *“C 
SOLUCIÓN 
En estado estacionario, con una única corriente de entrada y otra de salida, el balance de energía 
total es: 
2 
Ah +8 + Áz + a 4 La = desp + Wosp siendo: desp = EE Wosp = td 
2 a m m 
(a) Circulación isoterma (se anula el término de entalpía): T,:=T, Ah=0 
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Como no hay ningún compresor se anula el término de trabajo: Wo; =0 
En gases, el termino de energía potencial se suele despreciar: g+ Az = pequeño 


Para calcular el calor intercambiado, debemos calcular el término de energía cinética y 
aplicar el balance de energía 


Caudal volumétrico y velocidad a la entrada: 


PA 7.8107 mo iS 9 =38 25 
vl 1 
1 Ss S s 
Caudal volumétrico y velocidad a la salida: 
F+«R-T zi 
Q,):= 2=052 py:= % 22609 ” 
P, s S Ss 


No puede comprobarse el régimen de circulación, calculando el Reynolds, porque no se 
dispone de la viscosidad. 


A partir del balance de energía: 


2 2 
q:= - -m=12.9 kW q=(4.6-107) 
2-0) 2-0; 


(b) Circulación adiabática (no hay intercambio de calor): q=0 


Como no hay ningún compresor se anula el término de trabajo. 
En gases, el término de energía potencial se suele despreciar. 


Para calcular la temperatura de salida, debemos calcular el término de energía cinética y 
aplicar el balance de energía: 


Cp" (TT) +8-2+ «(»2 —v,)=0 


2-0 
La velocidad a la entrada es la misma que en el apartado anterior. 


La velocidad a la salida depende de la temperatura, por lo que no se puede calcular 
directamente. 


Y _0, 0) Mm ” _P,:-PM 

2 Ss v2 Pa 2 R-T, 
É m»«(R=T F-R-T 
sustituyendo en y) : E a A 2 


PE) = cy (LT) +» 
2-0 


T,=320.1 K 


Comprobándose que la expansión del gas produce su enfriamiento. 
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Problema 4.4 


Calcular la variación de entalpía que experimenta el agua entre las secciones (1) y (2) del sistema que 
se esquematiza en la figura. 

Las secciones libres de paso de líquido en los puntos antes indicados son S1=0.1 pie” y S2=0.5 pie” 
respectivamente. Una bomba hidráulica comunica al sistema una potencia de 1 CV y el cambiador de 
calor que se observa en la figura elimina 10,000 Btu/h. El agua circula a 258C con un caudal constante 
de 300 L/min. El sistema opera en régimen estacionario. 


Datos y notas: 
Propiedades medias del agua: p =1 kg/L , u = 0.001 kg/(m-s). 


Cambiador de calor 


z, = 85 pies 


DATOS 
Fluido: Agua 
Densidad agua: p:=1000 Ea 
m? 
Viscosidad agua: 1:=0.001 kg 
mes 
Temperatura entrada: T:=25 *“C 
Altura de salida: z,:=85 ft=25.908 m 
Altura de entrada: z¡:=10 ft=3.048 m 
Sección de salida: S),:=0.5 fé? =0.046 m? 
> , 4.5, 
Diámetro salida: D),:= = 0.243 m 
T 
Sección de entrada: S,:=0.1+f1? =0.009 m? 
4.5 
Diámetro entrada D¡:= L =0.109 m 
T 
ei L m* 
Caudal volumétrico: O,:=300 =0.005 — 
min s 
Caudal másico ic bsp=3 EE 
Ss 
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Potencia bomba: W:=1 hp=7457 W 


Potencia específica: 


BTU _ 


Cambiador de calor: q:=—10000 
a 


Potencia específica: 


2931 W 


E 149 E 

m kg 
O 7 
m kg 


La energía aportada al sistema es positiva, mientras que la energía eliminada del 


sistema es negativa. 


SOLUCIÓN 


Planteando el balance de energía en régimen estacionario: 


2 


Ah +80 40404 == Qesp + Wesp siendo: q. = 
104 
Se calculan las velocidades y el régimen de circulación: 
-v,D 
pyi= % 0.538 2 Re =P" """L=5853+10* Re,> 10% 
S; Ss 12) 
.v..D 
Yyi= 2 0.108 ri Re,>10%* 
S) Ss 12) 
con lo que podemos calcular el cambio de entalpía específica: 
2 2 
v v 
Ah := Qesp FWesp — 8 * (2221) + E A =-661 za 
2.0) 2+0; kg 


y el cambio de entalpía: 


AH:=4h=.m=-3307 W 
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- w 
Wesp => 
m 
QA = 1 
Turbulento 
0) :=1 
Turbulento 
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Problema 4.5 


En una conducción horizontal de 5 cm de diámetro interno, rodeada por una camisa de vapor, se 
caliente COz desde 21 “C a 200 “C, con un caudal de 25 kg/min. La presión a la entrada es de 
4,86.105 N/m? y a la salida de 4,73.105 N/. 


Calcular: 


a) La velocidad del gas en m/s a la entrada y a la salida de la conducción. 
b) Calor transferido del vapor al CO» en kJ/kg. 


Datos y notas: 


Calor específico medio del COz entre 219C y 2000C Cp = 0,921 kJ/kg-K. 
Admítase que el CO» en estas condiciones puede considerarse como un gas perfecto. 
Viscosidad media del CO» entre 21 %C y 200 *C = 0,015 cP. 


DATOS 
Fluido: CO» (Gas perfecto) 
ua J 
Calor específico: cy:=921 —— 
kg»K 
Viscosidad: 1:=0.00015 poise=(1.5+ 1075) HE_ 
m>».s 
Peso molecular CO»: PMco,:=44 año 
mol 
Diámetro conducción D:=5 cm=0.05 m 
Superficie conducción: $ D' =0.002 m? 
Temperatura entrada: T¡:=21 *C 
Temperatura salida: T,:=200 *C 
as . EN 
Presión entrada: P,:=4.85+10” — 
m 
A ed s N 
Presión salida: P),:=4.73 +10" — 
m? 
Caudal másico: 0. :=25 8 0.417 
min s 
SOLUCIÓN 
(a) Al ser un gas, el caudal másico se mantiene constante, pero no el resto de condiciones: 


Densidad entrada: 


Densidad salida: 
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(b) 


Velocidad salida: 


Om ? 


Caudal volumétrico entrada: O, :=*"=0.048 Y 
P1 a 
o Om Ai 
Caudal volumétrico salida: Op :=*2=0.079 == 
P2 8 
Velocidad entrada: v¡i= a =2432 2 
Ss 


»y:= 272 40.112 E 
S 


Ss 


Para calcular el calor transferido debemos plantear el balance de energía: 


Ah+4e, +46; = Dep + Way 


Si recordamos que en gases el cambio de energía potencial es despreciable y que no hay 
ningún compresor: 


Je, =0 Wosp =0 
el balance queda: 
2 2 
va v; 
Cp : (T, T;) e e Le 


2-0) 2-07 


Para calcular el incremento de energía cinética se debe comprobar el régimen de circulación: 


.v,D 
Reynolds entrada: Re; A Re,>10* aj:=1 
17) 
Ñ A) 
Reynolds salida: Re, sra 0410 Re,>10* d=l 
17) 


Sustituyendo todos los valores conocidos en el balance de energía, se calcula directamente el 
calor transferido: 
2 2 
v 
Oipi=o (1 =T) +22 165368 E 
2-4) 20; kg 


q = Qosp * O =68.903 KW 
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Problema 4.6 


En una planta desodorizante de aceite crudo vegetal, el aceite se precalienta en un cambiador de 
calor de flujo continuo en contracorriente que funciona en régimen estacionario. El medio calefactor 
es agua caliente obtenida de vapor de deshecho de la planta. El gasto másico de agua caliente a 
través del cambiador es de 2500 kg/h. El agua entra en el sistema a 100 *C y lo abandona a 43*C. 
El aceite crudo fluye con un caudal de 5000 kg/h. Si el aceite entra en el cambiador a 21 *C ¿a qué 
temperatura lo dejará? ¿Y si el rendimiento en la transferencia de calor fuera del 80 %? 


Datos y notas: 

Los tubos de entrada y salida del cambiador están a la misma altura y son del mismo diámetro. 
El sistema está perfectamente aislado de modo que no recibe ni pierde calor. 

Calor específico del aceite 2.09 kJ/(kg-K). 

Calor específico del agua 4.18 kJ/(kg-K). 


DATOS 
Agua (w) Aceite (a) 
Caudal másico m,,:=2500 de m,:=5000 de 
hr hr 
Temperatura entrada T,,¡:=100 *C T,¡:=21 *C 
Temperatura salida Ty/:=43 *C 
Calor específico Cpw:=4180 0 om Cpa+=2090 E 
kg. K kg. K 
SOLUCIÓN 


Planteando un balance de energía en regimen estacionario a cada circuito del intercambiador: 


El calor cedido por el agua es: y ¿=P * Cy * (Ty 9 —T,, 1) ==1.655 +10? W 
Por lo que el calor ganado por el aceite es: da“= q, = (1.655 . 10*) W 
y la temperatura de salida del aceite: a = Qma* Cpa * (Tas —Ta1) 
a A A 
ALIÓ ll AE 5 78 C 
Mo * Cog 


Si el rendimiento en la transferencia de calor fuera del 80%, significaría que solo una fracción del 
calor cedido por el fluido caliente se emplea en calentar el fluido frío: 


Rendimiento: 7 :=0.8 
Por lo que el calor ganado por el aceite es menor que el anterior: 
q.0:=-9..9n=(1.324- 10%) w 


y la temperatura de salida del aceite también será menor: 


dal o 
4az = ma * Cpa* (T,3— Ta) To E T,¡ + ——=066.6 C 
Ma * Cpa 
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Problema 4.7 


En un intercambiador de cilindros concéntricos se introduce una corriente de agua de 5 L/min a 5 “C 
que se quiere calentar hasta 40 %C. Para ello se utiliza otra corriente, también de agua, con una 
temperatura de 80%C, que se introduce en contracorriente y que sale del sistema a 50 “C. Calcula el 
caudal de agua caliente que se necesita si el rendimiento obtenido en la transferencia de energía 
térmica es del 80 %. 


Datos y notas: 
Calor específico del agua: 4180 J/(kg*K). 


DATOS 
Fluido: Agua 
Densidad agua: p:=1000 Ke 
m? 
4 ; J 
Calor específico agua: C,:=4180 —— 
kg. K 
Pr na. L =5 m? 
Caudal volumétrico de agua fría: Oyi=5:—= (8.333 +10 ) — 
: min s 
ria kg kg 
Caudal másico: mp:=0 yy" p=0.083 = my=5 —= 
: Ss min 
Temperatura entrada agua caliente: T./:=80 *C 
Temperatura salida agua caliente: T.,:=50 *C 
Temperatura entrada agua fría: Tp :=5 *C 
Temperatura salida agua fría: Tp :=40 *C 
Rendimiento: y :=0.80 
SOLUCIÓN 


El calor necesario para calentar el fluido frío es: 


gy:="mp= cy" (To Ta) =(1.219- 10%) 


por lo que teniendo en cuenta el rendimiento, el calor cedido por el fluido caliente es: 


q. —L=-1.524+10* W 


1 


Para llegar a calentar el fluido frío, se necesita que el fluido caliente ceda más calor que el ideal. 
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Planteando el balance de energía al fluido caliente: 


y despejando, se calcula el caudal másico de fluido caliente: 


m.= =012% — m=7290 
C,* (T.— T.1) s min 
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Problema 4.8 


Un depósito cilíndrico situado directamente sobre tierra se llena por completo con una disolución de 
un compuesto de Co” y se tapa inmediatamente. La desintegración del Co*% libera energía 
(de modo uniforme en toda la masa líquida) a razón de 812 Won”, proceso que incrementa la 
temperatura de la disolución progresivamente. Desde el instante en que la temperatura de 
la disolución supera a la de los alrededores, tiene lugar una pérdida de calor desde el depósito hacia 
los alrededores (suponer que sólo se pierde por las paredes y la tapa). La disolución está 
perfectamente agitada, de modo que puede admitirse que en cualquier instante posee una 
temperatura uniforme. Si la cantidad de calor por unidad de tiempo que pierde la disolución hacia 
los alrededores puede determinarse por la ecuación: 


q=U+A+ (T- Tam) 


donde: 
q es la cantidad de calor transmitida por unidad de tiempo (vatios). 
U es el coeficiente global de transmisión de calor (w. m0). 
A es el área de intercambio. 
T es la temperatura de la disolución. 
Tamb es la temperatura de los alrededores. 


(a) ¿Cuál será la temperatura de la disolución, a las 24 horas de haberlo llenado, si la 
temperatura en dicho instante era 20C? 
(b) ¿Cuál será la máxima temperatura que puede alcanzar la disolución? 


Datos y notas: 

Admitir que el valor del coeficiente global de transmisión de calor desde la disolución a los 
alrededores, por las paredes y la tapa del depósito, es U = 5.8 W- m?.K*. Dimensiones del 
depósito: 3 m de altura y 2 de diámetro (suponer despreciable el grosor de las paredes). La 
temperatura de los alrededores es constante e igual a 20%C. Calor específico de la disolución: 
5.6 kJ. kg *. Densidad de la disolución: 1040 kg- m?. Considerar despreciable la energía aportada 
por el agitador. 


DATOS 
Generación: Gf¿:=812 Bd 
mi 
Coeficiente global de transmisión de calor: U:=5.8 pe 
m*-K 
Calor específico: C,*=5600 e 
kg. K 
Densidad: p:=1040 kg 
m 
Altura y diámetro del depósito: l:=3 m D:=2m 
Li + D? 3 
Volumen depósito: Vi=1+ = (9.425 «10 ) L 
Masa en el depósito: M:=V¿+p=(9.802.10*) ke 
y . ES + D? 2 
Area de intercambio: Ay=[(7+D+1)+ =22.0 m 
(lateral + tapa) 
La temperatura incial es: T,:=20 *C 
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y la temperatura ambiente: T 


Ecuación de intercambio de 
calor con el ambiente: qe 


SOLUCIÓN 


Estableciendo el balance de energía total en estado no estacionario: 


A(m+h)+ aun) =q+ w+ | 6, av 
E E a dM 
teniendo en cuenta que no hay entradas ni salidas: A(m-h)=0 y ÓN 0 
t 
y que no se aporta energía como trabajo W=0 
LU q [gay 
dt 
Mm q | Gp dY 
dt dt 


Mec E —UA q (T— Tamp) + Gg+ Y 


dT _ dt 
(Gz+V+U-A¿* Tap) = UA¿*T Mc, 


a 


e integrando: - : 
dramas] ! di 
T, 0 


separando variables: 


(Gz+V+U-Ag* Tap) — UA ¿+ T 


a 


1 
—UA4 


an (Gp V+ U + Ag" Tam) — UAg" T = 1 1 
(Gz+ V+U-+A¿* Tan) — UA q * To 


y despejando la temperatura: 


—U- Ay 
Ae 
M-c 


(Gj+ V¿+U- Ag" Tam) — ((Gg V ¿+ U+ Ag" Temp) — U+ Ag" To) = e 
U-Ay 


T(s) := 


T(24. hr) =30.9 *C 


La temperatura máxima se alcanzará cuando se llegue al estado estacionario (ya que no existen 
restricciones). A partir de la ecuación diferencial: 


dT 


==0 0=—UAg* (Tos —Tamp) + Gg Y 
dt 
To i= Ss Ya 800 ec 
est*— Lamb === E 
U de Á d 
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Problema 4.9 


Un depósito está inicialmente lleno con 2300 kg de aceite a 15 “C. En un momento dado comienza a 
introducirse aceite con un caudal de 450 kg/h, también a 15 “%C y empieza a descargarse con el 
mismo caudal. Al mismo tiempo el tanque se calienta mediante vapor de agua a 130 “C que 
condensa por el interior del serpentín, estando el contenido del tanque perfectamente mezclado por 
medio de un agitador. La velocidad con que se transmite el calor desde el serpentín al tanque viene 
dada por la expresión: 


q=U:A- EA E pd 
Calcular: 


(a) El tiempo que deberá transcurrir hasta que la temperatura del aceite que descarga el 
tanque sea de 32 “C, 
(b) La máxima temperatura que puede alcanzar el tanque. 


Datos y notas: 

Energía aportada por el agitador: 1kW. 

Calor específico medio del aceite: 2 kJ/(kg K). 
Valor del producto (U-A): 158 W/K. 


DATOS 
Calor específico: Cy =2000 E 
kg» K 
Masa inicial: M),:=2300 kg Temperatura inicial: Tp:=15 *C 
Caudal entrada: m,:=450 E Temperatura entrada: 7T,:=15 *C 
r 
Caudal salida: m>,:=mM, 
Serpentin: = UA» (Trapo — 1) UA:=158 Y Trapoyi=130 9C 
p s q= vapor ES K vapor *— 
Agitador: W:=1000 W 
Temperatura final: Tina *=32 *C 
SOLUCIÓN 
Dado que los caudales de entrada y salida son iguales: 
Balance de materia: —- M:=M,=2300 kg 
t 
Balance de energía: m,*C,* Tm,» cy" T) Mec E T a+ W 
£ £ 
m=m M=O A a E 
dt dt 
M HL U-A-T U-A+T y T o TI+4W 
*Cp* q” «A. vapor — «A. —mM*C,* +me*C),* TE 
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A (W+U-A-T, +m+c,*T;) -(U-4A+m+c,)+T 


vapor 


M-+c,+* 


dT _ di 
(W+U-+A-+ Trapo +M*Cp*T¡)-(U-A+m-c,)*T Mc, 


T 
1 1 
dT= dí 
| (W+U-A- Trapo + mM+Cy*T¡) —(U-4+m-+c,) + T) | (M-c,) 
To 0 
l In ((7+U-A + Tuapor Em Cp" 11) —(U-A+m-cp)-1) | _ 1 .t 
—(U-A+m-c,) (W+U-A- Trapo + m+*C)*T¡)-(U-A+m+c),)*Tp)) M-c, 


—(UA +m+ Cp) we 
M-c, 


W + UA» Togpoy + m+Cp* T;) —((W+UA + Tyapoy + M+Cp* 1) — (UA +m-+cp) * To) e 
UA +m+c, 


T(t) = ( 


(7) := 


Mcg y [ XOVH VA Tiapor + m0) 71) — (UA + m- cp) +7) 
—(UA +m-+c,) 


((W+UA + Tapoy + m-+C)*T¡) — (UA +m-+c,) + Ty) 


¿(32 2C)=(5.064-10) s  1(32 %C)=1.407 hr 


e dT : 
Cuando se alcance la temperatura máxima, 2 =0, por lo que usando el balance de energía en su 
dt 
forma diferencial: 


m+*Cp* (max — Ty) =UA+ Mi Tia +W 


mc Eg UAT pig VAT aa tr CT + WE 


max 
(m + c, + UA) + Tia: = UA * Tyapoy + mM * Cp * Ty + W 


A 
NO A IS E A E 


de m+cC, + UA 
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Problema 4.10 


Una corriente de 8000 kg/h de una disolución acuosa de NaCl al 7 % en peso, se concentra hasta el 
40 % en un evaporador que funciona en régimen estacionario. La corriente introducida se calienta 
en el evaporador hasta 82 “C. El vapor de agua eliminado de la corriente y la solución concentrada 
salen a 82*C. Como elemento calefactor se utiliza vapor de agua sobrecalentado que condensa, 


suministrando a la disolución 1.71-10” kJ/h. 
Calcular: 


(a) A qué temperatura entra la disolución en el evaporador. 
(b) El caudal de disolución concentrada que se obtiene continuamente. 


Datos y notas: 

Calor de vaporización del agua a 82 “C 2300 kJ/kg. 

Admitir que el calor específico de las disoluciones de NaCl, de cualquier concentración, es igual al calor 
específico del agua (4.18 kJ/(kg-K)) y que la entalpía de las disoluciones de NaCl de cualquier 
concentración coincide con la del agua a 25 *C. 


DATOS 
Corriente entrada: m,:=8000 ze w,:=0.07 T, 
Y 
Corriente salida (vapor): m> w,:=0 T,:=82 *C 
Corriente salida (concentrado): M3 w3:=0.4 T¿:=T, 
Calor aportado por el calefactor: q:=1.71-107 el 
» 
ee ; z ! kJ 
Calor específico disoluciones: C,:=4.18 
kg. K 
OS , kJ 
Calor vaporización agua: L:=2300 Fa 
E 
SOLUCIÓN 


Es un sistema con una entrada y dos salidas. 


Aplicando los balances de materia, pueden calcularse los caudales de las corrientes de salida: 


BM total: M¡=M>,+M; 
BM NaCl: M¡*W¡ =M3*W)+M3*Ww; 
Mw 
EA 
wW; hr 


m,:=m,—m3=6600 l3 
hr 


Vicente Sanz Solana, Nuria Vicente Agut 


ISBN: 978-84-19647-45-0 DOI: http://dx.doi.org/10.6035/Sapiential90 


Índice 


Aplicando el balance de energía en régimen estacionario: A(m -h) =q 
Teniendo en cuenta que hay un cambio de fase en la corriente de vapor: 


m3+Cp* T¿+m)+ (c,* T,+L) =m,+c)*T,=q 


la temperatura de entrada de disolución es: 


m3*C,* T¿4+m)+ lc . T,+L) =Q 
NAT e AA o O 
Mm; Sl Cp 
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Problema 4.11 


Se dispone de un tanque aislado y perfectamente agitado que contiene un serpentín calefactor. Se 
introduce agua en el tanque con un caudal de 0.45 L/s y se extrae de él un caudal de 0.35 L/s. El 
volumen inicial de líquido en el tanque es de 2.8 nm y su temperatura inicial de 21 *C. 

Si el agua entra en el tanque a 65 *C, el elemento calefactor proporciona 50 kW y la potencia 
suministrada realmente por el agitador al líquido es de 4 kW, calcular la temperatura en el tanque 
en función del tiempo y la temperatura al alcanzarse el régimen estacionario. ¿Cuál sería la 
temperatura máxima si el tanque tuviera una capacidad máxima de 5000 L? 


Datos y notas: 
Calor específico del agua 4.18 kJ/(kg-K). 


DATOS 
Densidad agua: p:=1000 Eo 
m? 
Calor específico agua: Ccy:=4.18 + 10? E 
kg» K 
Volumen y masa inciales: V,:=2.8 m* Mo :=V,-p=2800 kg 
Volumen y masa máximos: — V,,,.:=5 m' Mamas :=V max + p =5000 kg 
Temperartura inicial: Ty:=21 *C 
Ea dea ] L kg 
Caudal volumétrico y másico de entrada: Oy,:=0.45 2 m¡:=0Qy,"p=0.45 
Ss Ss 
Temperatura entrada: T,:=65 *C 
os ES ad L kg 
Caudal volumétrico y másico de salida: Oy,:=0.35 2 m>:=0Q0yp=0.35 
Ss Ss 
Calefactor: q:=50000 W 
Agitador: w:=4000 W 
SOLUCIÓN 
Balance total de materia en régimen no estacionario: 
dM dM 
m,-m¡+—=0 UTE 
t 


separando variables e integrando entre las condiciones inicial y genérica: 


[rien M(1):=M,+ (m,—m)) «1 


M, M, 
max 0 
bag iS =22000 s 
MM) 
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Balance de energía en régimen no estacionario: 


d(M-c,+T) 
M7 + Cy * T=M* CyT += q+w 
dt 
derivando: m,+cy*T=m,=c)*T,+M-+c)" O Mdw 
dt dt 
dT 


sustituyendo 


m>+C)* T—m,+C)+T,+(M,+ (m,—m)) +1) =c . 


reordenando y dT_(q+w 
simplificando (M,+ (m, —m>) +1) + di al A +mT,)=m, T 
Pp 
dE t 
1 1 

separando A A A 
variables e A A, My + (m,—m)) «1 
integrando Es 0 


—m, 
mM, —M) 


r()=_. [Emo] (Lim 1] 1) (Etna 


m;¡ Cp Cp Mo 


Para un tiempo infinito, la función anterior daría el siguiente resultado: 


Dn q+w 


est* 


=9) “E 


m;¡*C, 


Sin embargo, el tanque tiene un volumen finito y se llenará antes de llegar al estado estacionario. 
Como consecuencia, la temperatura máxima alcanzable será la correspondiente al llenado del 


depósito. 


Tltrax) =88.4 9C 
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Problema 4.12 


Un tanque perfectamente agitado contiene 500 litros de una solución acuosa de NaCl al 10 % en 
peso a 20 *C. En un momento dado comienza a entrar una corriente de agua pura a 35 *C con un 
caudal de 10 L/min, al mismo tiempo que se extrae una corriente de disolución con un caudal de 
8 L/min. En este instante comienza a calentarse el líquido mediante una resistencia eléctrica que le 
suministra una potencia de 5 kW. Calcular la temperatura, composición y volumen del líquido en el 


depósito al cabo de una hora. 


Datos y notas: 


Considerar que a 25 “C la entalpía de cualquier disolución es igual a la del agua pura, y que la 
densidad y el calor específico de la disolución de NaCl a cualquier temperatura y composición son 
iguales a las del agua pura. Densidad = 1 kg/L, Cp = 4180 J/(kg-K). 


DATOS 
Densidad agua: p:=1000 E 
Calor específico agua: Cy:=4180 E 
kg * 
Volumen inicial: V,:=500 £ 
Masa inicial: M),:= V,-p=500 kg 
Concentración inicial: wp,:=0.1 
Temperatura inicial: Tp:=20 *C 
Caudal volumétrico entrada: O¡,:=10 En 
Caudal másico entrada: a O,,*p=0.167 kg 
Concentración entrada: w,:=0 " 
Temperatura entrada: T,:=35 *C 
Caudal volumétrico entrada: O;,:=8 — 
Caudal másico entrada: pe += 01, *p=0.133 kg 
Concentración entrada: W),=w " 
Temperatura entrada: T,=T 
Calefactor: q:=5 kWw=(3-10*) as 
min 
SOLUCIÓN 
Balance de materia: 
dM dM 
mM) A UA 
M(1):=M,+ (m,—m)) +1 M(1 hr)=620 kg 
y el volumen en una hora será: y: MU 11) _ 620 JL 
Pp 
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Balance de componente: 


ma 2 
dt dt 
m3 04 (Mp + (mm) +0) 40 (m,—m,) =0 siendo: W) =W w=0 
É 
separando variables e integrando: 
w E 
1 1 
dw= di 
(—m, -w) (M,+ (mm, m)) -1) 
wo 0 
A E A A a 
m; wo m¡—M) Ms 


—m, 


a 


w(1-hr)=0.034 


Balance de energía: 


dT dM 
m),»Cc,*T,=m,*C,*T,+M*cC,+ +T*c,+ =q 
2 "€p 2 1 "Ep 1 p di Pp di 
dT q , Ñ 
m>-T—=m,-T,+(M,+ (m,—m)) -1) —+T+ (m,-m)) == siendo: T,=T 
dt Cp 
separando variables e integrando: 
dls t 
z dT= z di 
[Ls moT,= m7) (M,+ (m, —m>) +1) 
C 0 
P 
Ta 
Em T ¿2 m,+T 
ll Cp s > a) 
eN E A My 
Cp 


1 
T(1) := . Es emyoTy [Lam T 47) 
mi; Cp Cp 


A E 
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Problema 4.13 


Un depósito, cuya capacidad máxima es de 800 litros, se encuentra inicialmente lleno con 200 kg de 
agua a 25 “C. En un instante dado se introducen en el tanque dos corrientes; una de disolución 
acuosa de A2Cr70, al 70 % con caudal de 65 kg/min a 47 “C y otra que lleva una disolución 
de KCrx0, al 45 % a 56 %C y con un caudal de 80 kg/min. Simultáneamente se extrae una 
corriente de mezcla a razón de 100 kg/min y se calienta el depósito mediante una resistencia 
eléctrica a razón de 3 kW. Calcular: 
(a) Una expresión que nos dé la variación de la masa total del depósito en función del tiempo. 
(b) La expresión que indica la variación de la fracción másica X¿Cr,0, en el interior del 
depósito con el tiempo. 
(c) Una expresión que indique la variación de la temperatura en el interior del depósito con el 
tiempo. 
(d) El tiempo de llenado del depósito. 
(e) La composición y la temperatura de la mezcla en el instante en que se llena el depósito por 
completo. 


Datos y notas: 
El depósito está perfectamente agitado y aislado. Las propiedades de las disoluciones de XCr,0, 


pueden considerarse iguales a las del agua pura: (p=1 kg/L; Cp=4200 J/(kg *C)). 


DATOS 
Densidad y calor específico de las disoluciones: p:=1000 ¿HE c,:=4200 + el 
m? kg» K 
Condiciones iniciales: — M,:=200-+ kg w,:=0 Ty+=25 *C 
Corriente entrada 1: m,:=65+ ME w,:=0.7 T,:=47 *C 
min 
Corriente entrada 2: m,:=80+ EE w»,:=0,45 T,:=56 *“C 
min 
Corriente salida: m3:=100 E W3 =w T¡=T 
min 
Límites del sistema: Vera: += 800 + E Mar := V max + p=800 kg 
Calor aportado al sistema: q:=3 kW= (1.8 . 10%) e 
min 
SOLUCIÓN 


Planteando el balance total de materia en régimen no estacionario: 


daM dM 
m3 (m,+m,) +=0 ——=M,+M),—Mj 
dt dt 
M t 
dM= (m,+m,—m3) + dt f1am= (m,+m,m) + | 1as 
Ma 0 
M—M),= (m, +m,—m3) + M(1):=M,+ (m,+m,-m3) -1 
Mina: Mi 
a OS lina =13.333 min M (tinas) = 800 kg 
M¡+M)—M; 
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Planteando el balance de componente en régimen no estacionario: 


ia 
dt 
dw dM 
M¿"W—(M¡-W¡+M,*W>) +FM*-——4+we—= 
3 ( JW; 2 2) di dí 


myo0— (m0, m0) + (Mp + (1 +m,=m3) +0) 40 (m, +m—m,) =0 


(1 001+m39 0) + (Mp+ (1, +m,=0m) +1) 400 (m,+m,) =0 


(My + (mmm) +0) (m,-w,+m3-w>)) — (Mm, +m)) -w 
dw dt 


(m, -w, +m,+w>) — (1, + m3) -w _ Mo+ (m,+m,—m) ++ 


w t 
1 1 
dw= di 
jara il saaacay 
0 0 


AO 


My+ (m,+m,—m3) +t 
(m,+m, (m, -w, +m3+w>) 


al 
m, + m,—m3) M, 
=1 


1 
md Near) E (enema) 


(m, -w,+m)+w>) Mo 
M, +(m,+m>,—m>) +t (1, +m>— m3) 
mM], +m)+W— (MM, +W, +m>-w)) + 0 ( 1 2 3) 
M, 
w (e) == 
mM,+M), 
w (tana) =0.556 
Planteando el balance de energía en régimen no estacionario: 
CE 
di 
m3*C,* T— (m)+cyT, + m0 cy T) +Mocy 4 T-c,> 3 
t t 


m3z+T— (1 7,+m3-7,) + (Mo + (m,+m,m,) 0) 4 T- (m,+m,—m3) =4 
Cp 


(m7, +m3+7,) + (Mo + (1; +m,m3) +1) E T-(m,+m,) =-L 
C 


o] 
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a — (m, + m)) + T 
dt E 
dw - dt 
(mii em>r3+2 — (m,+m)) + T My+ (m,+m,=m3) «1 
Cp 
w 1 
a dT= - di 
(m7, +m=73+2 — (m,+m)) + T ? My + (1, +m)=m3) ++ 
Cc 
J Pp 
-T PE, E de PR «T 
y , pa ¡EM da+ a (m,+m)) _ 1 e ] 
(m,+m) (m7, +m=73+2 — (m,+m)) -T, (m,+m,=m;) MG 
En 
-1 
(1, + m>) 
A O ap a 


A 
M, 
(m7, +m>7,+2) =(m,+m)=1 a 
c 
Pp 


M,+ (m,+m,=m3) «t 


—(m, + m>) 


momen) 


mpTimpTo Le (mis Temor 
Cp Cp 


M) 


] (m, +m>,- m3) 


TOS 
M¡+M) 


M1 RO 
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Problema 4.14 


Se dispone de una disolución concentrada de un polímero al 50 % en peso. Para poder proceder a 
su venta es necesario diluirlo al 40 % en peso, lo que se hace mezclándolo con agua en un tanque 
perfectamente agitado que tiene una capacidad máxima de 1000 litros. 

(a) Si el tanque contiene inicialmente 100 litros de agua pura y se alimenta la disolución 
concentrada con un caudal de 10 L/min, calcula el tiempo necesario para alcanzar una 
concentración del 40 % en el tanque y el volumen que contendrá en dicho instante. 

(b) Una vez alcanzada la concentración anterior (40 %), se introduce, además, una corriente 
adicional de agua pura en el tanque que debe mantener la concentración invariable durante 
el resto del llenado. Calcula el caudal de agua que hay que introducir a partir de dicho 
instante, así como el tiempo total de llenado. 

(c) Por otra parte, para facilitar el envasado de la disolución resultante, se ha instalado en el 
tanque un elemento calefactor de 30 kW que se conecta cuando empieza a introducirse la 
segunda corriente. Calcula la temperatura final alcanzada si la temperatura inicial de todas 
las corrientes y del agua contenida en el tanque es de 30 %C. 


Datos y notas: 


Considerar que todas las disoluciones tienen una densidad de 1 kg/L. 
Calor específico de las disoluciones: 4.18 kJ/(kg K). 


DATOS 
Mama := 1000 kg  Mp:=100 kg wp:=0 wy:=0.4 m,:=10 He 
min 
e,:=4180 E w,:=0.5 
kg»-K 
SOLUCIÓN 


Etapa 1. El sistema tiene una corriente de entrada (disolución concentrada) y no tiene ninguna 
corriente de salida. En el tanque hay una cantidad inicial de agua pura. 


Planteando el balance total de materia en régimen no estacionario: 


d d 
—m¡+—M=0 —M=m, 

di di 
M t 
fram=m,- f 11 M-M,=m;+t 
Ma 0 


Planteando el balance de componente en régimen no estacionario: 


mit (MA =b —m,w¡+M-+ q w|+w+ E =0 
di dí di 


d 
or (os) 
Í 


1 1 
dw= | -—————di£ 
| (m, -w,—m,-w) | (M,+m,-1) 
0 


A 


lm..waA 
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Mo Mo w; 
t):=w,-— «(w— t(w):= —M, 
1) ds M+m;,*t ve 10) 0») e ) m; 
1(0.4) =40 min M(40 min)=500 kg 


Etapa 2. En esta etapa se añade una nueva corriente de entrada (agua pura). 
El caudal de agua que debe añadirse para mantener constante la composición será: 


M¡*W; 


a E 


m,+=w,+m,-0= (m,+m>) 0.4 m>,:= Es 


Las condiciones iniciales de esta etapa son las condiciones finales de la etapa 1: 
Mp, :=M(40 + min) =500 kg 
El tiempo de llenado restante es, aplicando el balance de materia: 
—M,—M)3+ E 0 
dt 
M)(t):=My,+ (m, + m>) *f 


Mmnax 7 Mg 


0 min 
m¡+m) 


Planteando el balance de energía en régimen no estacionario a la etapa 2: 


T,:=30 *C Tp :=30 *C O :=30000 -W 
dT dM 
—m,*C,*T,—-m>,*C,*T,+M+c,*— +c,+T- =0 
di ds Pigp E dt 
dT 
—(m, +m)) Cp Ti+ (M,>+ (m,+m,) -1) O (m,+m)) =0 
dT 
(M,,+ (m, +m)) +1) - =2 + (m,+m3) -T,— (mm) +T 
dt Cc, 
t 
: dT= | E di 
[+ (m,+m):7,=(m,+m)-7] (Moz+ (m;+m7) ) 
Cp 0 
To 
1 Mo) 
T(0):= L + (m,++m>) -T,- (24 (m,+m)) =T,—(m,+m>) +To 
m¡+m), Cp Mo + (m,+m)) +1 Cp 
T(40 min) =47.2 *C 
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Problema 4.15 


Desde un lugar elevado circula agua a través de una conducción hasta una turbina, y a continuación 
fluye por otra conducción de igual diámetro que la anterior. En un punto situado a 92 m por encima 
de la turbina, la presión es de 210 XNW- m?, mientras que en un punto situado 3 m por debajo de 
ésta, la presión es de 120 XAN- m?. El agua circula con un caudal de 900 kg/s y la potencia que se 
obtiene de la turbina es de 740 kW. Sabiendo que la turbina suministra el 90 % de la energía cedida 
por el fluido calcular: 


(a) La pérdida de energía mecánica en el circuito. 
(b) El aumento de temperatura que experimenta el agua si no hay intercambio de calor con los 
alrededores. 


DATOS 
Densidad agua: p:=1000 EA 
m? 
Calor específico agua: — c,:=4180 e 
kg-K 
Altura inicial: z/:=92 m 
Altura final: z3=3 m 
Diferencia alturas: Az:=23-Z/¡=-—95 m 
Incremento energía potencial específica:  g+-4z=-—931.632 - 
8 
Presión inicial: P,:=210 Ea (2.1-10%) Pa 
m? 
Presión final: Psw= 120 SA (1.2.10%) Pa 
m 
P,-—P 
Incremento energía de flujo específica: si). E 1 
p kg 
Caudal: m:=900 He 
S 
Turbina: Wiom==740 kW 1 :=0.90 
Wnom 5 Wieal J 
Wal := =-8.222.10” W Wesp:= =-913.58 L 
1 m kg 
SOLUCIÓN 
: Does 1 Y AP 
Aplicando en balance de energía mecánica: pnl sa . +—= Was =y 
a p 
P,-—P 
despejando: ey=W esp 8 * (2321) — 271 108.052 L 
p kg 
Las pérdidas de energía mecánica son: Ey:=ey-m= (9.725 . 10*) W 
y el calentamiento sería: ey=C,*4T AT:= Y =0.026 49C 
Cp 
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Problema 4.16 


Mediante una bomba se hace circular de modo continuo agua a 15 “C a través de una conducción 
de 5 cm de diámetro interno. Un dispositivo registrador indica un caudal de circulación de 10800 L/h 
con una carga de 55.2 cm de agua a través de un diafragma medidor. 


(a) ¿Cuál será el diámetro del diafragma? 

(b) ¿Cuál será la pérdida neta de carga debida a la presencia del diafragma cuando el caudal de 
circulación es de 7200 L/h? 

(c) La pérdida neta de carga en un venturímetro es 1/8 de la del diafragma, pero costaría 120 
euros más. Si esta diferencia se tuviese que pagar en 5 años, el rendimiento total de la 
bomba es del 60% y el coste de la energía es 0.06 euros/kWh, ¿que resultaría más 
económico, instalar un venturímetro o un diafragma? 


Datos y notas: 
Viscosidad del agua a 15 %C: 1 cP. 


DATOS 

Densidad agua: p:=1000 Ea 

m? 
; ; . A kg 
Viscosidad agua: 1:=0.01 poise=0.001 
mes 

Temperatura : T:=15 *C 

Diametro conducción: —D,:=5 cm 

Rendimiento bomba: 9 :=0.6 

Caudal volumétrico: Oy:= 10800 e =0.003 Y 

hr Ss 
E kg 
Caudal másico: m:«=0y.p=3 2 
Ss 
Lectura manométrica: — 4h,,:=-—55.2 cm 
Diferencia de presión: 4P:=p+g+4h,,=-5.413+10* Pa 

Unidad monetaria: EUR:=1 
SOLUCIÓN 
(a) Por ser un líquido: Y = 1. Se supone C = 0.6 y luego se confirmará si Rez > 10000 

-2+p- AP 
C:=06  Y:=1 DN A 
(1-f*) 
Resolviendo la ecuación: D,:=0.039 m 
mn Bs D, 
y la relación de diámetros: a, de 
1 
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Comprobando si la suposición de que el régimen de circulación en el estrechamiento es 


turbulento 

S,=%.D? =0.001 m* 2 2,511 A 

4 S> s 

p*v2*D, 4 mia 
Re, := “2 ===9.794-+10 Suposición de C=0.6 correcta 
17 
3 
(b) Nuevo caudal: O), 7200 =0.002 2 mit E 
Y Ss Ss 


Despejando: — 4P,:= =-2.452+10* Pa 


m> a 
75) 2) 


Se comprueba nuevamente el valor de C: 


veD 
y, 22 1.674 Edd E 
S> s 17 
Suposición de C = 0.6 correcta. 
(c) 
Pérdida neta de carga en el diafragma: AP verm.dia *= (1 8) *4P,=-—960.216 Pa 
ed , a e Aita _ 
Pérdida neta de carga en el venturímetro: — 4APeym vent: === 120.027 Pa 
8 
Diferencia de pérdidas: AP erm.dia — APperm.vent=—840.189 Pa 
Á ig = AP, , 
Consumo_extra_energia :=—P 27240 Pm m,=-1.68 W 
1) 
Coste_extra_periodo := Consumo_extra_energia - 3600 A HE 365 20, S yr+ 2 0.06 ER 
hr day yr 7 KW-+ hr 
Coste_extra_periodo =—7.36 EUR 
El consumo extra en el periodo es menor que 120, por lo que no compensa el cambio. 
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Problema 4.17 


Por una conducción de 15 cm de diámetro interno circula aire a 40 %C con una velocidad media de 
1 m/s. Se ha instalado un manómetro de Hg en U de rama abierta que señala una diferencia de 
alturas de 15.1 cm. El aire pasa a través de un orificio medidor y se mide la diferencia de presión 
mediante un manómetro diferencial de dos líquidos ( CC/, y agua). 


(a) ¿Cuál ha de ser el diámetro del orificio para que la lectura del manómetro sea de 1.7 cm 

(b) Si circula aire a 20 “C y el manómetro de mercurio indica una diferencia de alturas de 
76 cm, ¿a qué velocidad ha de circular el aire para que el manómetro de dos líquidos 
siga señalando 1.7 cm? 


Datos y notas: 
Presión atmosférica: 760 mmHg. Densidad del mercurio: 13600 kg- m7? Densidad del CCl,: 1590 


Kkg- m7. Relación Cp/Cv para el aire: 1.4. Viscosidad del aire: 1.87 - 10? kg/(m's). 


DATOS 
Fluido: Aire 
Temperatura: T,:=40 *C 
Velocidad: v¡:=1 EE 
Ss 
Diámetro conducción: D,:=15 cm=0.15 m 
3 
Caudal volumétrico: O, :=v; gs D? =0.018 tn 
Ss 
, Es kg 
Densidad mercurio: Pug"= 13600 a 
m 
; . kg 
Densidad CCI4 : Pccig:= 1590 = 
m? 
z . kg 
Densidad agua: Pagua*= 1000 —= 
m? 
Peso molecular aire: Maire:=28.9 db 
mol 
Viscosidad aire: pi=1.87 107 F8_ 
m»s 
Razón calores específicos: y:=1.4 
Presión atmosférica: Pam =760 mm Hg 
Pam =P ein *PHa"E= 1.014-+ 10? Pa 
Presión estática conducción: Ah,,:=15.1 cm 
(despreciando la columna de Py3=P agp + Ap 98 *Pug= (1.215-105) Pa 
aire) 
Diferencia de presiones enel  4h,,,:=1.7 cm 
estrechamiento: AP:= —(pecrs —Pogua) -2+Ah,,,=-—98.361 Pa 
SOLUCIÓN JE 
(a) Para fluidos compresibles, como el aire, el caudal viene dado por: Q,=C-Y-4),- an a 
P]* (1 =P ) 
Densidad del aire en la conducción: pyi= PrMaire _ 1.349 kg 
R ld T; m? 
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Suponiendo que el régimen turbulento, el coeficiente de expansión es: C:=0.6 


Para orificios, el coeficiente de expansión es:  Y=1- - A -(0,414+0.35+p*) 0.2<8<0.8 
ET 
Sustituyendo: Q,=C:+ ñ IE Codi +0.35+ e] 1) PE 
Era E 
Se tantea numéramente o gráficamente: D),:=0.055 m 
Se comprueba la condición para el cálculo de Y: B=2=0367 (0.2<f$<0.8) 


1 


Y se comprueba la condición para el coeficiente de expansión supuesto: 


Presión en orificio: P,:=P¡+4P= (1.214. 10*) Pa 
P>-Ma; y 
Densidad aire en orificio:  p, ET ÓN 1.348 e (supuesto isotermo) 
*L7 
Caudstinásico: o 0.004 kg en gases se mantiene constante, 
Ape nmanlís: m:=0,+p,=0. A no así el caudal volumétrico 
3 

Caudal volumétrico y velocidad is 2 =0018 v):i= 4 =7.444 Y 
en orificio: P2 sd (5-0 3 d 

2 

4 
.vy.D 
Reynolds en orificio: Re), ¡= 227 12"22 =9.951-10* Re,> 10* 
17 


por lo que pueden aceptarse todos los cálculos realizados. 


(b) Las nuevas condiciones de circulación son: 
Temperatura: T¡,:=20 *C 
Lectura manométrica: — 4h,,¡,:=76+cm 
Presión absoluta en conducción: —P;,:=Pg + Alp" 2 *Pug= (2.028 . 10%) Pa 


Volviendo a aplicar el balance de energía mecánica: 


. Pip Main k 
Nueva densidad: dis E A de 
R dl Ti, m? 
4 
mm AP D, 
Factor expansión: Y:=1-— +[0.41+0.35+ =1 
py D; 
e ) —-2+ AP mi 
el caudal volumétrico será: O,yp:=C+Y+ (5-2) a =0.013 
Pr (1=f$*) a 
la velocidad: ás m 
y yy 21% 0737 
To. D,) Z Ss 
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Problema 4.18 


Por una conducción de 300 mm de diámetro circula agua con caudales no superiores a 0.1 m/s. El 
caudal ha de medirse con un venturímetro y la diferencia de presión entre la entrada del 
venturímetro y su cuello del mismo se medirá con un manómetro en U que contiene Hg como fluido 
manométrico. Si la lectura máxima del manómetro no ha de sobrepasar los 600 mmHg, calcular el 
diámetro mínimo que ha de tener el cuello del venturímetro, suponiendo que el coeficiente de 
descarga es constante e igual a 0.97. 


Datos y notas: 
Densidad del mercurio: 13600 kg/ nm. 


DATOS 
Fluido: Agua (líquido) 
Caudalímetro: Venturímetro 
Densidad agua: Pagua *= 1000 E 
m? 
; dad kg 
Densidad mercurio: PHg'= 13600 a 
m 
Diámetro conducción: —D,:=300 mm=0.3 m 
m? 
Caudal volumétrico: O,:=0.1 —= 
Ss 
Manómetro: en U con mercurio 
Lectura manométrica: — 4h,,:=600 mm 
Diferencia de presión: — AP:=—(Pxg —Pagua) "2 * AM ==—7.414+ 10% Pa 
Coeficiente de descarga: —C:=0.97 
Factor expansión: Y:=1 Fluido incompresible 
SOLUCIÓN 


Aplicando el balance de energía mecánica al venturímetro: 


2» AP 


D 4 
¡A Ma 
D) 


Pagua úl 


Se calcula el diámetro del venturímetro por tanteo: 


D,=0.103 m 
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Problema 4.19 


Por un gaseoducto circula metano con un caudal másico de 7.48 kg/s a 15.5 %C. Se emplea un 
venturímetro para controlar el caudal de gas que circula por la conducción. Si a la entrada del 
venturímetro la presión estática del gas es de 240 k/M- 1m *manométricos, calcular la diferencia de 
nivel entre las dos ramas de un manómetro diferencial de mercurio conectado al dispositivo 
medidor. 


Datos y notas: 
Diámetro interno del gaseoducto: 12 pulgadas. Diámetro del cuello del venturímetro: 6 pulgadas. 


Presión atmosférica: 1 bar. Cp/Cv = 1.3. Densidad del mercurio: 13600 kg- m7. Coeficiente de 
descarga: 0.98. 


DATOS 
Fluido: Metano (gas) 
Peso molecular metano: Mametano += 16 Pida 
mol 
: Ss kg 
Densidad mercurio: Pug*= 13600 a 
m 
Razón de calores específicos: y:=1.3 
Caudal másico: Qs =7.48 Kg 
Ss 
Temperatura: T¡:=15.5 *C 
Presión atmosférica: Pam i=1 bar=(1-10*) Pa 
54 iia kN 5 
Presión manométrica: P unan += 240 eii (2.4 - 10 ) Pa 
m 
Presión absoluta: P:=Poam + Pan = (3.4.10*) Pa 
Diámetro y área de la conducción: — D,:=12 in=0.305 m A; =D =0.073 m? 
Diámetro y área del venturímetro:  D,:=6 in=0.152 m A, ad] =0.018 m? 
, 6 D, 
Razón de diámetros: pB:=—=0.5 
D) 
Coeficiente de descarga: C:=0.98 
SOLUCIÓN 


El balance de energía mecánica en el venturimetro es: 


2.0,» AP 
O, =C-Y-4)» lid e dee 
=p 


en donde tenemos que calcular previamente la densidad. 
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: 2 P,-M, 
Densidad metano en conducción: Pmetano “= ——222 = 2.267 E: 
R-T, m? 
A A Le —AP 4 
Teniendo en cuenta la fórmula del factor de expansión: Y=1-— . (0.8 +0.8-f8 ) 
YE] 
—AP 2 AP 
O, =C-|1-—.(0.8+0.8-8*)|-4,- pa 
y*P; (1-f*) 
y resolviendo la ecuación: 
AP=-4.303-10* Pa 
de donde puede calcularse la altura manométrica: 
—AP 
AR, += = 0.323 m 
(Pig — Pmetano) "E 
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Problema 4.20 


Un rotámetro posee una escala graduada con 25 divisiones. El flotador empleado está elaborado 
con un material desconocido, y tampoco se conoce su volumen. Por el aparato se hace circular aire, 
su caudal es 0.696 L/s, y se obtiene los siguientes resultados: 


* Temperatura del aire: 33 %C 

* Lectura del rotámetro: 13 

+ Presión barométrica: 1 bar 

. Presión estática: 510 Pa manométricos 

» Peso del flotador en aire: 4.629 

» Peso del flotador en agua: 2.96 g 

+ Diámetro del rotámetro en la parte superior: 1.650 cm 
+ Diámetro del rotámetro en la parte inferior: 1.303 cm 
+ Diámetro del flotador: 1.270 cm 


Calcular el coeficiente Cr del rotámetro, suponiendo que Y = 1 a causa del reducido valor de Ap. 


DATOS 
Fluido: Aire (gas) 
Caudalímetro: Rotámetro 
Peso molecular aire: Maire "=28.9 is 
mol 
Escala rotámetro: H:=25 
Lectura rotámetro: h:=13 
Temperatura: T,:=33 *C 
Presión barométrica: — P,,,:=1 bar=(1-10%) Pa 
Presión estática: Prani=510 Pa P:= Pron + Pray =(1.005 + 10%) Pa 
E L -4 mi 
Caudal volumétrico: Oy:=0.696 2 =(6.96.107*) Y 
Ss Ss 
d A P Maire kg 
Densidad aire: pi= = 1.141 
R-T, mé 
Caudal másico: m:=0y*p= (7.943 . 107%) ES 
S 
Peso flotador en aire:  Peso,;,,:=4.62 gm 
Peso flotador en agua: — PesO.gy:=2.96 gm 
Densidad agua: Pagua*= 1000 HE. 
m? 
El volumen del flotador puede calcularse teniendo en cuenta el empuje en agua: 
Pes0 aire — Peso 
Volumen flotador: bad (Pesos asu0) — 0.002 L 
Pagua 
] PeS0O ire 3 kg 
Densidad flotador: Pfotador'= === (2.783 «107 ) a 
flotador m 
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Diámetro flotador: Dhotador*= 1.27» cm=0.013 m 


Área flotador: Aida Dorados” = (1.267+107*) m? 
Diámetro superior rotámetro: Dsotsup=1.65 cm=0.017 m 


Diámetro inferior rotámetro: Drotinf:= 1.303 cm=0.013 m 


Área inferior rotámetro: Aj= ED 2 =(1.333-10%) m? 


rotinf — 


Teniendo en cuenta consideraciones geométricas, el área del rotámetro a la altura del 


flotador es: E 
D —D rot. ¡ 
do EA Dit rot.sup rot.inf .hl. E (1.728 . 107*) m? 
4 H 
Área del anillo: Acnitlo= Az —Afotador = (4.616 - 107?) m? 
SOLUCIÓN 
Aplicando el balance de energía mecánica al rotámetro: 
2Zeo.p. Ve lpr— 
Me Cr iS TA y sat L es 2) 
Asilo 
Ay =p 
, A) 
de donde puede despejarse el coeficiente de descarga: 
Cp:= dl =0.598 
2 EPS Y otdor Ñ (Protador —P) 
anillo * 
Ág O ' al 5) 
'otador 
A 
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Problema 4.21 


Un rotámetro tiene un tubo de 12 pulgadas de largo, con un diámetro interno de 1 pulgada en la 
parte superior y 3/4 pulgadas en la base. El diámetro del flotador es de 3/4 pulgadas, su densidad 
es 4.8 g- cam? y el volumen es de 0.4 pulgadas cúbicas. Si el coeficiente de descarga es 0.72, ¿a 
qué altura se encontrará el flotador cuando se mida una corriente de agua con un caudal de 90 
galones/h? 


DATOS 
Fluido: Agua 
Medidor: Rotametro 
Densidad agua: Pagua*= 1000 E 
m? 
Longitud rotámetro: L:=12 in=0.305 m 
Diámetro superior rotámetro: — D,,,:=1 in=0.025 m 
Diámetro inferior rotámetro: Dar== in=0.019 m 
Área en la entrada: A,:=%.D,)¡ =(2.85-107*) m? 
4 
Densidad flotador: Prorador:=4.8 E%-= (4.8 +10?) Hg 
cm? m? 
Diámetro flotador rotámetro:  Daorador 7 in=0.019 m 
Área flotador: Añotador E “Dhotador =(2.85+ 107) m? 
Volumen flotador: — Vario :=0.4 in? =0.007 L 
Coeficiente de descarga: C:=0.72 
Factor de expansión: Y:=1 (por ser un líquido) 
3 
Caudal volumétrico: Qy:=90 > (9.464.107) Y 
Y Ss 
En k, 
Caudal másico: m:=Qy* Pagua =0.095 e 
SOLUCIÓN 


Para un rotámetro, el caudal másico viene dado por la siguiente expresión, donde el subíndice 1 
hace referencia a la entrada al rotámetro y el subíndice 2 al anillo entre el flotador y el tubo del 
rotámetro: 
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Sustituyendo los valores conocidos y resolviendo la ecuación: 


A... =9.96:10 7% m? 


anillo 


Conociendo el área del anillo y mediante consideraciones geométricas, podemos conocer la altura a 
la que se obtendrá dicho anillo: 


Area del tubo a la altura buscada: Asor*=Aanillo + Áfotador = (3.847 +10) m? 
.Z 4. Aro 
Diámetro del tubo a la altura buscada: — D,,,:= =0.022 m 
T 
Ei 0z nos Das Di 
Variación del diámetro con la altura en un tubo troncocónico: Dijo =D ¡+ h 
L 
; Dro — Dinp) *L 
Despejando: h:= (Dor —Ding) L — 0.148 m 
o a Dinf 
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Problema 4.22 


En una planta se bombea esencia de trementina desde la base de una columna de destilación hasta 
un gran tanque de almacenamiento con un caudal de 90 L/min. El nivel de líquido en el tanque está 
5 m por encima del nivel de la base de la columna. El tanque está abierto a la atmósfera, mientras 
que la columna opera a 0.4 atm de presión absoluta. 

La columna y el tanque se encuentran unidos por una conducción de tal forma que la pérdida de 
energía mecánica que experimenta el fluido a través del sistema es de 330 J/kg. 

Calcular la energía necesaria para impulsar el fluido mediante una bomba (en watios). 


Datos y notas 

Densidad = 0.86 g/ cn . Presión atmosférica = 1 atm. Dado el gran tamaño de la columna y del 
tanque de almacenamiento, considerar que la velocidad del fluido es despreciable en la base de la 
columna y en la superficie del tanque. 


DATOS 
Fluido: Esencia de trementina (líquido) 
Densidad del fluido: p:=0.860 E7_=860 € 
cm m? 
Diferencia de alturas: — 4z:=5m 
Presión de entrada: P,:=0.4 atm=(4.053 - 10*) Pa 
Presión de salida: P,:=1 atm=(1.013+10*) Pa 
Aa L m* 
Caudal volumétrico: O,:=90 —-=0.002 — 
min Ss 
Caudal másico: m:=Q,+p=1.29 E 
Ss 
Pérdidas de energía mecánica: e,:=330 - 
E 
SOLUCIÓN 


El balance de energía mecánica es: 


despejando, el trabajo específico será: 


PE J 
Wesp*= Quel) +2g+4z+e,=449.725 — 
p kg 
y teniendo en cuenta el caudal másico: 
W:=w.,) m=580.145 W 
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Problema 4.23 


Se dispone de un orificio medidor que tiene un diámetro de 0.75", situado en una conducción de 3" 
de diámetro interno. A través de la conducción circula agua a 15%C, La diferencia de presión se mide 
con un manómetro diferencial de mercurio bajo agua. Calcular el caudal que circula cuando la 
lectura manométrica es de 50 cm. 


Datos y notas 
Densidad del agua: 1 kg/L, Viscosidad del agua: 1 cP, Densidad del mercurio: 13.650 g/mL. 


DATOS 
Fluido: Agua (líquido) 
Medidor: Orificio 
Densidad agua: Pagua*= 1000 E: 
m? 
. . . 3 kg 
Viscosidad agua: 1:=10 7" Pa+s=0.001 
mes 

Temperatura: T:=15 *C 
Diámetro conducción: —D,:=3 in=0.076 m 
Diámetro orificio: D,:=0.75 in=0.019 m 

Área orificio: A) =D? = (2.85-107*) m? 

z “. D, 
Razón de diámetros: fp:=—=0.25 
1 
Manómetro: Diferencial de mercurio bajo agua 
Densidad mercurio: Pug"= 13650 pa 
me 

Altura manoméftrica: Ah,,:=50 cm=0.5 m 

Diferencia presiones: AP:=—Prig—Pagua) E * Am =—6.203 +» 10% Pa 

SOLUCIÓN 


Suponiendo que el régimen es turbulento, el coeficiente de descarga sería: C:=0.6 


Aplicando el balance de energía mecánica, y teniendo en cuenta que circula un líquido (Y=1): 


a 
m=C.4 | PP 1908 E 
=p 5 


3 
=0.002 Y 
AS) 


Pio z m 
el caudal volumétrico será: 


Oy:= 


Pagua 


Ahora debe comprobarse que el régimen de circulación en el orificio es, efectivamente, turbulento: 


.v..D 
E O Re, := Legua" Y2*22 1 976.10* Re,>10* 
2 A 2 2 


2 s 17) 
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Problema 4.24 


Se mide el caudal de circulación de agua en una tubería de 6" de diámetro interno, con un 
venturímetro que tiene la garganta de 2" de diámetro. Cuando la pérdida de presión entre la 
sección (1) y la garganta (2) del venturímetro es de 10.7 cm. de agua, el caudal de circulación es 
de 2.9 kg/s. ¿Cuál es el coeficiente de descarga del venturímetro para ese caudal y cuál es la 
pérdida de energía mecánica debida al rozamiento hasta la garganta?. Si la diferencia de 
presión entre la sección (1) y la sección (3) es de 1.52 cm. de agua ¿Cuál es la pérdida de presión 
permanente expresada como porcentaje de la diferencia de presión leída entre la sección anterior 


(1) y la garganta (2)? 


=-—1.049.10* Pa 


DATOS 
Fluido: Agua (líquido) 
Medidor: Venturímetro 
Diámetro conducción: — D,:=6 in=0.152 m 
Área conducción: Ap= DP =0.018 m? 
Diámetro estrechamiento: D,:=2 in=0.051 m 
Área estrechamiento: A ES =0.002 m? 
Razón diámetros: e 0.333 
D; 
Densidad agua: Pagua"= 1000 Ha 
m 
Caudal másico: m:=2.9 E 
Manómetro: Diferencial de agua 
Altura manométrica: Ah,,:=10.7 cm 
Diferencia presión: AP ¡3:=—Pagua 
SOLUCIÓN 


Aplicando el balance de energía mecánica, y teniendo en cuenta que circula un líquido (Y = 1): 


—2+ + AP 
O, pa C-A) a Pegue” 12 
5 
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Despejando el coeficiente de descarga: 


> *Pagua* AP ¡> 
1-P 


=0.982 


C:i=m= Á)* 


Para un venturímetro, la pérdida de presión permanente es, aproximadamente: 
AP moy =0.1 + AP, =—104.931 Pa 


perm* 


También puede calcularse la pérdida de presión permanente a partir del manómetro conectado 
entre las secciones 1 y 3: 


Ah,,3:=1.52 cm AP ¡5:=—Pagua* 2 * Ay 13 ==149.061 Pa 


que es muy próxima a la estimada por la fórmula general. 


La pérdida de presión permanente, expresada en la forma solicitada por el problema es: 


AP ¡3 


+ 100=14.206 
Pr 
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Problema 4.25 


Un orificio medidor se calibra para medir el caudal de una corriente de acetona líquida. Se utiliza un 
manómetro diferencial simple con un fluido manométrico de densidad 1100 kg/ m7. Los datos 
obtenidos experimentalmente son los siguientes: 


Caudal Lectura manométrica 
Q (mL/s) Ah (mm) 


(a) Para cada una de las lecturas, calcular la caída de presión a través del orificio, (-Ap) en 
mm de Hg. 

(b) El caudal a través de un orificio está relacionado con la caída de presión a través del 
mismo por la expresión: QO=XK-(-4p)” siendo: Q = ml/s y (Ap) = mm Hg. 
Determinar los valores de K y n que mejor se ajusten a los datos. 

(c) Si la lectura del manómetro en cuestión fuese Ah=17 mm, determinar el caudal d 
circulación de la acetona. 


Datos: 
Densidad de la acetona: 791 kg/ mM. 


DATOS 
Fluido: Acetona 
Manómetro: Diferencial con acetona y un líquido manométrico 
Densidad acetona: E 
m? 
Densidad fluido manométrico: Pan = 1100 Ea 
m 
: e kg 
Densidad mercurio: Phg:=13650 — 
m 
SOLUCIÓN 


(a) Para cada lectura, hay que calcular la diferencia de presiones (AP) y expresarla como altura 
de una columna de mercurio: 


AP (4h,,) == =( nun =P acciona) -2-4h,, AP= —Phg*8* Ah yg 
—AP 
0 a E AP(S mm)=-15.2 Pa di PA mm 
Ss PHg "E 
AP, 10 mm 
Om =P. AP, (10 mm) =-30.3 Pa Ahyg= al 
Ss PHg ES 
—AP3 (15 mm 
Oy :=424 mE AP, (15 mm) =-45.5 Pa diia 200 0.34 mm 
s Pug* 8 
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AP, (20 
Oys:=490 HE AP, (20 mm) =-60.6 Pa Aj AP. 0.45 mm 
s Phg*8 
AP, (25 
Op YE AP, (25 mm) =-75.8 Pa PM ad ata 0.57 mm 
Ss PHe*8£ 
AP, (30 
Oys:=600 HE AP¿(30 mm) =-90.9 Pa iy OI 0.68 mm 
s PHg*8E 


(b) Para calcular el valor de las constante de la ecuación, se toman logaritmos y se hace el 
ajuste por mínimos cuadrados de los datos: 


log(0) = log (K) +n- log (AP) 


Obteniéndose: EETDTAS n:=0.5 


(c) Para la lectura indicada, usando las mismas ecuaciones: 


—AP. 
Ah, :=17 mm AP.o:=—(Pman — Pacetona) * 8 * Ame ==51.514 Pa Ahygo'= “ =0.385 mm 


PhHg :g 
mL 
Oy.:=K+0.385" =451.588 Y 
Ss 
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Problema 4.26 


Por una conducción de 3" estándar de diámetro (diámetro interno = 3.068") fluye aire a una presión 
manométrica de 10.2 atm y 27 “C. En la conducción hay instalado un orificio medidor de 3.5 cm de 
diámetro. La presión del aire tras el orificio es de 9 atm manométricas. Considerando el aire como 
un gas ideal, calcular el caudal másico de aire que circula por la conducción, en kg/h. 


Datos y notas: 


Cp = 7 kcal/(kmol*K); Cv = 5 kcal/(kmol*K); u= 0.018 cP; Patm = 1 atm. 


DATOS 
Fluido: Aire (gas) 
Caudalímetro: Orificio 
Peso molecular aire: PM:=28.9 kg 
mol 
Viscosidad aire: y::=0.00018 poise=(1.8.107%) 4E_ 
mes 
Presión atmosférica:  P,,,:=1 atm=(1.013-10%) Pa 
Temperatura: T,:=27 *C 
de E , kcal 
Calor especifico a presión constante: cy:=7 
kmol + K 
Calor específico a volumen constante: Cy:=5 a 
kmol + K 
Razón de calores específicos. K=TE=14 
Cy 
Presión conducción: P,:=10.2 atm+ Pa =(1.135+10%) Pa 
Presión orificio: P,:=9 atm+Pa,, =(1.013 + 10%) Pa 
Diferencia de presión: AP:=P,-—P,=-1.216+10? Pa 
Diámetro conducción: — D,:=3.068 in=0.078 m 
Diámetro orificio: D,:=3.5 cm=0.035 m 
] e D, 
Razón de diámetros: f:=—=0.449 
D; 
Área orificio: Ar=%.p? = (9.621 107+) m? 
4 
SOLUCIÓN 


El balance de energía mecánica aplicado al orificio es: 


—2.+p,- AP 
m=C+Y+- A,» MEA ic 
1-8? 
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Previamente, hay que calcular: 


z P,+PM 
Densidad del aire en la conducción: p=— =13.14 kg 
R e T, m? 
o AP 4 
Factor de expansión: Yi=1l4 - (0.41 +0.35+8*)=0.968 
E 7 
Para orificios, suponiendo régimen turbulento: C:=0.6 
Con lo que el caudal másico es: 
-2+p,- AP 
m:«<=C+Y+-A)- EP" 3670 E 
15 hr 


Ahora es necesario comprobar si el régimen de circulación en el orificio es realmente turbulento: 


: , SEL P,+PM 
Densidad aire en orificio: P> i=2 2 =11.73 E 
R > T,; m? 
q io m m* 
Caudal volumétrico en orificio: Oy,:=—=0.087 —=— 
P2 s 
Velocidad aire en orificio: v> = 27 2903 EE 
Á) s 
ies .v,*D 
Reynolds en orificio: Ras AE Re>10* 
17) 
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Problema 4.27 


Para medir el caudal de circulación de un fluido incompresible que circula turbulentamente a 
través de un conducción, se instala en la misma un venturímetro cuyas dimensiones son tales que el 
area de la sección transversal de la conducción (1) es cuatro veces la de la sección de la garganta 
del venturímetro (2). 


Venturímetro 
— 


(a) Encuentra una relación entre las velocidades ví y vz en las secciones (1) y (2). 
(b) Escribe la ecuación del balance de energía mecánica para el sistema entre los puntos (1) y 
(2) y utilízala para demostrar que, cuando se desprecian las pérdidas de energía por 


rozamiento. 

(c) Calcular el caudal másico que circula por el sistema anteriormente descrito, si el fluido que 
circula es agua (p = 1000 kg/m3), el fluido manométrico es mercurio (pm= 13650 kg/m3), 
el diámetro de la garganta del venturímetro es de 2 cm y la lectura manométrica es de 


9 cm. 


D,=4+D, D,==— 


Fluido incompresible: 0a=05 


SOLUCIÓN 


(a) Las velocidades correspondientes en cada sección son: 


On _ Qv _0v_ Qn 
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y con el dato del problema: 


(b)  Planteando el balance de energía mecánica, sin pérdidas por rozamiento: 


P,-P y? y? 
A 
1) 2 


Suponiendo a1=a2=1 , reordenando y simplificando: 


2. (_ % Y_(_Q YP]_2040-0/' -p _ 2040-m? 


O 2 AS 2 
217 (D, Tp? 16+S) 16-p+S) 
414 4 
(c) Utilizando los datos del apartado (b): 
Densidad del agua: Pa:=1000 KE 
m? 

a 5% kg 

Densidad del mercurio: Phig'= 13650 
mE 

Altura manométrica: h,,:=9 cm 


Diferencia de presiones:  4P:=—(pjg—Pa) *2* Mm ==1.116+10* Pa 
Diámetro venturímetro: D),:=2 cm 


Diámetro conducción: D,:=4-D,=8 cm $=+ D? =0.005 m? 


—AP«16»p,.S,? 
m:= A Kg 
2040 s 


el caudal másico es: 
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Problema 4.28 


A través de una tubería de 150 m de longitud y 3 pulgadas de diámetro interno fluye agua con un 


caudal de 14.5 L/s. Al final 


de este tramo la tubería se ensancha hasta 4" de diámetro interno. 


Calcular la pérdida de energía mecánica de la conducción. 


Datos y notas: 


Densidad del agua: 1 kg/L. Viscosidad del agua: 1 cP. 


(2) 
] 
l=7- P,=2bar 
1 
] 
1 
| 
Az = 30 m 
(1) 
|] 
O_O e 
I L = 150 m 
P, = 9 bar 
DATOS 
Fluido: Agua 
Viscosidad agua: 1:=0.01 poise=0.001 8 
mes 
Densidad agua: pi=1000 E 
m? 
Presión entrada: P,:=9 bar=(9.10*) Pa 
Presión salida: P,:=2 bar=(2+10*) Pa 
Diferencia de alturas: — 4z:=30m 
Longitud conducción:  longitud:=150 m 
oz az j a+D¡ 2 
Diámetro y sección entrada: D,:=3 in=0.076 m S,:= a =0.005 m 
bd ez ; . a+ D, 2 
Diámetro y sección salida: D,:=4 in=0.102 m S):= a = 0.008 m 
As L m' 
Caudal volumétrico: O, :=14.5 2=0.015 
Ss Ss 
SOLUCIÓN 
El balance de energía mecánica es el siguiente: 
P,-—P 2 2 
E Lg Adz+ . AMES le =-£, 
P 2 14 01 
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Para calcular los términos de la energía cinética, debemos obtener previamente la velocidad y el 


régimen de circulación: 


Velocidad entrada: 


Reynolds entrada: 


Velocidad salida: 


Reynolds salida: 


O, 
S; 
Re) = dd 
_ Or 
S, 
Re, == 


y despejando, la pérdida de energía mecánica es: 


e, == 
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P,—P 1 [v 
22 + go 47+ . 
P 2 


8 longitud S 


2 


ea 
2 Y1%1 2.423-10* 


o 
20212 1817.10* 


A 


E, :=8,+p+Q,=5934 W 


90 
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39.6 — 


Re,>10* 0. :=1 


Re,> 10* ay:=1 
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BLOQUE Il! 


Balances de materia 
en sistemas con reacción 
quimica. Reactores Ideales 


Problema 5.1 


El antimonio se obtiene calentando sulfuro de antimonio, pulverizado previamente, mezclado con 
chatarra de hierro y drenando más tarde el antimonio por el fondo del reactor. La reacción que tiene 
lugar viene representada por el esquema: 


Sb2S3 + 3Fe 5 25Sb + 3FeS 


Si en una instalación de laboratorio (intermitente como el reactor descrito) se mezclan 0.6 kg de 
Sb2S3 y 0.25 kg de limaduras de hierro y se calientan, se obtienen 0.2 kg de Sb, calcular: 


(a) Cuál es el reactivo limitante. 

(b) El porcentaje de exceso de reactante que se encuentra en proporción superior a la 
estequiométricamente necesaria. 

(c) El grado de conversión alcanzado (referido al reactivo limitante). 

(d) El rendimiento. 


Renombrando los compuestos químicos: A, =Sb,S, A,=Fe Az=Sb A¿=FeS 
el esquema de reacción es: —A¡ —-3:4,+2:-A3+3-A,¿=0 
donde los coeficientes estequiométricos son: v¡:=-1 Vyi=-3  13:=2 v¿:=3 


Los pesos moleculares de cada especie: 


k k k 
PM,:=339.72 9 PM,:=55.85 PM,:=121.76 9 PM,:=87.92 9 
mo mol kmo mol 
Las cantidades iniciales en el sistema son: 
M,0:=0.6 kg Noi=24% 1.766 mol 
10 > 10 PM, , 
M»o:=0.25 kg Nao '= Moo — 4.476 mol 
PM, 
M3p:=0 kg Nz¿:=0 mol 
y la cantidad de producto: 
M; 
M3:=0.2 kg N3:= = 1.643 mol 
PM; 
SOLUCIÓN 
(a) Reactivo limitante: 
— | min 
v; 
N NÑN. 
19 | 1.766 mol 22 |—1.492 mol 
vi Va 
El reactivo limitante es el hierro 
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(b) Exceso de reactivo: 


vj 
Nio—Nro* PR 
Exceso_Aj= E 
Nr 
kO 
Vy 
Y 
Ni No" 
va 
Exceso_A; := =0.184 
vr 
Nao" 
va 


(c) Grado de conversión: 


A partir de la relación estequiométrica:: N,=N,¡ Ny: X- 


vr 
“ z va 
sustituyendo y despejando: N3=N3 Na X- 
vz 
N. . 
30 3 
ES =0.55 
V. 
3 
Nao ú 
vez 
(d) Rendimiento: 
moles_de_Aj_formados |1x 
7 moles de _Ak iniciales v; 
N3-N 
3 30 2 
B,:= .||=0.55 
N 
20 43 
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Problema 5.2 


En un reactor se produce la siguiente reacción: Ai+A2-A3=0 
Si la composición de la mezcla inicial es tal que Nio = 1 mol, Nz2o= 2 mol y N3o = 5 mol, indicar qué 
conjunto de valores de Nj podrían obtenerse por la reacción indicada. 


(a) Ni = 2 mol N2 = 2.5 mol N3 = 4 mol 

(b) Ni = 3 mol N2 = 4 mol N3 = 1 mol 

(c) Ni = 1.5 mol N2 = 2.5 mol N3 = 4.5 mol 

(d) Ni = 0 mol N2 = 1 mol N3 = 7 mol 
DATOS 


El componente clave es A3, ya que es el único reactivo. 
Los coeficientes estequiométricos son: v:=1 v):=1 v3:=—1 
Las cantidades iniciales son: Ni¡p:=1 mol N»py:=2 mol N3p:=5 mol 


En todos los casos tiene que cumplirse la relación estequiométrica: 


v; Nijo—N; 
N¡=Nj Ny X- despejando X= == 
Vr v; 
Nos 
Vr 
SOLUCIÓN 
(a) 
N,,.:=2 mol Na :=2.5 mol Nz, :=4 mol 
Nio-N;1 Nao —N» Nao N3 
La i= ” =0.2 Xq = ”=0.1 Xy, =0.2 
v D, LES 
N3o HER N3o El Ny Al 
vz Vz 3 
Imposible, el grado de conversión es distinto. 
(b) 
Ni¡,:=3 mol N»,:=4 mol N3,:=1 mol 
Ni =N Na =N: N3y=N 
ben Pl A XK ae IS pe y O 
vi v) 43 
N3o* N30*— N3o* 
vz vz vz 


Imposible, el grado de conversión es distinto. 
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(c) 


N,.:=1.5 mol N,,:=2.5 mol N3,:=4.5 mol 
Niy=Nj Na —N, N3-Nz 
X= -==0.1 X= -=0.1 A ==0.1 
Dv; V) DA 
N3p* N3o* N3o* 
vz vz vz 
Posible, el grado de conversión es el mismo. 
(d) 
N;¡¿:=0 mol N>¿:=1 mol N3¿:=7 mol 
Nio-Nia Na —Noa Ns —N3a 
X14 == -0.2 X24 == -0.2 X34 *= ————— = 0.4 
vs V, Vs 
Nyo* N3o* N3o* 
vz vz vz 
Imposible, no puede haber un grado de conversión negativo. 
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Problema 5.3 


Una mezcla líquida miscible, cuya composición inicial es co, reacciona a volumen constante (V1) 
hasta alcanzar un grado de conversión X. En ese momento se congela la reacción y se le añade un 
volumen (V2) de una mezcla con concentración c'í (en cada componente presente). Deducir la 
ecuación matemática que permite calcular la composición de la mezcla resultante (c”;), expresada 
como concentración molar de cada especie química presente. 


SOLUCIÓN 
La concentración de todas las especies cuando se congela la reacción es: 


DD: 


C¡=Cio Cro X+* 
V; 


El número de moles, Nj, procedentes de la mezcla inicial es: 


A 
N,=c,"V,= Cjo — Egg A 


-V, 
Ve 


El número de moles añadido es: 


¡SEN os 
N'y=0 ¿"Va 


El número total de moles en la mezcla es la suma de los moles aportados por la propia reacción y 
por la disolución añadida: 


v, 
/ J / 
N¡+N'¡= TA «Vi +cyVa 
k 


y la concentración será teniendo en cuenta el volumen total: 


D: 
Cjo—Cro*X- = «V +0" Va 
k 


” 


do VE 
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Problema 5.4 
En un reactor de flujo de pistón (RFP) se está desarrollando la reacción irreversible en fase gas: 


A> 3B 


El alimento es A puro, con una concentración inicial de 20 kmol/m?* . En la tabla adjunta se detalla 
el grado de conversión que se ha alcanzado a diferentes distancias desde la entrada del reactor. 
Calcular la concentración de A en cada una de estas secciones: 


% de L 0 12.5 25 50 75 100 
X 0 0.1 0.2 0.3 0.4 0.5 


siendo L la longitud total del reactor. 


DATOS 
El esquema de reacción es: —-1:-A+3-B=0 
y los coeficientes estequiométricos: vAi=-1 Ug:i=3 
) a kmol kmol 
Las concentraciones iniciales son: Ca0i=20 — Cgo:=0 —= 
m 
Como el alimento es A puro: Lap:=1 


Como no se indica nada, consideramos que el reactor trabaja a la misma temperatura y presión que 


la corriente de entrada: 
P=P, T=To 


El componente clave es A, por ser el único reactivo. 


SOLUCIÓN 


Se trata de una reacción en fase gas, con variación en el número de moles, por lo que la 
concentración variará con el grado de conversión según la siguiente expresión: 


Ne 
Cko* 0, XA. —= 


v, A 
e _ Ej = Tigo: 0 q e 
Py T Uy 
» (1 + Ex -X) 
PR 
Sustituyendo: 
P,*T A 
sl =1 Av:=1V¿+Up=2 Ó:= sl =2 EJ i= Tap 0=2 
P. To —VA 
CAo iS (1 —X) 
de (X) = 1+€,¿:X 
A” 
kmol kmol kmol 
ca(0)=20 == (e E ajena LE 
m m m 
kmol kmol kmol 
ca (03) 88 0 e) — 
m m m 
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Problema 5.5 


En un reactor químico continuo (sistema dinámico), se desarrollan simultáneamente las reacciones 
siguientes: 

3H2 + N2 S 2NH3 (1) 

H2 + CO2 S CO + H20 (2) 


si la mezcla reactante contiene el 68.4% de Ha, el 22.6% de N2 y el 9% de CO», y la corriente 
gaseosa a la salida del reactor contiene un 15% de NH3 y un 5% de H20 (en forma de vapor de 
agua en las condiciones de salida): 
(a) ¿Qué especie química es aconsejable elegir como componente clave, al que se debe 
referir el grado de conversión? 
(b) Calcular los porcentajes de H2 y N2 en la corriente que sale del reactor. 
(c) ¿Cuánto vale la conversión alcanzada por el Hz en la primera y en la segunda reacción 
consideradas? 


Notas: 


Todos los porcentajes son molares. Para realizar los cálculos supóngase que la mezcla reactante se 
introduce en el reactor con un caudal de 100 kmol/día. 


Asignaciones y coeficientes: 


A,=H, A,=N, A¿=C0, A¿=NH, A¿=C0O0 A¿=H,0 
V11:=-3 Vi9:=—1 Vi3:=0 Vi4:=2 Vi5:=0 Vi6:=0 
Va :=-1 Va :=0 Vo:=—1 Vo =0 Vosi=1 Voi=1 
Reacción 1: —3 A¡-4,+2 A4,=0 Ap :=0 ¡+2 a =-2 
Reacción 2: —A¡ —Azg+A5+A5¿=0 Apisi=0 
; kmol 
Base de cálculo: Fio:=100 da 
day 


Las fracciones molares de entrada y salida conocidas son las siguientes, de las que calculamos 
también los caudales de dichos compuestos: 


kmol 
2,0:=0.684 Fipi=Fioo 2 =68.4 Y 
day 
kmol 
La) :=0.226 Faoi=Fyo*La0=22.6 . 
kmol 
Z3p:=0.09 Fzy:=F9* L3p=9 Ta 
day 
x,¿:=0.15 
25:=0.05 


SOLUCIÓN 


(a) El componente clave es la especie química Ai, puesto que interviene simultáneamente en las dos 
reacciones como reactante. 


(b) Para determinar la composición de la corriente de salida del reactor hay que plantear las 
correspondientes ecuaciones estequiométricas para cada especie: 
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R 
para el componente clave: F,=Fso* l - Sx, 
i=1 


R " 
para cualquier componente: Fj=Fj—Fro* z |. 2) F¡=F,-2; 
i= ik 
y en la corriente de salida se cumplirá: 
Ss Ss Ss R Vi R Av, 
F¿= Y) Fj= Y Fio—Fro* ») »|X;" =Fo—Fror )|X;* 
j=1 j=1 j=1 i=1 Vi i Vik 


y sustituyendo se tiene el siguiente conjunto de ecuaciones: 


F,=Fp*(1-X,-—X;) 


Via 
F3=F29—|Fio*X;* 


Y11 
Ya3 
Et Pt 
Va] 
Fy=-—|Fi:X;* =F,* Lg 
11 


21 


Va 
TS | 


Va 
F¿=-—|Fio*X>* =F,*L6 
21 
Av; 
F¿=Fy —Fio*X;* 
V11 


que puede resolverse secuencialmente. Sustituyendo (7) en (4): 


v Av PT 
AS A EA A A X,:= a =0.286 
V11 Vi1 Av, Via 
Z4* Fi =Fig-= 
: o 11 vi 
sustituyendo (7), (6) y siguientes: 
Av kmol Fx 
Fi=Fo=Fio Xi ——=86.997 LE X,:= 2% — 0.064 
Vii ay Va 
—|F¡.— 
Var 
kmol v kmol 
Py :=Fip* (1—X,—Xy) =44.487 == [A OS A E Pi 
day 44 day 
a kmol kmol 
Fy=Fay | Fyp+ Xy: |=4.652 HO Fy=F¿+2,=13.043 42 
Va1 day day 
v kmol kmol 
Fa==|F «A |=a348 HS Fo:=F,-2.=4.348 22 
Var day day 
con estos resultados ya se pueden calcular las fracciones molares de la corriente de salida: 
F; F, 
(0) 2, :=—=51.16% 2, :=—=18.49% 
F, t 
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Problema 5.6 


En la fabricación del ácido nítrico se oxida el NH3 con aire para producir NO, que más tarde se oxida 
a HNOz3. El NH3 procede de un gasómetro y tiene una pureza del 100 %. El aire proviene de un 
compresor que lo suministra a la presión atmosférica, con un contenido en agua de 1.4 moles de 
H20/100 moles de aire seco. 

Admítase que en el reactor tienen lugar solamente las reacciones: 


2NH3 + 2.502 S 2NO + 3H20 
2NH3 + 1.502 5 N2 + 3H20 


y el 98 % del NH3 reacciona según el primer esquema de reacción y el 2 % restante según el 
segundo. 

Si la mezcla inicial contiene 11 1 de NH3 por cien de mezcla gaseosa exenta de agua, calcular la 
composición del gas a la salida del reactor expresada como fracción molar y expresada en kmol de 
un componente cualquiera por cada kmol de NH3 introducido en el sistema. 


Datos y notas: 
Admítase que los volúmenes molares de todas las especies químicas tienen el mismo valor. 
Refiéranse los cálculos a un caudal molar de mezcla reaccionante (exenta de agua) de 100 kmol/s. 


DATOS Y CÁLCULOS PREVIOS 


Asignaciones A¡=NH; A,=0, A¿=NO A¿=N, A; =H,) 
i=1 v11:=-2 11) :=-—2.5 Vi¡3:=2 DV¡4=0 Vi5:=3 
1=2 Vo] :=-2 Vo) i=—1.5 Voz :=0 Vo¿:=1 Va5:=3 

, kmol 

Base de cálculo: F,:=100 da 


Ss 


En la mezcla inicial, suponiendo comportamiento ideal, la relación en volumen es la misma que la 
relación molar. Entonces, para la mezcla exenta de agua, teniendo en cuenta la composición del 
aire, resulta: 


kmol 


£1p:=0.11 Fjo:=Fo*Ljp=11 


Ss 
kmol 


L7p :=0.89+0.21=0.187 Fap:=Fp* yy =18.69 


Ss 
kmol 


T 4) :=0.89-0.79=0.703 Fyo:=Fp* 249 =70.31 


Ss 
kmol 


Tgy:=0 F3p:=Fy* L3y=0 


el caudal molar de agua viene dado por la relación del enunciado: 


1.4 kmol 
Fsy:= (Fao +Fso) * =1.246 
100 Ss 
; kmol 
Por lo que el caudal molar total será: Pi:=Fi0+Fo+Fia+F»s9= 101.246 
Ss 
La conversión por cada reacción es: X,:=0.98 X,:=0.02 
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SOLUCIÓN 


Se plantean las relaciones estequiométricas para un sistema dinámico con dos reacciones: 


F¡:=F,p* (1-X,-—X,)=0 


Fy:=Fy— 
U13 

F3:=Fyy=| Fi: X;* = Fay a 

Pr=Fa-|Porx : 

Ea=Po=|Po-x , 


F :=F,+F)4+F3+F,+F;5=103.996 


Fi: X . 


kmol 


v 

z 5.0 kmol 
V) Ss 
v 

a y pl 

2 
v 

») —70.42 kmol 
V, Ss 

kmol 
V) Ss 
kmol 
Ss 


de donde pueden calcularse las fracciones y relaciones molares pedidas: 


F; 
Li :=—=0 

F, 
2,:=2=0.049 
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F; 
n; = —=0 
Fi 
n,:==0.459 
0 
F 
n3:=—=0.98 
0 
F 
ny:=—= 6.402 
0 
n5i=—=1.613 
Fi 
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Problema 5.7 


El NO se obtiene por combustión de amoníaco gaseoso con un exceso molar de oxigeno del 20 % 
sobre el estequiométricamente necesario, según la reacción. 


4NH3 + 502 > 4NO + 6H20 


En las condiciones de operación el grado de conversión del amoníaco es del 70 %. De la corriente 
producto se separa el amoníaco para reciclarlo de acuerdo con el esquema adjunto. 


REACTOR SEPARADOR 


Calcular: 
(a) Los moles de NO formados por cada 100 moles de amoníaco introducidos en la corriente de 
alimento fresco. 
(b) Los moles de amoníaco recirculados por cada mol de NO formado. 


En primer lugar se hacen las asignaciones correspondientes a las corrientes y a las especies 
químicas. 
Corriente de entrada al sistema: O Amoniaco: Ai 
Corriente de entrada al reactor: a Oxígeno: A>2 
Corriente de salida del reactor: b Monóxido nitrógeno: A3 
Corriente de salida del separador: Cc Agua: As 
Corriente de amoniaco recirculado: d 


La ecuación puede reescribirse como: —4-A¡—5-4,+4-A43+6-4,¿=0 
Siendo los respectivos coeficientes estequiométricos: — v,:=-—4 Vy:=-5 v3i=4  1v¿:=6 
Según el enunciado, el componente clave es el amoniaco, mientras que el oxígeno está en un 


exceso del 20 %. 
El grado de conversión es: X:=0.7 


SOLUCIÓN 


Para conocer los caudales molares de la corriente de entrada al reactor procedemos de la forma 
siguiente. En primer lugar, vamos a definir una base de cálculo, a partir de la corriente de amoniaco 
que entra al reactor: 


l 
FP. =100 2 


Ss 


. An p va 
El oxígeno estequiométricamente necesario es: Fi," 


y; 


por lo que si usamos un 20 % mas de dicha cantidad, el caudal de oxígeno a la entrada del reactor 


será: 
va mol 
Faa:=(1+0.2)+[F,,-—Ñ|=150 
v; Ss 
mol mol 
A la entrada del reactor no tenemos productos, luego: F,:=0 Fia=0 
Ss Ss 
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Con lo que ya conocemos completamente la corriente de entrada al reactor. 


Para conocer los caudales de salida del reactor, aplicamos las relaciones estequiométricas: 


Fyy:=F,,* (1-X)=30 


Fa :=Fz Fra 
Fa=Psa— Pra 
PP Pia 


14 
Xx. 


db 
WE 
3 
Xo 


XX 


mol 


Ss 


0 


Ss 


70 mol 


En las corrientes c y d, los caudales molares de cada especie son los mismos que en la corriente b: 


Fo :=Fyy Fy.:=Fyy 
En el punto de mezcla inicial, se cumplirá que: 


mol 


siendo: PFa0:=Foga=150 


Ss 


Calculándose la corriente de entrada al sistema: 


mol 
Fac :=F ay Fy¿:=F1,=30 A 
Pro =Pis+FPig 
F30'=F3a Fio:=Fsa 
mol 
Fio:=F1¿—Fig=70 a 


Los moles de NO formados por cada 100 moles de amoniaco fresco: 


(a) 


lo 


F 
A 


Los moles de NH3 recirculados por mol de NO formado: 


(b) 


Fa 
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Problema 5.8 


En un reactor discontinuo se planifica desarrollar en fase líquida la transformación de un 
componente A en otro B. La reacción que tiene lugar es irreversible y puede representarse por el 


esquema: 
A>B 


Calcular el tiempo de reacción necesario para que en cada carga la concentración molar del 
componente A descienda desde Cao = 1.3 kmol/ nr aca=0.3 kmol/ nr a la salida. 


Datos y notas: 


oz E ee . 3 
La relación entre la velocidad de reacción, expresada como kmol de A reaccionados por 7 y 


minuto, y la concentración molar de A, se detalla en la tabla siguiente: 


CA o.1 | 0.2 | 0.3 | 0.4 | 0.5 | 0.6 | 0.7 | 0.8 
(Kmol/ 177 ) 


1.0 1,3 2.0 


-RA = -va"r [| 0.1 0.3 0.5 0.6 0.5 | 0.25 | 0.10 | 0.06 


Admiítase que el reactor se comporta como un RDTA ideal. 


DATOS 
Todos los cálculos se refieren al compuesto A (componente clave): 
E pi e A 
mi m? 
Tabla de datos anterior. 
SOLUCIÓN 
a d 
El balance de materia en un RDTA es: —Ca= Ya r=R, 
di 
separando variables e integrando por trapecios: t= (7, dc, 
Á 
Cao 
Integral en 
Ca Ra 1/Ra intervalo 
kmol/m? | kmol/mé-min | m-min/kmol min 
0.1 -0.1 -10.0 
0.2 -0.3 -3.3 
0.3 -0.5 -2.0 
0.4 -0.6 -1.7 0.18 
0.5 -0.5 2.0 0.18 
0.6 -0.25 -4.0 0.30 
0.7 -0,1 -10.0 0.70 
0.8 -0.06 -16.7 1.33 
1 -0.05 -20.0 3.67 
1.3 -0.045 -22.2 6.33 
2 -0.042 -23.8 
SUMA 12.70 
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0.05 | 0.045 | 0.042 


t:=12.7 min 
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1/R, (m3/(kmolmin)) 


' 
N 
un 
o 


0,0 0,2 04 0,6 0,8 10 12 14 16 18 2,0 
Ca (Kmol/m*) 


NOTA: Existen diversos métodos para calcular la integral de una función entre dos límites. El 
método de los trapecios es el mas sencillo de utilizar aunque también es el menos preciso. 
Consiste en calcular el área del trapecio que se forma con cada pareja de puntos respecto al eje de 
abcisas. La fórmula general, para un intervalo (ij), sería la siguiente: 


Área; = (x, —x;) . Sex) = (Base) - (Altura_media) 


Se suman las áreas de todos los subintervalos, abarcando el intervalo global definido por los límites 
de integración. 
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Problema 5.9 


Se desea obtener acetato de butilo en un reactor discontinuo de tanque agitado (RDTA) que debe 
trabajar a 100 “C (isotérmicamente). La mezcla inicial estará constituida por butanol y ácido acético, 
en la proporción molar de 4.97:1. Asimismo, contendrá un 0.032 % en peso de H2SOx, que actúa 
como catalizador. Calcular: 


(a) El tiempo necesario para alcanzar un grado de conversión de un 50 %. 

(b) Determinar, para el grado de conversión indicado en el apartado anterior, la cantidad 
inicial de cada reactante por carga y el volumen útil del reactor necesario, si se desea 
obtener una producción de 50 kg/h de éste, utilizando un solo reactor y sabiendo que la 
carga, descarga y puesta a punto entre cada dos operaciones consecutivas, cuesta 30 
minutos. Se partirá de la mezcla indicada en el primer párrafo. 


Datos y notas: 
Para el sistema en cuestión, la velocidad de reacción es de segundo orden respecto al ácido acético: 
3 
r=kc?  k=029.10* 
kmol- s 


Densidades a 100 “C (kg/ mé Zi Ácido acético: 958, Butanol: 742, Acetato de butilo: 796. 

Pesos moleculares (kg/kmol): Acido acético: 60, Butanol: 74, Acetato de butilo: 116. 

Aunque la densidad de la mezcla reactante varía un poco con la conversión, el exceso de butanol 
utilizado reduce la magnitud del cambio. Por tanto, en primera aproximación, la densidad media del 


conjunto puede considerarse constante a lo largo de la reacción e igual a 750 kg/ Mm. 


DATOS 
Se asignan referencias a los distintos compuestos: 


Ácido acético: A1 (clave), Butanol: A2, Acetato de butilo: A3, Agua: A4 


El esquema de reacción será: —A¡ A) +HA3+A4=0 
siendo: v==1. vy3=1 me=l o vy¿=1 
z . kg 
Densidades: Pmedia =750 
s 
Aya YE pia PE pm EE pyi=1000 E 
m? m? m? m? 
Pesos 
ias rua E a E Pic. pue E 
MdIecHIares: kmol kmol kmol kmol 
Ecuación de velocidad: r=k>c? Constante de velocidad:  k:=.29.10* Y — 
N kmol + s 
Mezcla inicial: 2 24.97 
N, 
Contenido en ac. sulfúrico despreciable (catalizador). 
Grado de conversión: X:=0.5 
Producción buscada: m3:=50 E 
A 
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SOLUCIÓN 


Para calcular el tiempo de reación se tiene que plantear el balance de materia al RDTA. 


: 0d 2 d 
El balance de materia, en función de la concentración, es: —<C¡=v*r=R; 
dí 
ed d S 
Para esta reacción: —C =v r=v¡*k*C; 
di € 
: : 1 de 1 1 1 1 1 
Separando variables e integrando: — dt= . - t= de,= . =- 
v¡*k C; v¡"k Ey v¡*k C10 Ci 


Cl0 


La concentración inicial (ci) puede calcularse a partir de la densidad y el peso molecular medio de 
la mezcla. Antes, deben calcularse las fracciones molares iniciales (despreciando el ácido sulfúrico). 


aye sai 9708 
144.97 144.97 
ki 
PM medio *=X10* PM¡+xX29*PM3= 71.7 El 
kmo. 
Pmolar *= ENSES. =10.5 kmol 
medio m 


NOTA: No es necesario utilizar las densidades individuales. 


En donde ya se puede calcular la concentración inicial: Ci0"=Pmolar*X10=1.8 A 
m 
la concentración final: c¡:=Cjp* (1—X)=0.9 pr 
m 
. em 1 lo 1 3 
y el tiempo de reacción: A =—|= (1.97 105) s Creacción =32.8 min 
v¡*k NC  C 


Para calcular el volumen, se debe conocer el caudal volumétrico de producción y el tiempo total de 
un ciclo. 


Como el tiempo de trabajo también incluye la carga, descarga, etc., del reactor, el tiempo completo 


del ciclo será: 
L 


completo *— Ereacción 


+30-min= 62.8 min 


La producción esperada es de 50 kg/h, por lo que en unidades molares: 


m 
net A 
PM; min 
Aplicando la relación estequiométrica a A3, su concentración será: 
v 
Czp:=0-+ kmol C3*=C3p—= Cp X+ z =0.9 kmol 
m v; m' 


y el caudal volumétrico puede calcularse a partir del caudal molar y la concentración: 


3 
A o a 


C3 Ss 
Finalmente, el volumen buscado es: V:=Qy" teompleto=514 L 
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Problema 5.10 


En un reactor discontinuo de tanque agitado (RDTA) cerrado de 1 L de capacidad, se hacen 
reaccionar 2 mol de hidrógeno con 1 mol de yodo, a 823 K, para obtener yoduro de hidrógeno de 
acuerdo con el esquema de reacción: 


Ha + 12 S 2HI 


La reacción es reversible y de segundo orden total en ambas direcciones. A la temperatura de 
reacción antes mencionada, las constantes de velocidad de las reacciones directa e inversa valen 


nr 


respectivamente 1.25" 10* y 0.25" 10 ==y 
kmol- s 


¿Qué tiempo se requerirá para obtener 1 mol de HI? 


DATOS 
Se realizan las siguientes asignaciones: Ai=12 (clave), A2=H2,  A3=HI 


De este modo, el equema de reacción queda: -A1 - A2 + 2A3 = 0 


siendo: vy:=-1 v):=—1 v3:=2 
Volumen reactor: V:=1+L=0.001 m' 
Cantidades iniciales: Nip:=1 + mol N»p:=2 + mol N3p:=0 + mol 
Concentraciones iniciales: apa iaa Lat Cro ¡220 7 tmol C30 ¡22% y Hmol 
V m? V m? V m? 
Cantidad y concentración final: N3:=1-mol C3 == 1 Era. 
m 


Teniendo en cuenta que se trata de una reacción elemental, reversible y de segundo orden en ambas 
direcciones, la ecuación de velocidad es: 


E) a > 2. 2 
r=k¿»C¡*C03=k;,"C;3 
3 3 
m m 


donde: ky:=1.25+10%.%— k,:=0.25+10%.%— 
kmol + s kmol + s 


SOLUCIÓN 
Para calcular el tiempo de reacción debemos aplicar un balance de materia al RDTA: 


d 
—c¡=v/"r=V¡" (hy» 0,03 ki"c3?) 
de 


en donde puede apreciarse que la velocidad de reacción depende de las concentraciones de todos 

las especies. Para poder integrar esta ecuación deberiamos saber la evolución de las 

concentraciones con el avance de reacción, lo que viene dado por su relación estequiométrica. 
c¡=C;p* (1—X) de, =—cC¡p* dX sustituyendo y simplificando: c,=1-—X 


sustituyendo y simplificando: —c,=2-X 


C3= Cp Cip* AX 


v] 
Y3 ] ; ds 
C3=C3p—Cjp* X= sustituyendo y simplificando: —c¿=2+X 
Y] 
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Sustituyendo estas relaciones estequiométricas en la ecuación del balance de materia: 
y 


Ar —V] (cu 0) (eu=co a TS Vv3 


dí Ci0 


v; 


1 


Antes de integrar es necesario determinar el límite superior de integración. A partir de la 


concentración final: 


v C3p—C 
C3¿=C3p—Cip* X+ o 1 Ale 
1 3 
E Cip* 
Y] 
Cc 
t:= 0 AX =(3.59-10?) s 
Ya Y3 
tatoo (10): [e=enrx ” TS ' ] ] 
1 1 


0 
t=0.996 hr 


COMENTARIOS 
La integral se puede resolver analíticamente, ya que el denominador es una ecuación de segundo 


grado y puede descomponerse en sus dos raíces (a y b): 
—In (x— a) + In (x—b) 


ln” le 


0 


X 


La integral también puede resolverse, de forma aproximada, por trapecios. Para ello, se define la 


función a integrar: 


F):= 


C10 


v v 2 
talon 0): [en=cix 2) eye > | 
1 1 


Se calcula el valor de la función para varios valores de X, preferiblemente en intervalos regulares: 


F(0.1)=(4.7.10%) s F(0.2)=(5.682.10*) s 


F(0)=(4-10%) s 
F(0.4)=(9.615-10*) s 


F(0.3)=(7.156-10*) s 
Se calcula el área de los trapecios formados en cada intervalo y se suman. Como los intervalos son 
regulares, su amplitud puede sacarse como un factor comun: 


F(0.5)=(1.455-10*) s 


a (edo ) a (pen o), al sH02) ) 3 ¡qen pelos) m ( HEl04) ) 


t=3.643-10* s 
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Problema 5.11 


Un vinatero aficionado quiere preparar vino espumoso en el sótano de su villa. Para iniciar su 
producción ha comprado unos botellones de 5 litros de capacidad, que rellenará con mosto, cuya 
concentración en azúcares es de 120 g por litro. Los botellones fermentarán en un soporte de 
invención propia, que los va girando lentamente a lo largo del día. La reacción de fermentación del 
azúcar se puede modelizar aproximadamente por una reacción irreversible de segundo orden, cuyo 


esquema de reacción es: 


3 


CoH1206 > 2 CH3-CH20H + 2C02 r=0.072 — 
kmol- h 


Como el vinatero desea obtener un vino espumoso con 50 g de alcohol por litro, determinar el 
tiempo de fermentación necesario para obtenerlo, suponiendo que los botellones se comportan 
como un RDTA. Desgraciadamente, el aficionado olvidó que el CO2 sólo se disuelve en el líquido en 
una proporción de 0.2 mol/L, y el resto pasa a la fase gas. Si en los botellones el volumen que 


3 ee / A e: A : 
queda un volumen de 100 cn7 entre el líquido y el tapón, y dichos recipientes sólo resisten una 
presión de 2.5 atmósferas sin estallar, ¿llegará el vinatero a catar el producto de sus desvelos? 


Datos y notas: 


Suponer que la temperatura del sótano es estable e igual a 10 “C. 


DATOS 
Hacemos las siguientes asignaciones: 


Con lo que el esquema de reacción es: 


siendo 


Los pesos moleculares son: 
El volumen de las botellas: 
La concentracion inicial de glucosa: 
La concentración inicial de alcohol: 


La concentración final de alcohol: 


Dado que se trata de una reacción elemental, irreversible y de segundo grado: 


SOLUCIÓN 


(a) El grado de avance para alcanzar la concentración de alcohol buscada es: 


C¿= Cay c10-X- 


Para este grado de conversión, las concentraciones de glucosa y CO2 serían: 


A =C4*(1-X)=0.123 


Ai: Glucosa,  A2: Etanol, A3: CO2 
-A1 + 2A2+ 2A3=0 
v¡:=—1 v):=2 v3:=2 
PM,=180:-2_  pm,:=46. 9 PM, :=44."9 
kmol kmol kmol 
V:=5+L=0.005 m* 
C10mási kmol 
Cl0másica *= 120 + dd Ci0 A GET ma 
1 m 
kmol 
Ca :=0. a 
a. 
gm Comási kmol 
Camásica *=50 + Ca a” 1.087 ex 
r=k»c,? k:=0.72 + 107?. == — 
kmol - hr 
14 Con —=C 
z y, E 
vi Va 
Cip* 
y 
kmol kmol 
a esi=0,=1.087 0 
m' m' 
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El tiempo de reacción puede calcularse mediante el balance de materia al RDTA: 


d 
—e =1 r=1,+k.c/? 
dt . 
z - de; 1 1 
separando variables e integrando: dt = ti= * | —-dc,=38 day 
v¡«k.c,? Vik |c? 
y Les : mol 
(b) La cantidad máxima de CO2 disuelto es: Camas :=0.2 + 
El volumen de gas es: Vg:=0.1-L 
Temperatura : T:=283:K 
Los moles de CO2 generados en la reacción serán: N3:=cC3+V=5.435 mol 
Los moles de CO2 que pueden disolverse son: Na3mar'= C3max + Y =1 mol 


Por lo que los moles de CO2 que pasaran al gas serán (generados - disueltos): 
N3s =N3¿— Na3mar =14-435 mol 


Considerando el volumen superior, Vg, y la ley de los gases perfectos, la presión que se alcanzaría 
sería: 


P:=Ny,+R- 7 =(1.03-10*) atm 


Por lo que no soportaría la presión 
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Problema 5.12 


En un reactor continuo de tanque agitado (RCTA) se planifica desarrollar en fase líquida la 
transformación de un componente A en otro B. La reacción que tiene lugar es irreversible y puede 
representarse por el esquema: 
A>B 
(a) ¿Qué tiempo espacial y qué volumen de reacción se requerirán para alcanzar una 
conversión del 75 % si el componente A se alimenta con un caudal de 1000 kmol/h y su 
concentración es de 1.2 kmol/ m7 ? 
(b) ¿Cuál será el tiempo espacial y el volumen de reacción necesarios si se duplica el caudal de 
alimentación, manteniendo constante la composición del mismo? 
(c) ¿Qué tiempo espacial y que volumen de reacción se necesitarán para que la concentración 
de A pase de 2.4 kmol/ mn” 20.3 kmol/ Mm , Si dicho componente se alimenta con un caudal 
de 1000 kmol_A/h? 


Datos y notas: 
La relación entre la velocidad de reacción, expresada como kmol de A reaccionados por má y 
minuto, y la concentración molar de A, se detalla en la tabla siguiente: 


CA 0.1 0.2 0.3 0.4 0.5 7 
(kmol= m7) 


0 
Qro Pos] 05 Pos Jo 


El esquema de reacción es: -A+B=0 siendo vw i=-1 vg:=1 


Tabla de datos anterior. 


SOLUCION 
El balance de materia para un RCTA es: Co CF ter=0 
Aplicado al reactivo A, que es el clave, resulta: Cova rer=0 
a 5 ll CA —= CA0 
Despejando el tiempo espacial: (¿a KA 
V¿r 
La concentración de A y el grado de conversión estan relacionadas: €4=Cyo" (1—X) 
kmol kmol 
(a) Caso 1 X,:=0.75 F 194 *= 1000 + Cs00=1.2 + 
hr m? 
. z kmol 
Las concentracion final de A es: €4:=Ca404* (1 —X,)=0.3 
m? 
o , Sr . kmol 
Para esa concentración, la velocidad de reacción (tabla) es: r4:=0,5 +. 22 
min - m* 
ó A (E E 
Luego el tiempo espacial es: A o A 
Va"Fa 
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El volumen del reactor se puede obtener a partir del caudal molar, la concentración y el 


tiempo espacial. 


F 3 
Oya :=%=0.231 42 V,:=0y.1=25 m* 
CA0a s 
(b) Caso 2 X,:=0.75 isa EA ji 
hr m? 
En este caso se duplica el caudal molar, pero se mantienen las concentraciones de entrada y 
salida y, por tanto, también el grado de conversión, X. 
El tiempo espacial es el mismo: 1=1.8 min 
EF 3 
Oy := 2% =0.463 Y V,:=0yp-1=50 m? 
CA0b s 
Si se ha duplicado el caudal, la única forma de mantener el tiempo espacial y, como 
consecuencia, el avance de la reacción, es duplicando el volumen. 
(c) Caso 3 F40.:= 1000 + a Cao :=2.4 + pr 
4 m 


En este caso se aumenta la concentración de entrada, pero la concentración de salida debe 


seguir siendo la misma. 


El nuevo grado de conversión será: 


y el nuevo tiempo espacial: 


Xx 4074 0875 


c 
CA0c 


Ca Cao, E 
7:= E =4.2 min 


G 
Vi*Fa 


El volumen del reactor lo calculamos como en el caso anterior: 


F ; 3 
Qy,.:=“%£=0.116 
C40c s 


V.:=0Qy,*1,=29.167 m? 


En este caso, el volumen del reactor ha aumentado ligeramente respecto al primer caso. 
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Problema 5.13 


En un RCTA ideal isotermo (funcionando en régimen estacionario) de 120 litros de volumen, tiene 
lugar la reacción (en fase líquida) cuyo esquema se describe a continuación: 
ki 
A+BSR+S 
ka 
Al reactor se alimentan dos corrientes distintas con el mismo caudal volumétrico. En una de ellas 


(exenta de B) la concentración de A es 2.8 kmol de A/ mM, y en la otra (exenta de A) la 


.r 3 A .r A 
concentración de B es de 1.6 kmol de B/ mM. Si se desea obtener una conversión del reactivo 
limitante del 75 %, ¿cuál será el caudal volumétrico de cada una de las dos corrientes alimentadas? 


Datos y notas: 
Admítase que la densidad molar de las dos corrientes alimento y de la corriente resultante tienen 
prácticamente el mismo valor. 
Experimentalmente se ha deducido que la velocidad de reacción puede expresarse mediante la 
ecuación: 

r= Ky» Cq* Cg— Ko* Cp* Cs 


nP 


kmol- mín i 


E 


m 
— a la temperatura a la que se desarrolla el proceso. 
kmol- mín 


siendo ki = 7 


Segun el enunciado, el compuesto B esta en menor proporción, por lo que será el reactivo 
limitante. 


El esquema de reacción es: -A-B+R+S=0 


siendo: — v¿:=—1 vgi=—1 Vp:=1 vg:=1 
La velocidad de reacción es: r=k,»C¿*Cg—=Kk>*Cp* Cs 
siendo: k¡:=7 mm k,:=3 m 
kmol « min kmol « min 
Volumen del reactor: V:=120 L 
Avance de la reacción: X:=0.75 


El alimento está compuesto por dos corrientes que se mezclan en la entrada del reactor. Ambas 
tienen el mismo caudal pero diferentes concentraciones. El caudal será la suma de caudales, 
mientras que las concentraciones se calculan según se muestra. 


Corriente 1: Oy C41:=2.8 dls 
m? 
. . kmol 
Corriente 2: O0v7=0yw  Cg:=1.6 
m? 
y Cir? Csjo C C 
Concentración de A en entrada: — c¿y= “a On Ca Ori - Car je =la di 
Oy +0vw  2*Qy, 2 2 m 
Concentración de B en entrada: — czy= ca Qv _ cur 0r - a1 diia gran 
Ov+0v  2*Qy 2 2 m 
Las concentraciones de productos en el alimento son nulas: — czy:=0 EA Esp = un 
m m 
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SOLUCIÓN 


El balance de materia al RCTA es: Cro A+ vz r=0 


Para expresar la velocidad de reacción en función del grado de conversión se tienen que obtener las 


relaciones estequiométricas: 


v 
CA = Cao Ego X* 4 800 Mol 
Vz m? 
A E pa 
mm 
v 
epa 0 PA 
B m? 
Y 
Cs =Csp— Cgp* X- y A 
V m? 


Sustituyendo en la ecuación de velocidad: 


Va Vr Vs mol 
ri=k;* Cao Cgo*X-* + (cg (1—X)) —k+ Cro—Cgp*X-* + Cso—Cgp*X-* =0.667 5 
Vg Va V ms 
Ñ E —Cgp* X 
puede calcularse el tiempo espacial: T:=—=—=900 s 
VgrY 
V L 
y el caudal de la corriente alimento: Oy:===8 — 
T min 
Por lo que le volumen de cada corriente será la mitad: Oy =0y,=4- = 
min 
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Problema 5.14 


Se desea llevar a cabo la reacción: 2A1 5 A2 + A3 en uno o más reactores continuos de tanque 
agitado conectados en serie. El alimento se introduce con un caudal de 3.5 m / h. El alimento está 
exento de Az y A3, conteniendo Ai con una concentración de 48 kmol / nm: 


(a) Si las concentraciones finales de A2 y A3 han de ser iguales al 85 % de las de equilibrio, 
Calcular el tamaño del reactor que se requiere. 
(b) Si los recipientes de que se dispone tienen un volumen igual al 20 % de la capacidad del 
reactor único, calcular el número de RCTA en serie que se requerirán para alcanzar el 


grado de conversión deseado. 


Datos y notas: 


La constante de velocidad de la reacción directa, a la temperatura de reacción (que debe ser igual 
para todos los reactores), vale 2"10-4 ná | kmol"s y la constante de equilibrio es igual a 16. 


El esquema de reacción es: 
El caudal de entrada: 
Las concentraciones de entrada: 


Datos cinéticos y de equilibrio: 


-2A1+A2+A3=0 


siendo: v¡:=-2  v,:=1 v3:=1 
m? m? 
O,0:=3.5 E = (9.722.107) E_ 
hr s 
kmol kmol kmol 
C1p:=48 A Cap = a C3p:=0 a 
m m m 
3 
kys=2+ 107% Y K,:=16 
kmol » s 
C=0.85+C7,,  C3=0.85+C3, X=0.85+X, 


SOLUCIÓN 


El primer paso es calcular la conversión y las concentraciones en el equilibrio, para ello, hacemos 


uso de las relaciones estequiométricas: 


=Cjp* (1—X,) 
=C2p—C1p*X¿* a e ci0*%e 
V 2 
=C39—C1p*X,* dE = C10*%e 
Y 2 
Cio*X, C10*X, 
Kg _ Cog*Czp _ 2 2 w Cro X _ XxX. 
Ki Ci (cjo* (1—X,))* 4=cp «(1-X,)?  4=(1-X,)” 


con lo que se tiene una expresión que relaciona la constante de equilibrio con la conversión en el 
equilibrio. Dado que conocemos Ke, resolvemos para Xe, ya que se trata de una ecuación de 


segundo orden. 


EAS solve, X, 
—_— 


8 
A 7 
4. (1-X,)? 3 
9 


De las dos soluciones matemáticas, solo tiene sentido la segunda, ya que el grado de conversión 


debe ser menor que la unidad: 
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Con lo que también pueden conocerse las concentraciones en el equilibrio: 


-X, 
Cje3=Cj0* (1 —X) =5.333 nal aj o FA esa A 
m? 2 m? m? 
Dada la condición impuesta: X:=0.85+X,=0.756 


las concentraciones finales serán: 


4 
ai O e ia 
) 3 


m m m 


e 


Aplicando el balance de materia al componente clave, se podrá obtener el tiempo espacial y el 
volumen del reactor, para el caso de un reactor único: 


Balance de materia: Co X + vpo T-r=0 
.r s 2 A k 8 m? 
Ecuación de velocidad: r=kysc? —li«039 03 siendo : k,:=2%= (1.25.107*) 
K, mol + s 
Sustituyendo, despejando y operando: Cy X+vj+ te (gc? —Kk;c7+c3)=0 
O (75 s 1=0.215 hr 


y," (loc) —k,*C7*C3 


T=— V,:=1+ Op =0.753 m? 


En el caso de instalar una secuencia de RCTA iguales, el volumen y el tiempo espacial de éstos será: 


3 Ve; 
V,,:=0.2+ V,=0.151 m 1,:=—=0.043 hr 
vO 
El balance de materia, para cada reactor de la serie es: 
Cro" (A, — Xi) + Vo Tp 7. =0 siendo n las condiciones del reactor y m la conversión 


a la entrada, proveniente del reactor anterior 
sustituyendo y reordenando: 


Co (XX) +7 +59 (hoc —k;"c3c3)=0 


. . 2 
Cip* (XX) + Vo Tj* l (co (1 — Xi) )' —k;* (222) =0 


ey Epa lt (cjo* (1 —X,))? —Kk;* pu 


—V¡*T; 


2 C10* A ] 
Ar A = S ka" (ejo= (1 —X)) —k;* > 


C10 


V]*T; 2 Cj0* A, Ñ 
Ay = A, + ? Ka" (Cro (1 —X,)) —k;* 
Cro 2 


expresion que nos permite calcular la conversión a la entrada del reactor, si conocemos la 
conversión a la salida. 
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Ahora el cálculo del número de ractores puede hacerse de varias formas: 


Opción 1 (Desde el final hacia el principio), hasta que la conversión de entrada sea inferior a cero: 


Di*T; Cy + X. 2 

n=1 Xi =X=0.756 — X,=X,+ “tanteo x,))? al a y as 
Cio 2 
DVy"T, (co XV? 

n=2 X,=0.604 Xz:=X2+ Claro (19) 0 [2] Joox00 
Ci0 2 
Vy*T; Cia XaV* 

n=3 X,=0.156 Xy =X¿+ Cltarlens(1 xy" ao [ME )+-1.059 
Cio 2 


Hacen falta tres reactores. 


Opción 2 (Desde el principio hacia el final), hasta que la conversion de salida sea superior a la 
buscada: 


-T, > 
Xmm1:=0 po Eoltatons( Xn))' ta [Su y Xin =0 


Vi*T; 3 Cro" Xi Y 
n=nl1 X,,:=0.1 Xp :=r00t |X,1 + . ka" (Cro (1 X»1)) k;* Amr XAn1 [50.56 


7, e. 
v o ltastono(1= 29) 44 [200] Jm] 073 


n=n2 Xu =Xn  Xn¿:=root . + 


«T, «X.Y? 
Miri lt. (cro* (1 Xm9))? ki» — ] me) =0.812 


n=n3 Xms=X23  Xn“:=root . + 


Matemáticamente es mucho más sencilla la opción 1, ya que no hay que calcular la raíz en cada 
reactor, pero no nos da la conversión de cada reactor. 
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Problema 5.15 


En un reactor cuyo comportamiento fluidodinámico puede asimilarse al de un RCTA se desarrolla la 
reacción gaseosa reversible de primer orden: 
A1 S A2 
Operando a 300 K el volumen de reactor necesario para obtener una conversión del 60% es de 
0.1 mM. 

(a) Utilizando el mismo caudal volumétrico de alimento, calcular el volumen de reactor 
necesario para obtener el mismo grado de conversión (X=0.6), si la temperatura de 
operación es de 400 K. 

(b) ¿Cuál será la temperatura de operación que haría mínimo el volumen del reactor, para el 
caudal de alimentación y el grado de conversión anteriores? 


Datos y notas 
A 300 K la constante de equilibrio referida a las concentraciones y la entalpía de reacción toman los 


valores respectivos de K = 10 y AHok = -33490 J/mol. 
La constante de velocidad de la reacción directa varía con la temperatura según la ecuación: 


k=10% «exp (22 57 


El calor específico molar de las dos especies químicas que participan en la reacción puede admitirse 
que tiene el mismo valor para cualquier temperatura. 


DATOS 
El esquema de reacción es: -A1+A2 =0 
siendo : v¡:=-1 v,:=1 


Comportamiento del reactor a una temperatura de 300K: 
T3p0 += 300 K X300*= 0.6 V, “== 0.1 mi? 


Datos termodinámicos a una temperatura de 300K: 


- mo 
Datos cinéticos: —2400 Ej 
k,=A,y" ex, =Á¿" ex 
d=Ag »| T ) d ol 
siendo:  4¿:=1000+s”! E,s=2400+K+R=(1.995-10%) E 
mol 
—E k 
Ka 300 =4,o| ed Jos A k; 300=2% = 0.034 1 
- R- T300 s Ñ Ko 300 s 


Temperatura de operación del reactor: — T¿yy:=400 K 


SOLUCIÓN 


El primer paso para resolver el problema es deteminar el caudal volumétrico de alimento, que se 
extrae del balance de materia panteado a 300 K 


. V, 
Balance de materia: Co X+vporoTt=0 siendo r=— r= ky»c¡—k,*C) 
Oro 
A E ha v 
Relaciones estequiométricas: c¡=Cjp* (1—X) A O CE IÓ 
v 
Sustituyendo las relaciones estequiométricas en el balance de materia y operando: 
v, 
Cp X+v]* (% ¿+ c,—Kk,C7) «—=0 
vo 
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C10*X + v+ (Kyo (e70+ (1 —X)) —K; + (070 X)) * y =( 
VIA 
Cro" | X+ ln 40)] =e 
vO 
v, 3 
Oo += ar (Ka 500" (1 —X500) — ki 300* X300) =0.019 4 
—A300 s 


(a) Para calcular el volumen de reactor necesario a 400K, primero hay que determinar la 
constante de equilibrio a dicha temperatura 


—4H 
K. =K'.*exp|—— 
R-T 


Ed : K 
El término preexponencial K' se calcula con los datos a 300K: e O ÓN 1.476-10* 
exp 
R 7) 
A , —4H 
y la constante de equilibrio a 400K será: Ko 40 =K"* exp =0,349 
"400 


Recalculando las constantes cinéticas a 400K: 


—E 
haa, | pl Ja. Keaoo j 
*L400 c_400 


Ka 40 =7110 1 


Se puede replantear el balance de materia a 400K, podrá obtener el volumen de reacción en las 
nuevas condiciones 


a «X. 
V,:= 0" X00 =-3.482 L No tiene sentido un volumen negativo 


j y" (a 400 (1 —X300) —E_400* X300) 


Podemos comprobar, que la conversión máxima a 400K es la siguiente: 


K Cde _  Cip*Xze ui Ko 400 
q ATEN e 
1+K, 400 


= 0.258 
Cle  Cio* (1 —X) 


(b) Para calcular la temperatura óptima, hay que poner todas las variables en función de la 
temperatura y calcular el máximo de la velocidad de reacción 


—E¿ 
_E A" exp 2D 
¡=400p|7 2) (eus(1=99) - er Cp X 
K".+ exp a ) 
—E ÁA¿*Cip*X —E¿+4H, 
reordenando r= Ay» cio" (1—X) - exp e a a (E) 
A . 
d E —E, AgrCipo"X (E¿— 4H —E¿+4H 
—r=AguCip* (1=X): 2. a * exp 1 _ A : e . » exp E A =(0 
dT R T R=T E R RAT 
operando: 
di Toi. = 3290 
R 1 (1 —X00) “Eg 
B X 
BA (E,—4H) 
Ko 
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Problema 5.16 
Se quiere desarrollar en condiciones isotérmicas, la reacción en fase gaseosa: A> 3 B 


En la tabla adjunta se detalla la forma en que varía la velocidad de desaparición (por reacción 
química) del componente A, en función del grado de conversión, a la temperatura de 150 CC a la 
que se desea que opere el reactor en el que tendrá lugar la transformación, a la presión de 
10 bar abs. La mezcla inicial consistirá en una mezcla equimolecular de A e inertes, a la temperatura 
y presión de trabajo del reactor. 


ESA! [0 E A A E LA E 


a 53 52 so [as |a30 3 [25 lis [is 
(kmol/m? s) 


(a) ¿Cuál será la forma más conveniente de representar los datos de la tabla con vistas a 
calcular el volumen de reactor necesario para alcanzar una conversión dada, en un RCTA, y 
en un RFP? 


0.85 


10 


(b) ¿Qué información adicional sería necesaria en su caso, para calcular el tiempo de reacción 
requerido para alcanzar una conversión del 70 % en un RDTA? 


(c) Si el alimento se introduce al reactor con un caudal volumétrico de 2070 L/h. ¿Cuál será la 
concentración molar de A a la entrada (cao)? ¿cuánto valdrá el caudal molar de A (Fao)? 
Supóngase comportamiento ideal. 


(d) Cuál es el volumen de reactor necesario para alcanzar una conversión del 80 % en: 
- un RCTA 
- un RFP 


(e) Considérese un sistema formado por dos RCTA conectados en serie. Si a la salida del primer 
reactor se alcanza una conversión de A del 40 % ¿cuál será el volumen total de reacción 
(suma de los volúmenes de los dos reactores) para que la conversión final de A sea del 
80%? 


(f) Si se desea alcanzar una conversión de A del 80 %, ¿cuántos m3 de volumen de reactor se 
ahorran si se utilizan dos RCTA montados en serie, con una conversión intermedia del 
40 %, en lugar de un solo RCTA, para obtener una conversión global del 80 %? 


(g) Repetir el apartado anterior pero comparando el volumen de reacción necesario cuando se 
utilizan dos RFP en serie, y un solo RFP. 


(h) Caso de que se pueda disponer de un RCTA y de un RFP, con los volúmenes que se desee, 
y se monten en serie, de modo que la conversión intermedia sea de 0.7, ¿cuál deberá 
montarse en primer lugar, el RCTA o el RFP, para un volumen global de reacción más 
pequeño? Demostrar la decisión con cálculos. 


(i) Repetir el apartado anterior para una conversión intermedia de 0.5. ¿Qué reactor se situaría 
primero? 


El esquema de reacción es: -A+3B=0 siendo vyi=—1 vg:=3 


Dado que solo tenemos un reactivo (A), éste será el componente clave. 


Temperatura y presión del reactor: T:=423 K P:=10 bar= (1 10% ) Pa 
J 
Constante de los gases perfectos: R:=8.313 —=— 
mol» K 
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Fracción molar del alimento: 


Velocidades de reacción 
utilizadas en los cálculos: 


X4p:=0.5 
sto EA 
d 3 

ms 
at A 
úl 3 

m Ss 


stos E 
: 3 
m .s 
jotas PL 
Ñ 3 
m Ss 


SOLUCIÓN 


(a) El volumen de un reactor puede calcularse a partir del balance de materia en un RCTA y en un 
RFP, respectivamente: 


1 


—vV4 Y 


«X 


RCTA: V=F¡p" 


X 
RFP: rar | 


0 


En ambos casos es conveniente representar la inversa de la velocidad de reacción, referida al 
componente clave, frente al grado de conversión. 


dx 


—Vy*Y 


-Ra 1/-Ra 1000 
X | (kmol/m?-s) | (m3-s/kmol) 900 
0 0.0053 189 800 
0.1 | 0.0052 192 o 700 
0.2 0.005 200 2 600 
0.3 | 0.0045 222 > 500 
0.4 0.004 250 A 40 
0.5 | 0.0033 303 E 
0.6 | 0.0025 400 o en 
0.7 | 0.0018 556 = 200 
0.8 | 0.00125 800 100 
0.85| 0.001 1000 0 


0 0,11 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 
X 
(b) Tiempo de reacción, en el caso de un RDTA para X = 0.7. (No considerando que este tipo de 
reactores no se suelen utilizar para reacciones en fase gas.) 


Al tratarse de una reacción en fase gas, con variación en el número de moles (supuestos P y 
T constantes), debe tenerse en cuenta la variación del volumen al progresar la reacción y 
cambiar el número de moles. 


V= Vo" (1+8,-X) 


_Nio dx Neo dX Co dx 


dt = II — — — _ _—_——————————————A E > 0x->-——* 
Vo —v pr  Vg- (1 +e,X) —Vp Y (1 +8, X) — VÍ e Y 
Xx 
1 
t= 8 => dx 
— Y r+ (1 + € X) 
0 
En el problema: A 2 o =xyr0=1 
n el problema: Vi=V¿+ VBg= = = E =3Xgp*0= 
—V y 
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X -Ra 1/Ra(1+X) | Trapecios La concentración inicial se puede calcular: 
(kmol/m*-s) | (m*-s/kmol) | (m*-s/kmol) _P=xg 0.142 kmol 
0 0.0053 189 AO 3 
m 
0.1 0.0052 175 18.2 
0.2 0.005 167 17.1 Calculamos la integral, |, por trapecios, desde X = 0 
0.3 0.0045 171 16.9 hasta X = 0.7: 
0.4 0.004 179 17.5 a 
0.5 0.0033 202 19.0 1:=140.1 
0.6 0.0025 250 22.6 mol 
0.7 0.0018 327 28.8 
0.8 0.00125 444 t:=Cyp"1=19.9 s 
0.85 0.001 541 
| SUMA 140.1 | 


(c) La concentración de A a la entrada ya se ha caculado. A partir de ella y el caudal 
volumétrico se puede calcular el caudal molar: 


3 
Oy0:=2070 > (5.75-107%) — Oca) = 01080 mr 


(d) Si se utiliza un RCTA o un RFP desde X = O hasta X:=0.8 


da 
En un RCTA: PA EA 


—TVaA*Fos 


3 
mes 


X 
1 
En un RFP: Vomp =F yo" | ———— dX 1,:=261.7 
RFP A0 E ) d nal 
0 


—v¿*Y 


1000 ? 


1/Ra (m3*s/kmol) 
u 
o 
o 


0 01 0/2 0,33 0,4 0,5 0,6 0,7 0,8 0,9 
X 


El volumen necesario en un RCTA es muy superior al de un RFP. 
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(e) Si se utilizan dos RCTA en serie, con una conversión intermedia de X = 0.4: 


Foo X 
RCTA1- X,,:=0.4 Viena SL 
0.4 
Fao* (X7. —X0) 
RCTA2 Xo :=0.8 Vrecraz == =26.2 L 
0.8 
1000 , 
> 800 4 
le) 
É 
7 600 
E 
m 0 VrozaL'= Vgcrar + Vrcra2=34.3 L 
E 
se RCTAZ En consecuencia, el volumen requerido 
á pato para llevar a cabo la reaccion es menor 


cuando se utilizan dos RCTA en serie. 
0 01 012 0,33 0,4 0,5 0/6 0/7 0,8 0,9 


ys 


(f) El voumen ahorrado al utilizar los dos RCTA en lugar de un solo RCTA es: 


V, Y. 
Ahorro :=E14__" TOTAL — 34 49% 


Vec TA 


(g) Si se utilizan dos RFP en serie, con una conversión intermedia de X = 0.4, los volumenes serían: 


Xlg 
3 
RFP1: X¿:=04 — Vip=Fao [ar 1834.25 PO 
(—v4+7) kmol 
0 
Koa 
3 
RFP1: 0 Pas A 1,1784 HA 
(—v4-1) kmol 


Xig 


Vrorag*=Verer1 + Vrer2=21.4 L 


X -Ra 1/-Ra  [Trapecios| SUMA 1000 7 
(kmol/m'-s) (m*»s/kmol) RFPi , 

O | 0.0053 189 _ $00 

0.1| 0.0052 192 19.0 3 

0.2 | 0.005 200 19.6 E 600 

0.3 | 0.0045 222 21.1 ES 

0.4 | 0.004 250 23.6 83.4 E 400 

0.5 | 0.0033 303 27.7 3 

0.6 | 0.0025 400 | 352 + 200 

0.7 | 0.0018 556 47.8 

0.8 | 0.00125 300 67.8 178.4 0 

0.85 0.001 1000 0 0,11 0,72 0,3 0/4 0,5 0,6 0,7 0,8 0,9 

SUMA | 261.7 261.7 X 


El volumen total es el mismo, indistintamente de si se usa un RFP o dos reactores RFP en serie. 
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(h) Si se dispone de un reactor de cada tipo y la conversión alcanzada a la salida del primero ha 
de ser de 0.7, vamos a comparar las dos configuraciones posibles para determinar la mas 
favorable: 


(h1) RCTA + RFP 


Fs X, 
RCTA-H1: Xi ¡:=0.7 Vicia =—*2318 E 
_ - 10.7 
Xu > 
1 m+s 
RFP-H1: Xu ,:=0.8 Verp =Fap* =——— dx Ly :=67.8 Verp m3 =Fao* Ly =35.5 L 
> (—y e r) kmol Ñ 
Xu 
Ñ VroraL m1 =Y reza m1 + Vree 137.3 L 
(h2) RFP + RCTA 
Xul 
1 m+s 
RFP-H2: Xy 1:=0.7 Verp=Fo* dX Ip := 194 Verp m =Fao* lp =15.9 L 
a —v¿* r) kmo Ñ 
0 
For (Ximo 2 AX, 
RCTA-H2:  X;,, ,:=0.8 nds ¡2 For (na 2 Xm0.1) 65 E 
o B 70.8 
VroTaL_H2 E Verr m2 ap Vecra_m2 =224 L 
1000 ? 1000 , 
= 800 5 800 
E E 
+ 600 600 
a E 
€ mM 
0 E 400 
[2 RCTA o 
E ez < 
200 Y 200 E 
RFP EZ 
0 0 l 
0 01 012 03 0/4 0,5 0,6 0,7 0,8 0,9 0 01 0,12 0,3 0,44 0,55 0,6 0,7 0,58 0,9 
X X 
La mejor opción es la H2 (RFP - RCTA). 
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(MD) Si la conversión intermedia es de 0.5: 
(11) RCTA + RFP 


Fo" X, 


RCTA-l1: Xp /:=0.5 Vacra y 2 ALL 12,4 L 
- Po.5 
Xu 2 
3 
RFP-11: A O E A A A 
Ñ (—v e r) kmol a 
Xo1 
Vrora_n == Vecran ap Verrn =24.7 L 
(12) RFP + RCTA Xo 
1 3% 
REP-12: Xp /:=05 Vesp=Faoo A —dX lg ld 2 Vamp pi=Fsorlp=9.1 L 
e (—v o r) kmol Ñ 
0 
For (Xp »—X, 
RCTA-12: Xy 3:08 ros p:=0.00125 et a Ca > Xn 1) =19.6 L 
mos P0.8 12 
VroraL_n = Verp 1 ap Vecra_n =28.7 L 
1000 ? 1000 9 
800 a 800 
3 3 
S 600 É 600 
7 EN 
E 400 E 400 
Eo 2 
=> 200 Y 200 
RCTA RFP RCTA 
0 0 
0 01 02 0,3 0,14 0,55 0,6 0/7 0,8 0,9 0 01 0,12 0,3 0,4 0,55 0,6 0,7 0,8 0,9 
X X 
La mejor opción es la 12 (RCTA - RFP). 
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Problema 5.17 
En un RFP se está descomponiendo acetaldehído según la reacción: 
CH3CHO > CHa + CO 
El reactor funciona a temperatura constante e igual a 791 K y a la presión de 1 bar abs. 


El volumen del reactor y el caudal volumétrico de alimento (acetaldehído puro), son tales que resulta 
un tiempo espacial de 1/8 h (referido a 273 K y 1 bar abs). 


.s . 2 
La ecuación de velocidad es de la forma: — r=X-c, 


siendo: 
c1 = concentración molar de acetaldehído (kmol/ 77 ). 
3 
k=0.33 LP. 
kmol- s 


Determinar el grado de conversión alcanzado en el reactor. 


DATOS 
Haciendo las siguientes asignaciones: Ai = Acetaldheido, A2 = metano, A3 =Monóxido de carbono 
El esquema de reacción es el siguiente: -A1+A2+A3=0 
siendo vy:=-1 v):=1 v,:=1 
El alimento es acetaldehido puro: Xx1p5=1 
Tiempo espacial: == hr=450 s referido a las condiciones: 7T,:=791K P:=1 bar 
m? 

La constante cinética es: k:=0.33 

kmol + s 
SOLUCIÓN 


Como el tiempo espacial esta referido a una temperatura diferente a la del reactor, hay que 
recalcularla teniendo en cuenta la ley de los gases perfectos: 


Q,0(273) Q,o(791) 791 


== ——— Ovo (791) 0 Ovo (273) - 


273 791 273 
PO PRAT a 1 L=15531 5 
dl 0027) —= y 
273 273 273 


Ahora ya se puede plantear la ecuación del balance de materia al reactor, teniendo en cuenta que es 
una reacción en fase gas con cambio en el número de moles. 


NOTA: Po y To son las condiciones inciales en un RDTA o las condiciones del alimento en un RCTA. 
En este caso, T=To y P=Po: 


y; 
Cro" 0d 


Ve Cs Av 


c¡= siendo: 0 == ó= Ex = Xyp * O 

Po T Ch0 — Ve 

RPP (1 +ej" Xx) 

P T Ci0 v; + 14) + v; 

0 0|=Ú=1 ó:=_"_£_ “1 E ¡=3X¡p"9= 
C10 —=V1 
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Teniendo en cuenta la ley de los gases perfectos, la concentración inicial puede calcularse: 


P 0.015 kmol 
“Ñ di 


Cp 


Con lo que la concentración de acetaldehido variará con el grado de conversión según la expresión: 
_ Cp” (1 —X) 
pe PEEL A AS 
1+-X 


Planteando la ecuación del balance de materia en un RFP, resulta: 


X XxX XxX 
1 1 1 
T= Cro" | ———4AX=C 10" | ——————dX=Cyp* dX 
vr (Xx) —y ¡+ koc;? A 10" (1 —X) 2 
Vj,eko 
, 0 ! 1+X 
XxX 
2 
as 1 e 
—v¡*k*Cip (1 — Xx 
0 
que puede resolverse analíticamente mediante la expresión: 
UH ar prin 1H 
(1-X)* X-1 
A O A (PA E 
—v¡*k=Cip X-—1 0—1 


A X+4+ In pS a 
—v¡*k*Cip —1 X-—1 


que se cumplirá cuando: 
1 


—v,*k-Cip 


x;=root| «(x 4+ 7 A 


NOTA: Cuando la incógnita es un límite de integración, la integral no puede calcularse. Si no se 
pudiera resolver analíticamente, puede hacerse un tanteo, por trapecios, con límites crecientes hasta 
que se cumpla la igualdad. 


Sea la función: F(X):= 


2, F(0 prElO) ¿0,1]=130:459 5 
E «ko. Ce 


a AMA +F(0.1) ) 0174001) +F(02) y =93.153 s 
pa vj+keCjo *C10 a 


Ñ F(0 plo) -0.1+ 


e ek» eN 


-0.1+ 


F(0.1)+F(0.2) F(0.2)+F(0.3) 
e) 3 


. 01] =36.364 s 


ón Fl0)+F(01) 01 dE F(0.1) +F(0.2) dd F(0.2) +F(0.3) e F(0.3) +F (0.4) o 
v ko» AN 2 2 2 


por lo que el grado de conversión estará entre 0.3 y 0.4. 
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Problema 5.18 


En una fábrica de ácido nítrico el gas que sale de la planta de oxidación de amoníaco se enfría 
rápidamente, hasta la temperatura ambiente, para eliminar la mayor parte del vapor de agua que 
contiene. La mezcla reactante, una vez enfriada, contiene un 9 % de NO, 1 % de NO» y 8 % de O»2 
(porcentajes molares). Antes de introducirla en las torres de absorción donde se obtiene ácido 
nítrico, conviene oxidar la mezcla a fin de que la relación NO2/NO en la misma sea 5:1. 

Si dicha oxidación se realiza en un reactor tubular, con una refrigeración adecuada para que la 
temperatura de la mezcla se mantenga constante e igual a 293 K, ¿cuál será el volumen de reactor 
necesario, cuando se introduzca un caudal de alimentación de 10700 m7 /h (medido a 293 K y 
P = 1 atm), si la presión de operación es 1 atm? 


Datos y notas: 
Se sabe que la reacción 2:NO + O2 5 2:'NO2 es una reacción gaseosa homogénea prácticamente 
irreversible, que se comporta como una reacción elemental. 


La constante de velocidad a 293 K vale 1.4: 10* referida al NO, cuando las concentraciones se 
expresan en kmol/ m y el tiempo en segundos, es decir: 


—R yo) = 1.4. 101 Ep” *Co> 


Despréciese la posibilidad de que se pueda formar algo de N204, 


DATOS 

Haciendo las asignaciones: A1=NO A2=02, A3=NO2 

El esquema de reacción es: -2A1-A2+2A3=0 

siendo v¡:=-2 v):=—1 v3:=2 


Las presión y temperatura de operación del reactor y del alimento son las mismas: 


P:=1 atm T:=293 K 
A 4 mé 
La constante de velocidad es: k:=1.4.10% —=— 
kmol? +5 
3 
El caudal de alimentación es: Oyp:=10700 3 
> 
y la composición del alimento: x1p3=0.09 Xp 5= 0.08 X3p:=0.01 
donde el componente clave de la reacción es el NO: 210 =0.045 2210.08 
vi v> 
z . N 
Razón molar final: == 
N, 


SOLUCIÓN 


Se trata de una reacción en fase gas con variación en el número de moles y las condiciones (P,T) 
del reactor son las mismas que las del alimento: 


peas 
C;50 == C; A+ 
Ve : , Av 
e = ——— siendo: /v:=v,+v,+v3=-—1 d:= =-0.5 Ep=0+x¡p=—0.045 
1 + Ex . X Y] 
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Para relacionar las concentraciones con las fracciones molares, hacemos uso de la ley de los gases 
perfectos: 


N, P, 
qa Ls Pp -x; PP 41599 mol 
V RT RT R-T m? 
l P ti" (1-X 
C0:= -«210=3.743 22 e, (X):= Ed at, 
R-T mi RT 1+€,-X 
v) 
1 P có UD 
mo 
Cap i= == 387 — ca (X):= . . 
R-T m* R-*T 1+€,:X 
va 
1 P Lay Tip X+ a 
is «Lgy=0.416 22 Cy (X) :2 + 
R-T m* R*T 1+€,:X 
Para hallar la conversión final hay que partir de la condición impuesta a la composición de salida: 
v3 
Lgy Lip A+ 
N C v x 5 
A FRA - x= 2 0.815 
N; Ci Lp (1—X) vz 
Lip 5 * 10 
XxX 4 x 
, Vv 1 1 
Planteando el balance de materia: LE dX= dx 
Cro Fo Up er — Aix 
0 
X 
mol 1 
Fi=0n "6 =11:126 V=Fio* | ——=— dx 
Ss k+c¡% *Ca 
0 
X 
1 3 
VE y IX =1107.7 mM 
lo 
Loy Tip A+ 
P L10* (1-X) e P vz 
R-T 1+€,:X R TS Ep, XxX 


0 


Integral resoluble numéricamente (Simpson o trapecios): 


F(X):= - 


b P Lio + (1—X) A P UN 
R*T 1+€,:X 


La solución por trapecios sería: 


3 3 3 
F(0)=1.532 2 F(.2)=2.625 F(4)=52 2 
mol mol mol 
3 3 
F(0.6)=13.313 27 F(.815)=73.889 22 
mo mo 
V,:= F(0)+F(0.2) 7, de F(0.2) +F (0.4) 0d F(0.4) +F (0.6) ba F(0.6) +F(0.815) 7715 Fo 
2 Z 2 2 
V,=138.2 m? 
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BLOQUE IV 


Balances de energía 
en sistemas con reacción 
quimica. Reactores Ideales 


Problema 6.1 
La reacción de deshidrogenación del etanol para obtener acetaldehído: 
C2H50H(g) > CH3CHO(g) + H2(g) 


se lleva a cabo en un reactor con calefacción para evitar que la temperatura disminuya demasiado 
rápido y por tanto se pare la reacción cuando la conversión todavía es reducida. 

Si la composición del alimento es de un 90% de etanol y un 10% de acetaldehído y su temperatura 
300 “C, la temperatura de salida de los gases del reactor es 253 %C y en la misma hay 41 moles de 
acetaldehído/100 moles de alimento. 

Calcular el calor adicionado por cada 100 moles de gas alimentado. 


Datos y notas: 

Cp(etanol) = 0.110 kJ/mol'K, Cp(acetaldehido) = 0.080 kJ/mol'K, Cp(H2) = 0.029 kJ/mol*K 
AHr(etanol) = -56.12 kcalmol, AHr(acetaldehído) = -39.76 kcal/mol 

Suponer que las variables termodinámicas son independientes de la temperatura. 


DATOS 
Asignando : A1: Etanol A2: Acetaldehido A3: Hidrógeno 
Esquema de reacción: —-A1+A2+A3=0 donde v¡:=-1 wi=1 Y i=1 


Las entalpías de formación y los calores específicos (independientes de la temperatura) son: 


AH y :=—56120 a O €p/:=110 — 
mol mol mol +K 
A A A cp7i=80 — 
mo mo mol» K 
J J 
AH 3:=0 C,3:=29 
di mol de mol» K 
La temperatura del alimento es: T,:=300 *C 
La temperatura del reactor es: T:=253 *C 
La composición del alimento: x]p:=0.9 X2p:=0.1 X3p:=0 
El componente clave es el etanol, por ser el único reactivo: 
Suponiendo una base de cálculo: Fi):=100 se 
Ss 
mol mol mol 


El caudal del alimento es:  F,,:=F,9+x¡9=90 


F3y:=F9*x2p=10 F3y:=F9*x39=0 —= 
Ss Ss 


41.  moles_A2 _ mol 


100  moles_alimento s 


En la salida del reactor: salida = 


SOLUCIÓN 


Se trata de un reactor continuo, donde se lleva a cabo una reacción endotérmica con aporte de calor. 


TT, =HX=0 
—AH, F, AH, 
A C 0 
siendo: J=_ 0=——L— y = 44,2% 
Sa) PE) Ero Y 
Jj “pj k0 Jj pj 
j=1 j=l 
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Para plantear el balance de energía es necesario conocer la entalpía de reacción y la conversión 


alcanzada. 

La entalpía de reacción se calcula a partir de las entalpías de formación a la misma temperatura: 
4H,= Y) (9: 4Hy) AH, +=," AHy +v3+ AH +v3 + 4H = (6.85 +10*) — 

j mo 


La conversión se calcula a partir de la información de la corriente de salida, utilizando las relaciones 
estequiométricas: 


FF 
Fy=Fyp=F ip" xo despejando: — X:i==20"2=0.344 
Y] va 
Fip* 
v] 
Calculando J: 
az 
v; mol 
Ss 
Fi Fx J 
C= ), (0,* cy) Ci=L0 cy ++ Cp = 118.889 
j=1 Fo 10 mol» K 


donde únicamente se tienen en cuenta los componentes del alimento: 


_—4H, 
C 


J: =-576.135 4*C 


resolviendo el balance: 
O:=T-—T,—J.X=151.446 4*C 


y teniendo en cuenta la definición de Q, el aporte de calor será: 


q:= (Fo C)-0=(1.62-10%) w 
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Problema 6.2 
La reacción irreversible de primer orden en fase líquida: 
A>R+S 


se lleva a cabo en un RCTA hasta alcanzar una conversión del 90 %. El reactante se alimenta al 
reactor a 20 %C y los productos se extraen a 150 %C, 
Determinar: 


(a) El volumen del RCTA necesario para tratar 1000 kg/h del componente A. 
(b) La energía que ha de intercambiar el reactor con el exterior. 


Datos y notas: 

La constante de velocidad a 130 “C, 170 “C y 190 “C vale 0.067, 1.713 y 7.010 (1/h), 
respectivamente. 

Calor de reacción a 25%C = -52250 julio/mol. 

Densidad de reactantes y productos 900 kg/ mM. 

Peso molecular de A = 100 g/mol. 

Calores específicos medios: Cpa = 50.6 J/mol'K, Cpr = 45.4 J/mol*K, Cps = 60.5 J/mol'K. 


DATOS Y CÁLCULOS PREVIOS 


Esquema de reacción: —-A+R+S=0 donde v¿:=—1 Vp:=1 vg:= 1 
La temperatura de entrada al reactor es: Tp:=20 *C=293.15 K 
La temperatura de salida del reactor es: T:=150 *C=423.15 K 
La conversión en el reactor es: X:=0.9 
El caudal másico de entrada al reactor es: m¿:=1000 + Ñ =0.278 kg 
r Ss 
3 
Teniendo en cuenta la densidad: p:=900 a O :=".= (3.086 . 107*) Ei 
m 12) Ss 
: kg my mol 
Teniendo en cuenta el peso PM,:=100+ Fjo:= =2.778 
molecular: kmol 1 Ss 
e J J 
Calores especificos, constantes  c,,:=50.6+ Cpr:=45.4 + Cps:=60.5 + 
. l-K mol» K mol + K 
con la temperatura: mo 
La entalpía de reacción a 25%C: AH,,y:==52250 » > Tyg=25 "C=298.15 K 
: mo 
La entalpía de reacción a la temperatura del reactor será: 
y 
AH, = AH, ¿+ | Y (Yy* Cp) AT =4H yg + 2 (Y Cp3) * (TT 3) 
Ñ J 
Tre 
AH, :=AH ey + (va *Cpa VR" Cor + VS" Cps) . (T— Tres) =-4.534-10* — 
mo 
Constante cinética: 
E 0.067 1.861.107] 403 
¡=4:0p| 2 k:=|1.713|.—=|4.758.10*| = T,:=|443 |.K 
R+T 7.010] % |1947.107| * 463 
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Para poder calcularla a cualquier temperatura es necesario hacer el ajuste del modelo: 


a —10.892 ,  [o.o0248] , 
In(k)=In (4) -— In(k=s)=| —7.65 —=|0.00226 | — 
R=T 6241 DT loo0216| E 


de la pendiente se obtiene la energía de activación: 
WE 


pendiente := slope E In[k- ») =-1.446-10* K E¿:==—R+ pendiente = (1 202 +10? ) a 
mo 


r 


de la ordenada en el origen se obtiene el término preexponencial: 


ordenada := intercept la ln (9) =24.996 Ayi=2P (ordenada) _ (7173.10%) E 
r Ss K 
y la constante cinética a la temperatura de reacción es: 
—E 
k(T) =Ag* exp | —L K(T)=1.028-10* L 
R-T s 
SOLUCIÓN 
. Fo. X 
Planteando el balance de materia: V, === 
— VÍ" Y 
siendo: Fio= Q,*Cx0 r=k:C, cy=Cyp* (1 —X) 
"CN X 
sustituimos : ya (Orta) Xx 
—Vp* (k (cxo* ( —X))) 
el volumen del reactor será: 
.X 
(a) y,= 90 =27.018 m' 


Y planteando el balance de energía: 
g= Y (Esc (177) +" AH, 2) 
j 


La energía necesaria por unidad de tiempo será: 


4 
(b) q:=F¡p*Cpa* (TT) + Fo" 4H, + X=-=9.507-10% W 
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Problema 6.3 


Una mezcla que contiene dos reactantes A y B, en igual concentración (6.67 kmol/ nm ), se alimenta 
a un RCTA, adiabático, donde reacciona de acuerdo con el esquema siguiente a T = 120 %C, 


A+B>C 
3 


La reacción es irreversible y de segundo orden, y su constante de velocidad, expresada en mal 
mol. s 


puede calcularse por la expresión: 


k=5.5-10 a 7) 


¿Cuál deberá ser el caudal volumétrico y la temperatura de entrada del alimento, si el reactor tiene 


un volumen de 1 litro y se desea alcanzar un grado de conversión del 50 %. 


Datos y notas: 
Calor de reacción AH = -80000 kJ/kmol (a 120%C). 


Calor específico medio de la mezcla 2720 kJ/ ne K (a las condiciones de operación). 


DATOS 
Esquema de reacción: -Ai -A2 +43 = 0 siendo v¡:=—1 v,:=—1 vs:=1 
Calor de reacción a 120%C: AH, :=—80000 e 
mol 
Calor específico referido al volumen: Cpy=2720 + ed 
m? -K 
7 r nm? 
Parámetros de la constante cinética: Agi=5.5+107 «—L— 
kmol + s 
E,=10000 K+R=(8.314+10*) LE 
mol 
Concentraciones inciales de reactivos:  c,,:=6.67-+ al Cp :=C¡p=6.67 pal 
m m 
Volumen del reactor: V,:=1+L 
Conversión en el reactor: X:=0.5 
Temperatura del reactor: T,:=120 *C 
SOLUCIÓN 


Para calcular el caudal volumétrico tenemos que recurrir al balance de materia: 


v, 
Ci X+v¿eror=0 Ci XA+v¡+k=C,*C7* *=0 
vO 
= —E, 7 m? 
siendo : k(T):=4¿+exp k(T,) = (4.934 + 107”) 
R=T s- mol 
C¡=C3=Cy* (1 —X) 
, Y, 
sustituyendo: Cc X+v,+k(T) + (c10-(1-X))? «E =0 
vO 
k(T 1-X))' el 
v,* ( ) "(er l > )) a 
—lP..o 
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LU IATA O 


3 


O,o:=v1 k(T,) * (co (12M)? A 0 
AC 10* ol 


v, 
Tí= =607.701 s 
vO 


Para determinar la temperatura de entrada del alimento, se ha de recurrir al balance de energía de 
un RCTA adiabático: 
O=T-=T, IX O=0 T=T.=J-X 


Para Calcular la J] no podemos usar directamente los datos del problema, ya que son 


dimensionalmente inconsistentes: 
—AH, AH 
Ja an Gt E 
—v; mol 


La entalpía de reacción está referida a los moles, mientras que el calor específico está referido al 
volumen. Se puede solucionar de la siguiente forma. La energía necesaria para calentar el alimento 
es la misma, independiente de su cálculo: 


Ss 
Fjo* Y (0y* cy) = Oo * Co las unidades son: kJ/(K-s) 
j=l 
Ss 
operando: (0,0 y C10) á yn (0, Ñ Cri) = Do d Cy 


5 
c 
_— pv 
2 (0, 1) = 
j=1 Co 


—AH, 
de =196.176 K 


C10 


Ahora, ya es posible calcular la temperatura de entrada del alimento: 


Tp:=T,—J+.X=295.1 K 


Índice 
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Problema 6.4 


(a) Calcular entre O y 130 %C la conversión de equilibrio de Ai para la reacción en fase acuosa 
A1 S Aa2 para la cual AG298 = -3375 cal/mol y AH298 = -18000 cal/mol, dando el resultado 
en forma de una gráfica donde se represente la conversión frente a la temperatura. 

(b) Qué restricciones se tendrán que imponer a un reactor que opere isotérmicamente, si se 
desea obtener un grado de conversión del 75 % o superior. 

(c) Si el reactante Ai se introduce continuamente a 20 “C en un reactor que opera 
adiabáticamente, el cual es suficientemente grande para considerar que se alcanza la 
conversión de equilibrio, determinar: 


- La ecuación de la adiabática de conversión. 
- La temperatura de salida del reactor. 
- La composición de la corriente de salida. 


Datos y notas: 
Suponer que los calores específicos de las disoluciones son idénticos a los del agua. Suponer que la 
disolución es ideal, por lo cual K = Kc. 


DATOS 
Esquema de reacción: -A1 + A2 =0 siendo: v¡:=—1 v,:=1 
Datos termodinámicos: T,¿:=298 K 
AG,.f:==3375 A EE Fr 
mol mol 
cal a y 
AH,,y:==18000 L=-7.536+10* — 
mol mol 
C,:=75 Es Constante e igual a la del agua 
mol + K 
Temperatura mínima: Tu =0 PE 
Temperatura máxima: Fi=130*6 
SOLUCIÓN 


(a) Cálculo de la curva de equilibrio: 


La constante de equilibrio depende de la energía temperatura y la variación de energía libre: 
—A 
K.= exp pta 
R-T 
Teniendo en cuenta que: 4G=4H-—TAS, también puede espresarse mediante la variación de 
entalpía: JH 
Ko 0 Ko * EXP | —— 
R*T 
Si conocemos la variación de energía libre a una temperatura, podemos calcular la constante de 
equilibrio a dicha temperatura: 


=AG jej 
Ko ref:= exp | | = 299.869 
o Trej 
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Conociendo la constante de equilibrio a una temperatura y la variación de entalpía, podemos calcular 
el factor preexponencial y utilizarlo para calcular la constante a otras temperaturas: 


Kos —AH ef 
Ko i= 2 1.848 +10"! K¿(T) := Kg exp | 2 
—AH jej | R=+T 
exp | ———= 
R ld To) 


A partir de la constante de equilibrio y su expresión en función de las concentraciones, podemos 
obtener su relación con la conversión de equilibrio: 


c,=C jp" (1—X) c¿=Cp"X 
k=2= Cio*X, = Xo 
cy Cipr(1-X) 1-X, 
despejando Xe: K.(7) 
X,(1):= 
1+K.(7) 

X.(273 K)=1 X,(293 K)=0.998 X.(313 K) =0.986 
X.(333 K)=0.925 X.(353 K)=0.724 X,(373 K)=0.398 
X,(393 K)=0.161 X,(403 K) =0.098 


—9— Curva de equilibrio 


0.6 
x< 05 ———RCTA isotermo 
0.4 RCTA adiabático 
0.3 
0.2 
0.1 
0 
273 293 313 333 353 373 393 413 
T (K) 
(b) Reactor isotermo: 
“Y no Xo, 
Para alcanzar una conversión máxima de: Xo A > Kc:=—==3 
1-Xo, 
AH 
Teniendo en cuenta: K. = Ko » exp 
R*T 
La temperatura máxima a la que se podrá trabajar es de: 
A A 100 
isotermo “— RN > isotermo — 
C 
R«In|— 
Kop 
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(c) Reactor adiabático: 


Temperatura de alimentación al reactor: Tp:=20 *C 
Balance de energía al RCTA, adiabático: T=T,+3*X 
Ñ —AH os 4H 
cálculo de la J: AH,:= E 75410? Ji= L=1005 4%C 
v; mol 6, 


La conversión máxima alcanzable está en el corte de la recta adiabática con la curva de equilibrio: 


Recta adiabática: Toa =T, +3 og 
Aa í AH. 

Curva de equilibrio: 9d =K.y"exp El 

l—X ad RTad 

xx —AH o 

Combinandolas: e =K.y" exp dl 

1 —X ad R+(T,¿+J+ Xq) 

X, t K, Arg |x| =0.107 
ad*=Y00 —Kop* Xp | ——————|. =U. 
dl E ARANA 
Xo4=0.107 TE AU a 
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Problema 6.5 
En una instalación de producción de hidrógeno se obtienen 10 Tm/día. El proceso se basa en la reacción. 
CO(g) + H20 (g) >C02(g) + H2(g) 


que se realiza haciendo fluir una corriente de CO y vapor de agua por un lecho catalítico, a presión próxima a la 
atmosférica y en régimen estacionario. El reactor está refrigerado mediante una camisa por la que circula agua. 
La mezcla alimento se introduce en el reactor a 150 “C y el producto se extrae a 480 “C. Para conseguir que la 
conversión del CO sea total, se emplea un 50% de exceso de vapor de agua sobre el estequiométrico. 

Calcular los kg/h de agua de refrigeración que serán necesarios si su temperatura no ha de aumentar más de 
20 “C respecto de su valor de entrada. 


Datos y notas: 
Capacidades caloríficas molares medias entre 25 y 1500C (en J/mol*Kk): 
Cp(CO) = 29.30 Cp(H20) = 34.16 (gas). 
Capacidades caloríficas molares medias entre 25 y 480%C (en J/mol*K): 
Cp(CO) = 44.79 Cp(H20) = 35.58 (gas) Cp(H2) = 29.30. 
Entalpías de formación a 298 K y 1 atm (en J/mol): 
AHf(CO2) = -393702 AHf(CO) = -110577 AHf(H20) = -241942 (gas). 
DATOS 
Establecemos las siguientes asignaciones: A1=AK=CO A2=H20 A3=CO2 A4=H2 
Esquema de reacción: -A1-A2+A3+A4=0 siendo: — v,:=-—1 vy:=—1 v:=1 vy=1 
Producción másica de hidrógeno: 00. 115.741 EL 
24:-3600 s s 

Producción molar de hidrógeno: PM,:=2+£% Py e a A 

mol PM, s 
Exceso de agua: 50% 
Grado de conversión: Xp:=1 
Temperatura de entrada: T,:=150 *C 
Temperatura de salida: T,:=480 *C 
Temperatura de referencia: Trefp:=25 “C 


Entalpías de formación: 


AH y :=—110577 + =d AHp := —241942- de AHg := 393702 + AH pa :=0+ 
] mol : mol mol mol 
Calores específicos medios: 
. J J 
Corriente entrada: c,,,:=29.3+ C,.p:=34.16+ 
pe mol» K e mol» K 
eiii E cy3,1=44.79 Coq 29.3. 7 
Corriente salida: i mol» K 7 mol» K E mol+«K 
Datos de la refrigeración (agua): 
Cambio temperatura: AT,:=20+4%C 
Calor específico agua: Cpai=75 * -b 
mol» K 
Peso molecular agua: PM, :=18+. € 
mol 
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SOLUCIÓN 


Planteando el balance de energía: A(F-h)=q 


Los caudales molares se obtendrán a partir de las relaciones estequiométricas: 


F 
PO fs F¡:=F jp" (1-X)=0 Ec 
Xy a S 
Fap:=Fyp* (1+0.5)=86.806 "2 E O CI 
Ss == V] Ss 
y O A 
== V] Ss 
E F,=F,=57.87 D2L 


Ss 


La entalpía molar específica de cada especie hay que obtenerla a partir de la ecuación de Kirchoff: 


T, 
Ah =4hy+ | e, dT=4h,+c): (T-T,) 


To 


Las entalpías de las especies presentes en la corriente de entrada son: 


hyp3= AH y +Cy19* (Tp —T) =-1.069+ 10% E 
: Í mol 
h30=4Hp + Cpa0* (Ty Tre) ==2.377 + 107 o 
i 5 mol 
Las entalpías de las especies presentes en la corriente de salida son: 
has :=4H pg +Cp95* (T, Te) =-2.258+10% LE 
: mol 
h3,:= 4H g + Cp3s * (T, — Tref) ==3.733 + 10* zA 
Ñ mol 
hio=4 Hg Cas (T,—T,) =(1.333+10*) 
j h mol 


Ya se puede calcular la cantidad de calor transmitido por unidad de tiempo: 


qu= (F)+ has + Fo hz, HF ¿oh 5) — (Ejo* hijo + Fzo* 20) ==5.466 + 107 W 


Este calor se ha retirado por el agua de refrigeración: 
q=Fy Ci AT, 
luego el caudal de agua necesario será: 


q 1 


Fy=%L_=364.4 9 


Cra* AT, s 


EM oa 
Ss 
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Problema 6.6 
La reacción endotérmica en fase líquida 
A+B>2C 


tiene lugar, de modo casi completo, en un RCTA en régimen estacionario. Determinar la temperatura de 
operación suponiendo que la agitación dentro del reactor es perfecta y teniendo en cuenta los siguientes datos: 


- El volumen del reactor es 475 litros. 

- El área de intercambio de calor es 1 m”. 

- La temperatura de condensación del vapor de calefacción es 185.5 %C, 

- El coeficiente global de transmisión de calor es 850 W/(/m? :k). 

- La potencia del agitador es 18.5 kW. 

- El calor de reacción es 46520 kJ/kmol (independiente de la temperatura). 
- Alimento 10 kmol/h de cada una de las especies A y B. 

- Temperatura del alimento: 25 “%C, 

- Calores específicos en kJ/(kmol*K): Cpa = 102, Cps = 102, Cpc = 110. 


DATOS 
Esquema de reacción: -A-B+2C=0 vyi=—1 vg:=—1 Vc:=2 
Entalpía de reacción (independiente de temperatura): 4H, :=46520 e 
m 
Calores específicos: 24102 ES Cpg += 102 AS Cpc:=110 e 
Temperatura entrada al reactor : Tp:=25 *C 
Volumen del reactor: V,:=0.475.-m?* =475 L 
Área del reactor: Aj:=1m* 
Coeficiente global intercambio calor: U:=850 lid 
m*»K 
Temperatura vapor del intercambiador de calor: — 7,,,:=185.5 *C 
Potencia aportada por el agitador: w:=18500 W 
Caudales: is A Fso:=Fro 
hr s 
Conversión: X:=1 


SOLUCIÓN 1 
Planteando el balance de energía en su forma extendida: 
Y (Ejo* Cp * (TT) +Fro* 4H, X,¿=q+w 


J 


El calor intercambiado se puede expresar: q=U-4,- (Trap — T) 


Sustituyendo y despejando: 
Fy9* Cp (TT) +Fog* Cp" (TT) +Fip-4H,*X=U+ A, (Tygp —T) + w 


2+F ig" Cp (TT) +F,0: 4H, X= U- Ay + (Tap — T) + w 
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w—F py 4H, 
Mi= 4 


SOLUCIÓN 2 


.X+ A AS par To 
2+Fip*cp4+U+A, 


=48 11 E 


También puede resolverse por la ecuación condensada, incluyendo wen el parámetro Q: 


L= To =J«X= Ocoregida 


Se calcula la J: 


—4H, 
AH, := 
Va 


= (4.652 + 10*) 


Se calcula la Q corregida: 


—AH 
e Ji= = 228.0 49C 
mol Fi + Fx 
. A el B 
Fi j 10 ] 
UA, (Trap = 1) +w o 
Oeoregida*= 2 1) =246.2 4%C 


Fi0* Cpa + Fog" Cog 


y, despejando en el balance de energía, se obtiene la temperatura del reactor: 


e O meca +J-X+ T, =43.1 *C 
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ANEXOS 


Formularios 


Índice 


Nomendlatura 


A lo largo del texto, debido al software 


estandarizado, como: 


utilizado algún símbolo de unidades difiere del 


Nombre de las unidades Unidades SI Unidades 
Horas h hr 
Días día day 
Toneladas t tonne 
Bloque 1 
bolo esla Nombre de la variable Unidades SI 
variable 
Aj Cualquier especie química que intervenga en el proceso 
Ak Componente clave 
E Entrada. Cantidad de materia que cruza la superficie kg/s o kmol/s 
límite del sistema de fuera a dentro, en un tiempo dado. 
S Salida. Cantidad de materia que cruza la superficie límite kg/s o kmol/s 
del sistema de dentro a fuera, en un tiempo dado. 
Acumulación. Diferencia entre la cantidad propiedad 
A existente dentro del sistema en un tiempo t + At y la 
cantidad de propiedad existente en un tiempo t. 
Generación. Cantidad de propiedad que aparece en el 
GeE sistema en un At sin estar presente inicialmente ni kg/(m3*s) 
haberse transferido a través de los límites del sistema. 
mi Caudal másico total de la corriente ¡ kg totales/s 
sa Caudal másico del componente A; en la corriente ¡ kg A; en la 
mn corriente ¡/s 
e Fracción másica del componente Aj en la corriente ¡ kg Aj/kg 
a totales 
F; Caudal molar total de la corriente ¡ mol totales/s 
Fs Caudal molar del componente Aj en la corriente i mol A; en la 
E corriente ¡/s 
mm Fracción molar del componente Aj en la corriente ¡ mol Aj/kg 
a totales 
PM; Peso molecular de la especie A; kg Aj/kmo! Aj 
Densidad: masa contenida en un volumen de control en kg/m3 
P un tiempo t 
pm Densidad molar mol/m? 
Paj Densidad de la especie A; kg/m? 
PMAj Densidad molar de la especie A; mol A; /m? 
Mo Masa total contenida en el sistema en el instante inicial kg 
M Masa total contenida en el sistema en un instante dado kg 
Vo Volumen de control inicial mi 
NA Volumen de control en un instante dado mi 
V máx Volumen máximo mi 
Fracción másica inicial kg Aj 
107) iniciales/kg 
totales 
t Tiempo Ss 
QLi Caudal volumétrico total en la corriente ¡ m/s 
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Bloque II 


Símbolo de la 


. Nombre de la variable Unidades SI 
variable 
Cp Calor específico J/(kg * K) 
Qy Caudal volumétrico m?/s 
z Altura M 
T Temperatura K 
q Flujo de calor W 
Qesp Flujo de calor específico W/kg 
Wominal Potencia nominal W 
Wieal Potencia real W 
Weosp Potencia específica W/kg 
n Rendimiento 
v Velocidad lineal media de la corriente m/s 
a Coeficiente adimensional dependiente del tipo de flujo Adimensional 
h Entalpía específica a una temperatura dada J/kg 
Epot Energía potencial J/kg 
Ecin Energía cinética J/kg 
Ctotal Energía total J/kg 
P Presión Pa 
. 8.314 
R Constante de los gases ideales J/(mol-K) 
S Superficie m? 
p Viscosidad dinámica kg/(m's) 
Re Número de Reynolds Adimensional 
D Diámetro de la conducción m 
U Coeficiente global de transmisión de calor W/(m?-k) 
C Coeficiente de descarga del caudalímetro Adimensional 
Y Factor de expansión Adimensional 
Bloque III 
Simbolo de la Nombre de la variable Unidades SI 
variable 
Uj Coeficiente estequiométrico de la especie Aj 
Mo Masa de la especie A; en el momento inicial kg Aj 
M; Masa de la especie A; en el momento t kg A; 
Njo Moles de la especie Aj en el momento inicial mol Aj 
N; Masa de la especie A; en el momento t mol A; 
x Grado de conversión Adimensional 
Cjo Concentración de la especie A; inicialmente en el reactor mol A¡/m? 
C Concentración de la especie Aj en el reactor mol A¡/m? 
Ck Concentración del componente clave kg Ax/m* 
R; Velocidad de reacción intensiva respecto a Aj mol A; /(m3*s) 
r Velocidad de la reacción única mol/(m3-s) 
k Constante de velocidad 
Ka, ki Constante de velocidad directa e indirecta 
AO Factores preexponenciales de las constantes de 
d, Ai . 
velocidad 
E Ei Energía de activación 
Ko Constante termodinámica 
e Tiempo espacial de un reactor continuo de tanque E 
agitado (RCTA) 
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Bloque IV 


Símbolo de la Nombre de la variable Unidades SI 
variable 


Cociente entre la variación de la entalpía de reacción y la 
J capacidad calorífica del alimento por mol de componente K 
clave, cambiando de signo. 
Cociente entre el calor de intercambiado por el sistema y 


Q la capacidad calorífica del alimento. Ñ 
o, Razón molar de cada componente respecto al 
componente clave 
AHk Calor de reacción por mol de componente clave J/mol 
Xx Grado de conversión 
To Temperatura a la entrada del reactor K 
T Temperatura a la salida del reactor K 
AG Variación de energía de Gibbs J/mol 
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BLOQUE 1. Balances de materia en sistemas sin reacción química 


Ecuaciones generales: 


Balance de materia total 
(unidades másicas): 


Balance de materia total 
(unidades molares): 


Sin reacción química o sin 
variación en número de moles: 


Balance de materia a un componente 
(unidades másicas): 
Á 


Balance de materia a un componente 
(unidades molares): 


Sistemas finitos: 
BM total (másico): Jm > 
t 
sl dV 
BM total (volumétrico), p=cte : 44 
t 
AS dM,; d(M- 
BM componente (másico): — 4m += = | 6,47 4 (m-0w;) + aio) = | 6 dy 
Pay” de) 
dN, 
BM total (molar): AF + z d = | cav 
EE 
dN, d(N,- 
BM componente (molar): ar —=L= [6 av A (F,-x)) + y ») = | 6, av 
t t 
V 14 


Sistemas finitos, sín generación (G=0), en régimen estacionario (A=0): 


AA AS 
xq 

E 
=_ A 


Un 
a |3, 
a 


SS 
E q 
SI 1 
a, ES 
LA 


PA 
ta 
UA 


A 
> 

3 

os SS 
SS 
E E 


BM total (másico): Am=0 Sm= Y m (2) 
Salida Entrada Ss 
BM total (volumétrico): A0y=0 >)0/= > 0, mi 
Salida Entrada A 
Ss 
E kg_4; 
BM componente (másico): 4Am,=0 >m= Y m, 
Salida Entrada Ss 
kg A; 
a(moj=0 Elmo) E moy EL 
Salida Entrada Ss 
BM total (molar): AF,=0 YN F= Y Es => 
Salida Entrada Ss 
kmol_A; 
BM componente (molar): AF,=0 »NF= Y) Fo a) 
Salida Entrada Ss 
kmol_A; 
A 
Salida Entrada Ss 
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BLOQUE II. Balances de energía total en sistemas sin reacción química 


Balance de energía total 


te — — d J 
Ecuación general: [nse,+er) plo n)a4+— | (u+ep+er) pav=q+w+ [Gp av e 
dt Ss 

A Vv Vv 


En sistemas sin generación (GE=0) 


[n+e,+e;) plo n) d4+— f (use, +e) pdaV=q+w 7 
t 


A Vv 


si además, las corrientes están definidas y el sistema es estacionario: 
J 
=q+w — 


si además, solo hay dos corrientes (m = Mmentrada = Msalida): 


2 
me-v 


A(m-n)+9:2(m)4h-a| 


1 y? : q w J 
Ah+g- Az+ E A = Qesp + Wesp siendo: — qesp= e Wesp = pa 5) 
si además, no hay cambios de fase: h7=h%+c,*(T—T*) 
1 v? J 
a En 
Balance de energía mecánica 
E HUA d J 
Ecuación general: (P-v+e,+ez) plo 1) 44+— | (ep+ex) pav=w+ | Gp av =, 
A di Vv Vv ds 
En sistemas con corrientes finitas, sin generación (Ge=0), en estado estacionario: 
1 ú : E J 
A(P-v)+9-Az+=-A4 == Wi Es siendo: e, =— — 
2 Q m kg 


L : E 
para tramos rectos: — E,=4+m>-f E siendo f el factor de fricción 
para accidentes: E,=K.v* 


si además, el flujo es ideal (ev=0) y no hay bombas ni compresores (w=0): 


2 
A(pu0) eg: a2h-a[ E) =o 


para fluidos incompresibles (Ec. Bernouilli): 


AP 1 2 P,—P 1 (vw.  v,? J 
19) 2 2 
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Caudalímetros 


Caudal másico: 


z : 4 2-p¡- AP kg 
Diafragmas, boquillas y venturímetros: meC0+Y Ayu 11 anda 
1-8 Ss 
Rotámetros kg 
Ss 
Caudal volumétrico: 
3 
Diafragmas, boquillas y venturímetros: Qy=C-Y-A)- s a E 
pr: (1 =8 ) Si 
Factor de expansión: 
Fluido incompresible: Y=1 
E P, Co 
Fluido compresible: r=— y=—= 
P; Cy 
A , P; —P, 4 
Diafragmas: Y =1-——2.(0.41+0.35-8*%)  0.2<8<0.8 
VE 
S % P; Py 4 
Boquillas y venturímetros: Y =1 - ———>» (0.8+0.8-8 ) 
y-Py 
Rotámetros: Y=1 
Coeficiente descarga: 
Diafragmas: Figura 4.1 C=0.6 Re> 10000 
Boquillas y venturímetros: Figura 4.2 C=0.95 Re> 3000 
Rotámetros: Figura 4.1 
Pérdida de presión permanente: 
Diafragmas: A _ (1 -p?) j AP broma (Pa) 
1-6? 
Boquillas: AP rmatente = A MP pino (Pa) 
1+8 
Venturímetros 2 A AA (Pa) 
Rotámetros: O (Pa) 
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Figura 4.1. Diafragmas y rotámetros 


1,00 
Er En e? Rotámetro con V 
0.90 al flotador tipo e 
0,80 IMA ia 
200 Wi 7%" VISIO 
0,70 
0,60 
0,50 
0,40 
0,30 
10 20 4060 102 10 10% 105 
V¿D 
e A 


Figura 4.2. Boquillas y venturímetros 
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- 


- 


o o ooooo 
u) K 00 O 00D 
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BLOQUE III. Balances de materia en sistemas con reacción química. 


Estequiometría 


Número de reacciones: 


Esquema de reacción: 


Conversión molar: 


Conversión intensiva 


Reactivo limitante: 


Exceso reactivo Aj: 


Razones iniciales: 


Grado de conversión, X: 


Sistema estático 


Sistema dinámico 


Sist. estático/dinámico 
con p=cte ó Av=0 
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=1 R>1 
> ("4 =0 > (14; =0 i=1..R 
j=1 
N;=Njo a 
7 N¡=Np+ vj" N/=Np+ 2 (v ¿) j=1..5 (mol) 
hi e 
C;—Cj0 E mol 
pa C= Cjo RV" X E= 0 + lv; Xx) j=l..S ; 
v; 1=1 m 
; F. C; X: 
| min Dl min | min | min 
Y y; Y y 
N.,—N, MESA FEL—F, ACA east e EE cede OL 
50 k0 0 k0 Cjo— Cxo Xjo — Xk0 
. Ve Ve Ve Vi 
A A A A 
v v, v v. 
Nao — Fig — Cx0 — Xk0 — 
Ve Ve Ve Ve 
9= No Eo 2 o 7% 
'j 
Ni Fi Cro  *Xw 
Ya Nio=Ne _ Fio=Fr_ Cock _ Yo =Xk Y X= 
Nuo Fi Cx0 Xk0 i 
R 
Ni=Ngo* (1 —X) N¿=N0* [ ES 2] (mol) 
i=l 
y; S Vi 
N=Ny=NypX-— N=N0=Mjr DA (mol) 
VE i=l Vik 
y, £ Vi 
N=Nio* [0 X-—= N=Nw* 0, 314, (mol) 
VE i=l Vik 
ds ol 
Fy¿=Fg* (1—X) Fi=Po: [1 Ex) E 
¡=1 
v E Vi mol 
ESF y-Fy 2 Fy=Fy—Fir 4H ES 
Ve ¡=1 Vik s 
v; S Vi mol 
F,= Fig" |0,-X- Fy=Fig*|0,- E |%: 
Ve i=1 Vik Ss 
Z mol 
cp= Ci * (1 —X) ca (1- Dx) 
i=1 m? 
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R 
y mol 
SS , amore Er 2) 7 
Ve i=l Vik m 
Cin10.—X+ di 
N F. ea v Av 
Sist. gaseosos con ss LS E siendo: £¿=xyy9 ¿J=ÍT- 
Av, PÓT variables Yo 0 [Po T] (1+2,-X) má 
k 
P+ To 
Po, To: presión y temperatura iniciales o del alimento 
v. 
9x2 
A Lyn O N; VE 
Sist.estático/dinámico  x,=—=xyp"|————Ñ 
N, PS 
ea ¡ el MN ll E-=Fs | 
Rendimiento: o =Moles_ de Ay formados (Vel Uy" tjo Ye] _Cytgo (Ye 
moles_de_Ak_iniciales | v; No v; Evo v; 
Selectividad: Es moles_de Aj formados Vel _ N¡—Njo VI _ F; Fi Vir 
E  E>_ ———————— — —| == —__— o. |l|l|= 
moles_de_Ak_reaccionados |vy| NN |vyl Fiw=Fx | vi 
Cinética Química 
Velocidad de reacción intensiva: y a PRA rICOneS 
: sm 
% R; mol 
Velocidad de una reación: r=Ñ 
0 ms 
dz ] SS j mol_A; 
Relación entre r y Ry: R¡=v,*1 R¡= D) (v;-r;) j=elo.s —— 
i=1 ms 
Ss Ss 
Infuencia de la composición r=kp" Tc —k,* Ile” 
ca — y. j ¡=1 j=1 
Pb E e e j Ñ , A Ñ , a 
ya APraduicras e=to- |op=c0: Lx) tr Tlso—0w: Lx) , ] 
j=1 ve j=l ve ms 
k Ss Ss 
Equilibrio (r=0) k.=2= [Tot 4= ITc? 
Tr ¡Al Jal 
—Ep 
(Ep — Ej) 4H" 
; Kkp _ Áp*e .” Áp RT RT 
Influencia de la temperatura: K === =_e K =K-e 
K; ==. Y 
Ae RT 
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BLOQUE III y IV. Reactores ideales 


Balances de materia 
N=1 


BM RCTA: Cjo CF Ye Tr ZE A + Voto r=0 


N>1 


jm 


siendo m=n-1 


C; — Cp E V¡* To * Ln = Cp? (Xi — Xi) EV Tn *Fn 


QSL FLOG : Va 0" Yn_ 40 _ cp: 4X 
Oo Fr YT Vir “06 Fro Y¡*Fy VE" Tn 
Cho siendo 4c;=C; —Cim » AX=X, — Xi 
Y V, Y, 
BMRFP:O 1220” a de; y — a Te , de; 
O Fi v¡*r ¡UN Fio v¡*r 
XxX A 
C C 
t= L- € T, = Y 
—V Y —YpeF 
0 Am 
€; XxX 
de, 0 
BM RDTA: Eye Pe 7 o ay 
dt dt Cro ver —Vp o Y 
t=1tp lp=t +HIR+l+l, Vi=0y*tr 
Balances de energía 
N=i N>1 
BERCTA: T-T,-J-X=0 Ta— Ti — + (XX) = O, 
siendo: O= L O, = de 
Fo" Y (0, Cpj) Fw Y (0, Cy) 
j=1 j=1 
ya 4H, 
TS 
2 (8, Cp) 
j=1 
F, AH, 
q=E AH¿==L 
" Fg VE 
BE RDTA: CA E, Y A siendo: 9 L_ 
dt dt s a 
No" 3 (0, * Cp) 
j=1 
BE RFP: A EA á 
dL dL 
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